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ABSTRACT 


Hall  effect  semiconductor  devices  and  titanium  dioxide  diodes  have  been  found  f 

to  be  highly  resistant  to  transient  X-ray  pulses.  Transient  X-ray  radiation  I 

effects  on  air  surrounding  resistive  "lements  as  a  function  of  pressure  were  J 

measured  and  a  theory  developed  for  \  effects.  Continuous  X-ray  speetrums  t 

have  been  calculated  from  X-ray  transmission  data  in  the  energy  range  from  I 

180  kilovolts  to  600  kilovolts  and  this  method  has  proved  to  be  satisfactory  | 

for  outputs  of  both  the  flash  X-ray  and  the  DC  X-ray  machines  ih  this  energy- 
range.  In  addition,  basic  radiation  effects  measurements  were  made  on 
dielectric  materials  using  the  Hall  effect  and  magnetoresistance  measurements 
in  a  steady  state  gamma  flux.  Hall  effect  magnetoresistance  effects  were 
observed  in  polystyrene,  polyethylene,  and  other  dielectric  materials.  A 
Schering  Bridge  method  was  used  to  determine  the  changes  of  capacitance  of 
components  during  flash  X-ray  exposure.  (This  bridge  is  much  less  responsive 
to  conductance  changes  than  to  capacitance  changes.)  And  a  Monte  Carlo  code 
has  been  written  which  describes  time -dependent  photon  transport  within 
finite  pieces  of  material.  The  energy  of  the  ejected  electrons  is  recorded 
and  photon  histories  are  traced. 
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SECTION  I 


TRANSIENT  GAMMA  RADIATION  EFFECTS 
ON  RESISTIVE  AND  INSULATING  MATERIALS 
By  W.  J.  Byatt  and  Harold  Cates 

1.  Introduction 

When  a  pulse  of  photons  is  incident  on  materials,  electrons 
are  generated  within  and  around  the  sample  by  photoelectric, 
Compton,  and  pair-production  interactions.  The  electrons  then 
move  so  as  to  produce  transient  effects  which  may  be  observed 
as  changes  in  conductivity,  current,  voltage,  or  charge. 

Experimental  findings  of  the  effects  of  transient  gamma 
irradiation  of  materials  are  oftentimes  analyzed  from  a  some¬ 
what  empirical  point  of  view.  For  example,  if  a  transient 
current-time  curve  resulting  from  irradiation  of  a  sample  is 
given,  the  ’technique  used  by  many  investigators  in  analyzing 
such  data  goes  roughly  in  the  following  way:  The  system  response, 
exclusive  of  the  component  under  irradiation,  is  calculated  or 
determined  experimentally  and  subtracted  from  the  data.  The 
remaining  current  is  then  plotted  versus  time  and  the  curve  is 
drawn  to  fit  the  data.  Often,  it  develops  that  a  series  of 
exponentials  will  fit  the  data  and  the  results  are  then  inter¬ 
preted  in  terms  of  a  trapping  model  having  the  observed  time 
constants . 

The  following  report  has  a  number  of  purposes.  First,  a 
number  of  papers  on  transient  radiation  effects  within  resistive 
and  dielectric  materials  will  be  reviewed.  Then,  the  results  of 
some  experiments  which  have  been  performed  on  resistors  will  be 
presented.  A  set  of  relations  which  can  be  used  to  predict  the 
observed  time-dependent  behavior  of  transient  voltage,  current, 
conductivity  and  charge  behavior  within  resistive  and  insulating 
materials  will  be  developed,  and  how  the  relations  can  be  used 


^i.e.,  the  response  of  such  things  as  coaxial  cable  to  gamma 
pulses,  or  the  measuring  system  response  time. 


to  analyze  not  only  data  of  this  experiment  but  data  of 
other  workers  in  the  field  will  be  shown. 

Before  embarking  on  the  above  program,  it  is  convenient 
to  discuss  some  basic  ideas.  It  is  known,  for  example,  that 
the  conductivity  of  a  gas  which  is  slightly  ionized  is  given 
by 


a(t) 


ne(t)q2r 


(1) 


where  n  (t)  is  the  number  density  of  electrons  of  charge  q 

and  mass  m  ,  and  r  is  a  suitably  averaged  relaxation  of  recom 
€ 

bination  time.  Now,  under  certain  conditions,  an  expression 
of  the  form  of  (1)  holds  in  metals  and  semiconductors,  it  is 
possible  to  write,  for  an  n-type  semiconductor. 


(2) 


with  n  (t)  the  time  dependent  behavior  of  n-type  impurities  of 

© 

charge  q  and  effective  (scalar)  mass  mg.  For  p-type  materials, 

V*>  -  (3) 

where  p(t)  is  the  concentration  of  holes  of  charge  q  and  effec¬ 
tive  scalar  mass  mg.  The  total  conductivity  can  then  be 

expressed  as 

°  ■  <{«**„  +  »*h]  w 


with  the  mobilities  defined  by,  for  example,  m  *  qr/m*.  The 
existence  of  a  number  density,  n#( t),  of  electrons  gives  rise 
to  the  possibility  that  if  the  electrons  move,  a  charge  density 
proportional  to 


p(t)  -  -qn^(t) 


2 


(5) 


will  be  transported.  Finally#  a  current  density  will  flow  if 

a  charge  is  transported.  It  will  be  given  by  either 
nt 

b 

j(t)  *  -qne(t)v  (6) 

or  by  solution  of  the  equation 

!£  +  ’  •  (j)  -  o  (7) 

It  is,  therefore,  obvious  that  a  knowledge  of  the  time-dependent 
behavior  of  free  electrons  is  of  importance  in  determining  the 
transient  behavior  of  materials  under  gamma  irradiation. 

2.  Development  of  Equations 

Once  again  the  expression  for  the  conductivity  (l)  should 
be  considered;  it  is 

,  .  n  (t)q2r 

<r(  t  -  - 

me 

The  concentration  of  free  electrons#  n  (t)#  after  a  gamma  ray 
pulse  has  passed  over  the  sample  under  test,  will  decrease  due 
to  recombination  of  the  electrons  with,  for  example,  ions  and 
traps.  One  can  write  that 

n0(t  +  At)  -  ne(t)  -  (At)  ne(t)  (8) 

which  states  that  the  number  density  (or  concentration)  of  free 
electrons  decreases  in  time  if  there  is  no  source  for  generation 
of  electrons.  The  factor  which  governs  the  decrease  is  called 
pabs  * n  the  relation  given  by  (8)#  and  is  tf*e  probability 

that  in  a  time  At  a  free  electron  will  disappear.  For  simplicity, 
it  is  assumed  that 

P.b.  *  Mt  (9) 


m- 


s 


>f 

Ls, 


c- 
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with  >,  a  constant.  Ori  passing  to  the  limit  At  0,  equation 
(8)  becomes 

dn 

3F“  =  "  ^  ne  (10) 

If  the  pulse  of  gamma  radiation  ceases  at  time  t  =  tQ,  at  which 
time  n  (t)  *  n  (t  ),  the  solution  to  (10)  is 

v  c  O 

ne(t)  “  ne(to)q~*t  ^  >  fco^  (xl) 

According  to  (l)  the  conductivity- time  curve,  after  the 
gamma  ray  intensity  has  gone  to  zero,  should,  if  assumption  (9) 
is  valid,  have  the  general  form 

a(t)  «  <r(t0)q“Xt  (t  £  tQ)  (12) 

It  is  of  interest  to  generalize  assumption  (9).  If  there 
are  several  time  constants  XT1  with  relative  weights  a.^,  a 
generalization  of  (12)  is 

-*,t 

<j(t)  -  or(tQ)  £  aie  1  (t  *  tQ)  (15) 

If  both  ai  and  depend  continuously  on  i,  then 

00 

•  0(0  \  a(i)«"A(i)tdi  (14) 

One  can  generalize  further.  If 

paba  *  x(t)  At*  0-5) 
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then  the  counterpart  of  equation  (ll)  is 

>(«)  <»« 

“,(*)  “  0  <l6> 

The  generalizations  of  equation  (13)#  (14)  and  (15)  then 
follow  naturally. 

One  could  just  as  well  have  been  deriving  equations 
for  the  current  density  or  the  charge  density  according  to 
equations  (5)  and  (6)  which  are  equivalent.  What  is  important 
is  a  knowledge  of  n  (t),  the  electron  concentration.  During 
this  discussion  it  is  assumed  that  the  gamma  pulse  is  square 
as  shown  in  Figure  1. 

The  above  relations  hold  only  for  t  *  t  #  with  t  the  time 

o  o 

at  which  the  gamma  pulse  amplitude  goes  to  zero.  It  is  easy 
to  remedy  this  shortcoming  by  including  the  effect  of  the  genera¬ 
tion  of  electrons  during  the  time  that  the  pulse  sweeps  over 
the  component.  Restricting  attention  to  the  time  behavior  of 

n  (t)#  one  writes 
© 

ne(t  +  At)  -  ne(t)  -  n*bs  +  n|en  (17) 

I 

The  last  two  terms  in  (17)  are  the  number  of  electrons  absorbed 
and  the  number  generated  per  unit  volume  in  a  time  At.  For  the 
number  absorbed  in  a  time  At#  it  is  postulated  that 

nebS  ■  »«(*>  eab*(Afc>  (18) 

stating  that  the  number  absorbed  in  the  time  interval  (t#  t  +  At) 
is  equal  to  the  number  present  at  time  t  multiplied  by  the  pro¬ 
bability  that  absorption  will  occur  in  time  At.  The  quantity 
pafcs(At)  is  proportional  to  the  time  interval  At,  with  the  con¬ 
stant  of  proportionality  being  the  concentration  of  absorbing  ' 


I 
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canters  of  various  kinds  and  their  cross  sections.  For  the 
latter 

K  -  (19) 

Then  equation  (18)  becomes 

n*b*  -  ne(t)  K  At  (20) 

The  number  of  electrons  generated  in  a  time  At  is  a  function 
of  the  time -dependent  incident  gamma  ray  intensity,  the  number 
of  electrons  available  for  generation,  and  the  cross  section 
for  the  interaction  of  a  gamma  ray  with  a  bound  electron.  The 
above  means  that 

n®®"  -  U  f[r  (t)]  At,  (21) 

where  the  quantity  n  is  proportional  to  the  product  of  the 
number  of  electrons  available  for  generation  with  the  cross 
section  for  the  interaction. 

When  (20)  and  (21)  are  substituted  into  (17)  and  the 
limit  At  -*  0  is  taken,  the  result  is 

+  K  ne(t)  -  Mf[r(t)]  (22) 

Equation  (22)  can  be  integrated  between  two  instants  of  time, 
tQ  and  t.  If,  at  tQ0  nc(t)  ■  the  solution  to  (22) 

is 

na(t)  -  n#(tQ)e  *  +  V  $  *[y(0>'K(t"T)  dr  (25) 

*o 

It  will  be  observed  that  if  tQ  is  the  time  at  which  the  gamma 
pulse  amplitude  goes  to  sero,  then  the  second  term  of  (25) 
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goes  to  zero  and  (2?)  is  identical  with  (ll).  The  generaliza¬ 
tions  following  equation  (ll)  then  hold  for  equation  (2j5)*  The 
expressions  for  the  current  density,  charge  density,  and  conduc¬ 
tivity  become 


j(t) 


-Kj  (t-t  )  f  r  -K.(t-r) 

j(t0)e  3  +  I ij  \  f j[^(T)]e  j  dT 


fco 


(24) 


(  r  .  “K«(t-r) 


q(t)  «  q(tQ)e  q  °  -q  ^  -q 


♦  Wq  ^  fq[y(r)]e  q  dr  (25) 


and 


cr(t) 


-K_(t-0  f  r  -k  (t-r) 

♦ "» i  ^(4, 


dr  (26) 


If  the  zero  reference  level  of  current  density,  charge  density 
and  conductivity  is  taken  just  prior  to  radiation,  and  if  t  ■  0# 
then  the  build  up  of  the  number  of  free  electrons  as  the  pulse 
passes  over  is  given  by 


f  r  -K(t-r) 

ne(t)  -  ^  dT 


0  <  t  «  t, 


(27) 


For  a  square  radiation  pulse  of  duration  t^  seconds,  such  that 
*[r(€)]  -  A  in  0  *  {  * 
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At  t  «  t1# 

-Kt 

n.fV  -  HS  (1  -  e  l).  (28) 

If,  in  (23),  tQ  is  now  set  equal  to  t^, 

-k( t  -  t.) 

ne(t)  -  ne(tx)e  1  t  *  (29) 

On  the  basis  of  the  very  simple  model  chosen.  Figure  1  gives, 
schematically,  the  time -dependent  behavior  of  the  number  of 
free  electrons.  The  . assumed  pulse  shape  is  also  shown.  From 
equations  (24),  (25)  and  (26),  the  behavior  of  the  current 
and  charge  densities,  together  with  the  conductivity  can  be 
found  if  equations  such  as(l)#  (5)  and  (6)  are  used  to  deter¬ 
mine  the  constants  of  proportionality  between  n  (t)  and  the 
quantity  being  computed.  The  equations  developed  in  this 
section  will  be  used  to  analyze  the  data  of  this  experiment 
as  well  as  the  data  of  other  workers  in  the  field. 


Figure  1.  Illustrating  the  Behavior  of  the  Number 
of  Free  Electrons  per  Unit  Volume  Under  the 
Influence  of  an  External  Rectangular  Pulse 
of  Duration  t^  Seconds 
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5*  Review  of  Previous  Work  on  Resistive  and  Insulating  Materials. 

Several  papers  will  be  reviewed  in  this  section.  The 
literature  in  transient  radiation  effects  is  now  so  enormous  that 
a  word  about  the  papers  reviewed  is  in  order.  In  part  4  of  this 
section  the  results  of  the  experiments  will  be  discussed.  The 
papers  reviewed  there  discuss,  in  general,  experimental  results 
which  may  be  compared  to  the  present  experiments.  This  fact  has 
influenced  the  selection  of  papers. 

Selections: 

(a)  Transient  Radiation  Effects  in  Electronic  Materials 
(Reference  1) 

The  object  of  the  above  report  written  by  D.  B.  Ebeoglu 
et  al.  is  to  determine  the  mechanisms  of  pulsed  radiation  effects 
in  electronic  materials.  The  sources  of  radiation  are  a  KEWB 
water  boiler  reactor  and  a  600  kev  flash  X-ray  unit.  Only  find¬ 
ings  from  the  latter  source  of  radiation  shall  be  reviewed. 

In  pulsed  radiation  studies,  air  ionization  effects  and 
cable  effects  give  rise  to  spurious  signals.  A  large  part  of 
the  report  is  devoted  to  a  discussion  of  methods  for  coping 
with  the  above  mentioned  effects.  It  is  to  be  noted  that  Ebeoglu 
et  al.  maintain  that  under  flash  X-ray  pulses,  the  use  of  epoxy 
coatings  is  not  so  effective  in  reducing  air  ionization  effects 
as  is  the  use  of  a  vacuum.  On  irradiating  polyethylene  at 
various  pressures,  it  was  observed  that  the  air  ionization  effect 
increased  until  the  pressure  within  the  test  chamber  was  about 
200  m  end  then  it  began  to  decrease.  The  decrease  continued 
until  a  pressure  of  about  10  m  was  reached,  after  which  there 
was  no  noticeable  change  (See  Table  IV  of  Reference  l). 

The  magnitudes  of  radiation-induced  signals  are  in  the  ratio  of 
12:8  5*5  for  atmospheric,  epoxy-coated,  and  vacuum  cases  respec¬ 
tively. 

The  width  of  the  pulse  of  X-rays  is  about  .1  M  second, 
while  the  decay  of  induced  conductivity  was  of  the  order  of  1 
u  second.  The  author  asserts  that  it  is  reasonably  certain  that 
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the  observed  decay  time  is  indeed  characteristic  of  the  poly¬ 
ethylene  since  care  was  taken  to  reduce  the  time  constant  of 
the  measuring  circuits  to  a  magnitude  significantly  less  than 
1  m  second.  It  was  necessary  to  eliminate  all  cable  conductor 
to  shield  capacitance  in  the  measuring  circuits.  This  was  done 
by  using  the  cables  as  wires  with  the  center  conductor  carrying 
the  signal  and  the  shield  being  left  unconnected.  With  such 
an  arrangement,  resistors  of  values  up  to  100  Kft  were  used  as 
test  components.  Above  such  a  value,  the  circuit  time  constant 
became  large  enough  to  dominate. 

The  authors  conclude  that  the  experiments  performed  on 
polyethylene  with  both  the  KEWB  reactor  and  the  flash  X-ray 
unit  point  to  the  excitation  of  trapped  carriers. 

In  the  above  report,  large  numbers  of  photographs  are  pre¬ 
sented  of  the  time  dependent  behavior  of  samples  under  irradia¬ 
tion.  These  are  displayed  together  with  the  irradiation  pulse 
shapes.  No  effort  was  made  to  fit  the  curves.  The  general 
form  of  the  observed  signals  is,  however,  as  shown  in  Figure  1. 
of  this  report. 

(b)  Some ftEf facts  of  Pulsed  Radiation  of  Electronic  Com- 
ponents-III.  (Reference  2) 

The  above  report  presents  the  results  of  tests  conducted 
at  the  Godiva  Reactor  at  Los  Alamos  in  January,  1959*  Only  the 
results  found  on  exposing  resistors  to  the  reactor  pulse  will 
be  discussed  herein.  With  the  test  circuit  used  by  the  IBM 
group,  the  change  in  the  resistance  under  test  is  given  by 


AR 


EAI 

I ( I+AI ) 


(30) 


where 

AR  »  the  resistance  change, 

E  -  the  applied  voltage 

I  -  the  original  current  through  the  test  resistor, 
AI  ■  the  change  in  current  through  the  test  resistor. 


It  is  pointed  out  that  the  resistance  varies  inversely  as  the 
flux  during  the  reactor  pulse  and  is,  therefore,  time  dependent. 
The  authors  present  AR^eaJc  in  tabular  form.  They  also  point  out 
that  changes  in  resistance  can  be  attributed  not  to  the  resis¬ 
tor  itself,  but  to  the  creation  of  a  shunt  leakage  path  about 
the  resistor  due  to  air  ionization.  In  Table  I,  Table  B1  of 
the  IBM  report  under  discussion  is  reproduced.  From  Table  I 
the  percentage  change  in  resistance  as  a  function  of  dose  rate 
can  be  plotted,  with  the  result  shown  in  Figure  2.  The  lower 
values  of  resistance  have  not  been  included  since  the  AR/R 
values  are  less  than  1 %  in  general.  From  the  data  of  Figure  2, 
it  is  clear  that  the  ionization  shunt  leakage  path  becomes 
increasingly  important  as  the  rated  value  of  the  resistor  under 
test  increases.  It  will,  therefore,  be  noted  that  if  circuits 
can  be  designed  having  low  values  of  pure  resistive  elements, 
at  least  the  changes  in  resistance  values  will  be  small;  i.e., 
the  circuits  will  be  more  radiation  resistant. 


?  rads  (tissue)/sec 


Figure  2.  Percentage  Change  in  Resistance 
as  a  Function  of  Dose  Rate 
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Table  I  Observed  Change  in  Resistance  at  Maximum 
Flux  (Taken  from  IBM  Rept.  59-816-916 f  Jan.#  1959) 


Resistance  Type 


Allen  Bradley 
Carbon  Compo¬ 
sition  2 -Watt 


Original 
Resistance 
( Ohms ) 


1.0553  x  10° 

1.0222  x  106 
1.0508  x  106 
1.0465  x  105 

1.0079  x  105 
1.0017  x  105 
1.0228  x  104 


AR 

( Ohms ) 


1.0240  x  10'  8.89  x  10C 


7.90  x  KK 
4.92  x  105 
5.47  x  105 
2.59  x  104 

I.985  x  10* 

1.55  x  104 
2.14  x  102 


Maximum  Gamma  Ray 
Dose  Rate  (Rads 
(Tissue )/sec) 

”  15  xlO6 

11  x  106 

6.4  x  106 

6.2  x  106 


10  x  10v 


5  x  10c 
6.5  x  10* 
6.4  x  10* 


1.0107  x  xO4  4.89  x 

1.0234  x  104  2.66  x 

1.0236  x  105 

9.9882  x  102  3.54 

1.0291  x  10? 

1.0244  x  105  8.99 

1.0228  x  105  12.63 

1.0215  x  102  .103 

1.0359  x  102  .024 

1.0170  x  102  .053 


102  10  x  106 

102  8.1  x  106 

8.8  x  106 

4.5  x  106 
16.0  x  106 

8.1  x  106 
14  x  106 
14  x  106 

6.5  x  106 

4.5  x  106 
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(c)  Transient  Radiation  Effects  on  Electronic  Parts 
(Reference  3) 

The  main  purpose  of  the  above  mentioned  report  is  to  pre¬ 
sent  experimental  data  on  radiation- induced  effects  in  coaxial 
cable.  In  particular,  the  radiation  induced  conductivity  in 
RG-58  cable  having  a  solid  polyethylene  dielectric  is  found  to 
vary  approximately  as  the  .55  power  of  the  dose  rate  in  the 
range  of  dose  rates  from  .36  to  7  x  10  rads/sec.  For  a  given 
dose  rate,  the  conductivity  is  independent  of  voltage  applied 
to  the  cable  in  the  range  -300  -  V  £  300  volts.  Secondary  elec¬ 
tron  emission  is  constant  over  a  measured  dose  rate  of  from  .37 
o 

to  10  rads/sec.  The  conclusion  is  that  gamma  photoconductivity 
is  important  at  low  dose  rates,  while  secondary  emission  is 
important  at  high  dose  rates . 

(d)  Transient  Radiation  Effects  on  Electronic  Parts 
(Reference  4) 

A  relatively  comp’ete  discussion  of  transient  radiation 
effects  in  a  Mylar  capacitor  is  given.  The  analysis  is  presented 
in  terms  of  the  equivalent  circuit  shown  in  Figure  3- 


Figure  3*  Equivalent  Circuit  for  a 
Capacitor  under  Irradiation 
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The  quantity  R  (t)  is  a  time-dependent  shunt  resistance  induced 
hy  the  incident  radiation.  The  quantity  R  is  a  series  resis¬ 
tance  representing  the  parallel  combination  of  1  megohm  resistor 
in  series  with  the  power  supply  and  the  1  megohm  input  impedance 
of  the  measuring  amplifier.  The  induced  conductivity  across  the 
capacitor  C  is  shown  schematically  in  Figure  4  together  with  the 
radiation  pulse. 


Figure  4.  Induced  Conductivity-Time 
Curve  in  a  Capacitor  due  to  the  Radia 
tion  Pulse  Shown. 


Rather  than  following  the  analysis  of  the  authors  of  the  report 
under  discussion  Figure  4  will  be  discussed  on  the  basis  of  an 
equation  developed  in  part  2  of  this  report.  Consider,  there¬ 
fore,  equation  (26)  for  the  conductivity.  It  is 

t 

-KL(t~t  )  ( 

tr(t)  -  °  *  uff  \ 
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In  order  for  (26)  to  conform  with  the  data  of  Figure  4,  it 
is  convenient  to  set  t  =  -  It  I,  and  to  take  the  range  of  t 

O  ^  /  V  / 

to  be  -|t  |  ^  t  >  ®.  Experimentally  it  is  found  that  a(to)  =  a(- 
P  .  y  * 

=  or  ,  an  initial  condition.  Further,  since  y(r)  exists  only 

in  the  interval  -  | t  |  ^  t  ^  0,  equation  (26)  becomes 


"p0 


-K  (t+|t  |) 


ff(t)  =  CT^e  CT'  P  +  fCT[y(T)]e  a  dT(-jt. 


t  *  o)  (31) 


-it 


and 


a(t)  = 


cr  e 
P 


-K  (t  +  1 1  |  ) 
or'  1  p  1  * 


t  Sr  0 


(32) 


The  most  primitive  assumption  regarding  fa|_'>'(T) 


is  that  it  is 


constant  in  - | t  |  ^  t  ^  0.  With  such  an  assumption, 

P 


or(t)  =  a  e  CT 
'  P 


-K_(t+|t  |  nAr  -K  (t+|t  |), 


'P‘+  ^1-e  "O' 


K  L" 
CT 


-|tp|  S  t  S  o  (75) 


and 


-K  (t  +  |t  j ) 

or(t)  =  CT  e  a  P 

P 


t  o 


(34) 


In  the  above  relation,  A  is  the  constant  amplitude  of  fCTj_y(T) 

in  - 1 1  |  ^  t  rs  0. 

P 

It  will  be  observed  that  as  t  =  0  is  approached  from  the 
negative  side, 


-Kit  I 
cr(t)  -  ct  e  P 

'  '  p 


U  A  -K  It 

Ml  "  e  a  P  ) 


K 


(35) 
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while  as  t  -*  0  from  the  positive  side. 


o(t) 


-K„|t  | 


V 


(36) 


From  these  two  forms,  it  is  seen  that  the  observed  discontinuity 
in  a(t)  at  t  =  0  is  predicted.  While  it  would  be  of  some 
interest  to  find  the  values  of  the  parameters  K^,  and 
which  sould  permit  a  fit  of  the  data  in  the  report  under  discus¬ 
sion,  it  is  not  possible  since  no  absolute  magnitudes  are  given, 
(e) 

1)  "Gamma-Ray  Photoconductivity  Decay  in  Organic 
Dielectric  Materials."  (Reference  5) 

2)  "Measured  Behavior  of  Gamma-Ray  Photoconductivity 
in  Organic  Dielectrics."  (Reference  6) 

3)  "Gamma-Ray  and  Neutron- Induced  Conductivity  in 
Insulating  Materials."  (Reference  7) 

The  above  three  papers,  by  staff  members  of  the  Sandia 
Corporation  of  Albuquerque,  New  Mexico,  represent  some  of  the 
more  definitive  work  done  in  the  area  of  transient  radiation 
effects.  The  papers  discuss  the  transient  gamma  ray  response 
of  a  variety  of  materials,  including  polyethylene,  polypropylene 
polyisolretylene.  Mylar,  Teflon,  cellulose  acetate,  reconsti¬ 
tuted  mica,  tantalum  oxide,  and  polyvinylchloride,  in  all  of 
the  quoted  materials  the  authors  find  that  the  response  of  the 
materials  to  gamma  radiation  can  be  split  into  three  regions 
which  they  call  A,  B,  and  C.  The  three  regions  are  shown  in 
Figure  5*  It  will,  first  of  all,  be  noted  that  the  logarithm 
of  the  conductivity  is  plotted  versus  time  in  Figure  5*  In 
region  A,  Harrison  et  al.  find  that  the  conductivity- time 
curve  varies  as 

-t/t 

q  -  o0  -  A(1  -  .  °)  (37) 


'4. 
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}  where  crQ  is  the  dark  conductivity  and  tq  is  a  function  of  the 
|  dose  rate.  Specifically,  it  is  found  that  rQ  can  be  determined 
■  from  a  knowledge  of  dose  rate  by  writing 

?  To  -  VU  (38) 

t 

|  where  kq  and  m  are  determined  empirically. 

•  In  region  B  equilibrium  is  reached,  at  which  time  genera- 

j  tion  of  Carriers  is  in  equilibrium  with  recombination.  In 
f  region  B, 

!  O  -  o0  -  A  *c  (39) 

| 

f  with  A.,  and  C  constants  to  be  determined  empirically. 

y 

?  In  region  C  it  is  found  that 

n  -t/r. 

<r  “  o  S  K.e  1  (40) 

eg  iml  i 

with 


"aq  ~  °o*  V'  (41) 


Figure  5*  The  Conductivity-Time  Behavior 
for  the  Gamma  Pulse  Shown. 
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It  will  be  observed  that  equation  (37)  is  of  the  form  of 
equation  (28)#  while  equation  (40)  is  identical  with  (13)* 

The  fact  that  region  B  exists  in  the  experiments  of  Harrison 
et  al,  is  due  to  the  time  duration  of  the  pulse .  it  is  of 
the  order,  in  many  of  their  experiments,  of  10  seconds.  This 
means  that  the  generation  of  carriers  by  th&  external  source 
can  come  to  equilibrium  with  the  recombination  of  carriers 
within  the  material.  In  the  experiments  which  we  shall  dis- 

_7 

cuss  in  part  4,  the  pulse  duration  is  of  the  order  of  2  x  10  ' 
second,  and  equilibrium  cannot  be  established. 

4.  Experimental  Work 

A  Fexitron  flash  X-ray  unit  capable  of  delivering  a  .2 
M  sec.  pulse  at  accelerating  voltages  up  to  600  kev  provided 
the  gamma  ray  source  in  the  experiments  to  be  described.  The 
gamma  ray  pulse  is  collimated  by  a  1  inch  hole  cut  in  about 
1/2  an  inch  of  lead.  A  Consolidated  Vacuum  Corporation  vacuum 
unit  with  a  glass  bell  jar  of  average  wall  thickness  of  about 
1/4  inch  is  placed  inside  a  screen  room.  The  bottom  plate  of 
the  vacuum  chamber  has  openings  for  a  number  of  electrical 
leads.  Lead  shielded  coaxial  cable  runs  from  connections  made 
outside  the  bell  jar  through  2  inches  of  lead  block  shielding 
to  a  Tektrona:rll21  A  amplifier.  The  outer  shields  of  the 
coaxial  cable  tsed  in  the  test  circuit  are  grounded.  A  coaxial 
cable  runs  frowv  the  output  of  the  1121  A  amplifier  to  a  551 
Tektronix  oxcil^frscope  with  a  type  CA  plug-in  unit.  The  cable 
mentioned  above  provides  a  matched  impedance  from  the  output 
of  the  1121  A  to  the  oscilloscope  input.  Figure  6  is  a 
schematic  diagram  of  the  circuit  used  for  testing  the  compon¬ 
ents.  Figure  7  ii  a  drawing  of  the  vacuum  system.  The  1121 
A  amplifier  and  the  remainder  of  the  circuitry  as  shown  in 
Figure  6  are  enclosed  in  a  lead  brick  shield  to  the  right  of 
the  vacuum  system. 

Certain  tests  made  to  determine  whether  there  was  any 
discernible  pickup  as  the  flash  X-ray  unit  was  pulsed  when  no 
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test  component  was  in  place  across  the  two  contact  posts  shown 
in  Figure  7  are  mentioned  before  describing  the  experiments 
themselves.  An  18  volt  potential  difference#  supplied  by  the 
power  supply  indicated  in  Figure  6  was  placed  across  the  bind¬ 
ing  posts.  If  air  ionization  effects  were  important#  if 
secondary  electron  emission  from  the  glass  walla  a  pro¬ 

blem#  or  if  spurious  signals  were  generated  within  the  cables, 
some  effect  would  be  observed.  None  was  seen.  When  a  resistor 
was  placed  between  the  binding  posts  and  a  2  inch  lend  block 
was  placed  inside  the  bell  jar  so  as  to  interrupt  the  photon 
beam#  no  effect  was  seen.  When  the  lead  block  was  removed#  there 
was  a  voltage  pulse  observed  on  the  oscilloscope.  Finally#  when 
a  bare  wire  was  inserted  between  the  binding  posts  and  the  x-ray 
unit  was  pulsed#  there  was  no  detectable  signal  at  the  Invelm 
at  which  we  were  operating.  On  the  basis  of  the  above  ramsults# 
it  can  be  assumed  that  spurious  signals  due  to  secondary  elec¬ 
tron  emission  from  the  bell  jar  wall,  cable  effects#  and  the 
like  are  not  a  factor  in  these  experiments. 

One  final  comment  should  be  made  before  the  results  to  date 
are  briefly  described.  Far  more  thorough  analyses  of  the  data 
will  be  given  in  the  theses  of  H.  T.  Cates  and  L.  C.  Meyer  of 
the  University  of  New  Mexico. 

Resistors 

Allen  Bradley  carbon  composition#  solid  core  resistors  have 
been  tested  with  rated  values  of  1  megohm#  1  watt#  mod  1  megohm 
1/4  watt.  The  tests  have  been  run  with  the  flash  X-ray  tube 
operating  at  4 00  kev  with  an  anode  to  test  sample  distance  of  8 
1/4  inches.  At  this  distance#  the  dose  rate  is  approximately 
10^  roentgens  per  second.  Since  the  pulse  duration  is  .2  x  10~^ 
second#  a  total  dose  of  approximately  .2  roentgen  was  delivered 
by  the  tube  at  the  position  of  the  sample.  Clearly#  the  .2 
roentgen  is  not  absorbed  dose.  The  vacuum  system  was  operated 
such  that  exposures  to  the  X-ray  pulse  took  place  at  pressures 
between  atmospheric  and  10  microns. 
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A  variety  of  exposures  was  taken.  In  some,  both  the  leads 
and  the  resistor  itself  were  not  potted.  In  some  experiments 
only  one  lead  was  potted,  and  in  others  both  leads  were  potted. 
Finally,  not  only  were  the  leads  potted,  but  the  body  of  the 
resistor  was  also  potted. 

The  experimental  findings  are  presented  in  a  series  of 

figures.  In  Figure  8  the  magnitude  of  the  peak  voltage  developed 

across  the  10  K  viewing  resistor  for  1  megohm  resistors  of  1/4 

watt,  1  watt,  and  2  watt  power  ratings  is  plotted.  From  the 

figure,  it  is  seen  that  the  peak  amplitude  of  the  voltage  is 

inversely  proportional  to  the  pressure  from  atmospheric  to  about 

150  mm.  As  the  pressure  is  further  lowered,  the  amplitude  of 

the  voltage  finally  becomes  independent  of  pressure  changes,  in 

previously  reported  tests  of  the  kind  discussed  here,  the  same 

general  shape  of  curve  was  seen.  The  detailed  behavior  of  the 

peak  voltage  versus  pressure  as  discussed  by  Ebeoglu  et  al. 

(Reference  1)  is  somewhat  different.  The  fact  that  the  present 

tests  were  at  about  5000  ft.  (the  altitude  of  Albuquerque) 

equivalent  pressure  and  at  low  humidity  can  account  for  the 

difference.  In  Figure  9  there  is  plotted  as  a  function  of  R , 

0 

the  original  resistance,  at  one  dose  rate.  It  is  seen  that  as 
the  value  of  RQ  is  increased,  the  percentage  change  in  the  resis¬ 
tance  increases.  This  result  can  be  explained  from  the  equation 
for  the  time -dependent  resistance 


R(t) 


R0Rs(t) 
Ro  + 


(42) 


where  R  (t)  is  the  shunt  air  leakage  resistance.  Thus  a  greater 
8 

percentage  change  in  total  resistance  would  be  expected  for  a 
higher  value  of  resistance.  This  result  is  in  keeping  with  the 
IBM  results  shown  in  Table  I. 

In  Figure  10,  the  value  of  AV  as  a  function  of  the  tube 
accelerating  potential  (or  dose  rate  when  the  conversion  is  made) 
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rOl  *  9 


Acceleration  Potential  (kev) 


Figure  10.  Maximum  Transient  Voltage 
as  a  Function  of  Accelerating 
Potential 


Reproduced  From 
2  5  Best  Available  Copy 


is  given  for  a  1  ohm,  1  watt  resistor.  The  results  from  Figure 
10  have  been  plotted  again  in  Figure  11  as  AR  versus  dose  rate. 
As  previously  reported  in  the  TREES  Handbook  (Reference  8)  a 
model  for  radiation  effects  observed  in  resistors  postulates  a 
shunt  leakage  path  given  by 

».<*>  “  FTtJ 

and  a  replacement  current  equal  to 


(44) 


where  ?(t)  is  the  dose  rate  in  r/sec#  A  is  a  constant  having 
units  of  ohm-r/sec,  and  B  is  a  constant  having  units  of  amp- 
sec/r. 

The  values  for  A  and  B  calculated  from  these  experiments 
at  atmospheric  pressure  (630  mm  Hg)  for  an  unpotted  Allen 
Bradley  1  met,  1  watt,  carbon  composition  resistor  are  the 
following! 


A  -  2.1  x  10 
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B  - 


7-8  x  IQ"”  amp-sec/r. 


Figure  11.  Maximum  Change  in  Resistance 
as  a  Function  of  boss  Rates 
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These  values  are  In  the  range  of  those  reported  in  the  TREES 
Handbook.  Figure  12  shows  the  relative  effects  of  an  unpotted 
resistor  at  atmospheric  pressure,  an  unpotted  resistor  at  a  low 
pressure,  and  an  epoxy-coated  resistor  at  atmospheric  pressure. 
Ebeoglu  et  al.  maintain  that  evacuating  the  test  chamber  reduces 
the  transient  effect  on  resistors  more  than  potting  them.  The 
results  of  these  experiments  are  contrary  to  their  findings. 
However,  it  has  been  found  that  when  different  thicknesses  of 
epoxy  coating  on  the  resistor  are  irradiated,  the  one  with  the 
heavier  coating  gave  the  smaller  transient  effect.  It  is 
possible  that  Ebeoglu  et  al.  had  thin  epoxy  coatings.  A  com¬ 
parison  of  potting  the  resistor  leads  and  potting  the  resistor 
and  its  leads  was  made  and  there  was  not  any  noticeable  difference 
in  the  transient  pulse. 

For  1  roentgen  per  hour  source  of  photons  of  energy  of  the 
order  of  .5  mev  the  gamma  ray  energy  flux  is  about  5  *  lO3  mev/cm 
sec.  The  number  flux  of  photons  is  twice  this:  namely,  N-,  = 

10  mev/cm  -sec.  The  dose  rate  delivered  by  the  flash  X-ray  unit 
is  about  10^  roentgens  per  second.  Therefore,  for  the  source, 

Ny  *  5600  x  1012  y's/cm2-sec  or 

Ny  3.6  x  10«  y* s/cm2-sec  (4 6) 

This  number  flux  of  photons#  when  multiplied  by  the  total  Klein - 
Nishina  cross  section  for  the  generation  of  electrons,  will 
yield  the  probability  per  unit  time  that  a  single  electron  is 
generated.  On  denoting  the  probability  density  by  p(t)«  the 
order  of  magnitude  estimate  is  obtained: 

p(t)  *  x  lO1*  x  .6  x  10"24  sec-1 

■  2.16  x  10"9  sec"1  (47) 

Since  the  gamma  pulse  duration  is  of  the  order  of  .2  x  10"^ 
second  the  total  probability  that  an  electron  is  generated  is 
given  approximately  by 

j 9  -  p(t)dt  •  .4  x  10-W  (48) 
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Then  the  number  of  electrons  generated  per  cubic  centimeter 
is  the  number  available  times  the  probability  of  generation. 
At  atmospheric  pressure  the  number  of  electrons  per  cubic 
centimeter  is  of  the  order  of 

N  M  4  x  1020  cnf5  (49) 

G 

The  number  of  free  electrons  generated  in  the  gas  for  the 
pulse  length  will  be 

n£  -  pN  m  105cm"5  (50) 

G  G 


The  conductivity  of  an  ionized  gas  can  be  calculated  frcm 

! 

J  a  * 

:)S 

.  l 

i 

4 

I  M  *  mass  ion 

I 

r  *  relaxation  time 

J 

when  '  ’mai.  At  atmospheric  pressure  the  dominant  process 
|  leading  to  the  recombination  of  free  electrons  is  the  attachment 
i  of  the  electrons  to  oxygen  molecules.  The  time  constant  for 
I  the  process  is  given  by 

l 

!  r"1  -  (52) 

u2 

where  X  ~  2  x  10“^°cm^/sec  and  Nq  is  the  square  of  the  number 
density  of  air  molecules.  The  2  calculated  value  of  r  is  of 
the  order  of 


2  f 

q  »e  T 

M 


(51) 


T  « 


.2  x  10“7 


seconds . 


(55) 
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Now  c r  can  be  calculated.  Ite  magnitude,  with  Mf  «  lO^cm"^, 

2  —20  2  —7  ® 

<3  «  25  x  10  (esu)  *  r  »  .2  x  10"'  second,  and 

m  «  9.H  x  10-27  gma  becomes,  on  converting  from  (statohm-cms)" 
“  -1 

to  (ohm-cms)"  ,  is 

a  «  10"^  (ohm-cms) (54) 

Since  the  shunt  resistance  will  be  inversely  proportional  to 
the  conductivity,  one  should  expect  to  find  a  peak  value  of 
R_  (the  shunt  resistance)  of  the  order  of  10^  ohms  if  t/k0  the 
length  to  area  ratio,  is  of  the  order  of  unity.  Calculations 
confirm  that  the  value  of  the  shunt  resistance  in  these  experi¬ 
ments  is,  indeed,  of  the  order  of  1  megohm. 

The  peak  voltage  developed  across  the  10  KO  viewing  resis- 
tor  in  these  experiments  is  about  5  x  10  volts.  The  current 

flowing  through  the  viewing  resistor  is,  therefore,  of  the 
2 

order  of  a  cm  ,  the  current  density  will  be  of  the  order  of 

j  «  5  x  10"7  amperes/cm2  (55) 

It  is  of  interest  to  compute  the  expected  velocity  of  the  elec- 
troti*  within  th*  gas  if  it  is  assumed  that  a  current  density  of 
the  same  ;  of  magnitude  as  that  given  above  flows  from  the 

gas  across  the  boundaries  of  the  irradiated  resistor.  Thus, 

j(nj  q)_1  -  v  (56) 

with  j  *»  5  *  10-7  ampere'  'cT.i  nf  w  lo^cro-^,  and  q  *  1.6  x  10“^ 

© 

amps -sec 


7 

v  m  10'  cm/ sec. 


(57) 


Now  this  velocity  is  what  one  would  expect  if  the  electrons  had 
come  into  thermal  equilibrium  with  molecules  within  the  gas.  It 
is  certainly  not  typical  of  electron  drift  velocities  within  a 
conductor  or  resistor.  These  are  of  the  order  of  a  few  cm/sec. 

We  may  interpret  the  findings  as  providing  further  evidence  for 
the  fact  that  the  transient  effects  observed  in  irradiating 
resistors  are  wholly  due  to  effects  within  the  gaseous  medium 
surrounding  the  sample  under  test. 

This  discussion  on  experimental  findings  will  be  concluded 
by  examining  the  behavior  of  the  time  decay  of  the  voltage  pulse 
observed  on  irradiating  resistors.  If  the  attachment  of  elec¬ 
trons  to  oxygen  molecules  is  the  dominant  process  over  most  of 
the  pressure  range  in  which  we  have  conducted  experiments,  it 
would  be  possible  to  correlate  the  theoretical  predictions  of 
time  constants  with  the  experimentally  observed  ones .  Figure 
13  is  a  plot  of  the  theoretical  time  constants  found  by  evaluating 


with 


x  n; 


6.023  x  10 
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2.24  x  10 


ip  («2)  §- 


(58) 

(59) 


The  number  of  molecules  in  22.4  liters  of  air  is  just  6.023  x 
2”5 

10  .  The  factor  .2  accounts  for  the  fact  that  only  one  fifth 

of  the  total  number  (approximately)  are  oxygen  molecules. 

Finally,  since  the  number  will  decrease  with  pressure,  the  fac¬ 
tor  p/pQ  is  the  ratio  of  pressure  p  to  atmospheric  pressure  p0. 

Superimposed  on  the  theoretical  curve  of  Figure  13  are 
experimental  points.  The  agreement  between  theory  and  experiment 
in  the  range  of  pressures  between  1  and  10  nans  of  Hg  ir.  within 
the  uncertainties  of  plotting  the  oscilloscope  traces.  Above 
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tf-aa  »»*•• 


m 

•I 

if 

ft 

|  100  mms  of  Hg  time  constants  wars  found  to  be  of  the  order 
I  of  2  x  10‘ 6  seconds  to  pressures  of  630  mms  of  Hg.  It  is 
|  felt  that  such  a  time  constant  is  characteristic  of  these 
§  circuits  and  not  of  the  decay  processes  in  the  gas.  Attempts 
|  were  made  to  lower  the  time  constant  of  the  circuit  of  Figure 

h 

I  6  by  replacing  the  10  Kfl  viewing  resistor  with  a  1  KO  viewing 
J  resistor.  Severe  noise  problems  resulted,  and,  as  a  conse- 

af 

|  quence,  the  construction  of  a  screen  room  about  the  capacitor 
|  bank  of  the  flash  X-ray  unit  has  been  recommended.  Below  pres- 
|  sures  of  10  mms  of  Hg  the  time  constant  of  the  decay  begins  to 
1  drop  towards  values  of  a  few  microseconds  -  once  again.  Since 
|  radiative  recombination  may  be  the  dominant  process  at  such  low 
V  pressures,  and  since  the  constant  for  such  a  process  is  of  the 

it  _7 

I  order  of  10  1  sec,  it  is  expected  that  a  drop  in  the  value  of 
f  the  time  constant  should  occur  and  that  it  should  be  masked  by 

I  the  circuit  time  constant. 

$ 

'  5.  Conclusions 

Experiments  have  been  performed  which  involved  irradiating 
carbon  composition  resistors  with  the  Kirtland  Air  Force  Base 
flash  X-ray  unit.  From  these  experiments  the  following  informa¬ 
tion  has  been  found x 


l)  The  magnitude  of  the  transient  voltage  observed  on  an 
oscilloscope  is  a  function  of  pressure.  The  voltage  increases 
until  a  pressure  of  about  130  mms  of  Hg  is  reached,  and  then  it 
decreases  until  a  pressure  of  1  mm  of  Hg  is  reached.  The  mag¬ 
nitude  then  remains  constant. 


2)  The  time  constants  of  the  voltage  decay  agree,  within 
experimental  error,  with  a  model  which  maintains  that  the  dom¬ 
inant  process  in  a  gas  loading  to  disappearance  of  free  elec¬ 
trons  is  the  attachment  of  free  electrons  to  oxygen  molecules 


over  the  pressure  range  from  10  to  100  mms  Hg.  Above  100  mms 


of  Hg,  time  constants  of  voltage  decay  are  of  the  order  of  a  few 
microseconds,  these  being  characteristic  of  the  measuring  circuit. 
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3)  A  theory  had  been  developed  In  which  the  time -dependent 
behavior  of  the  number  of  free  electrons  plays  a  dominant  role. 
The  theory  has  been  used  to  calculate  the  peak  conductivity, 
current  density,  and  pulse  shapes  of  the  present  work,  as  well 
as  that  of  other  investigators.  The  agreement  is  good. 

4)  It  has  been  found  that  the  value  of  resistance  changes 
as  a  function  of  dose  rate  and  the  percentage  change  in  rated 
resistance  value  depends  critically  on  the  original  value  of  the 
resistance.  Such  work  confirms  earlier  findings  which  are  also 
discussed  in  this  report. 

5)  Potting  of  resistors  leads  to  a  decrease  of  air  ioniza¬ 
tion  effects,  and  it  has  been  found  that  potting  of  the  leads 
only  has  a  major  effect  in  reducing  the  magnitude  of  the  observed 
transients.  This  suggests  that  heavier  coatings  cause  a  further 
reduction. 

6)  One  may  conclude,  as  has  been  done  earlier  by  other 
workers,  that  transient  resistance  changes  under  gamma  irradia¬ 
tion  are  entirely  due  to  air  ionization  effects.  No  bulk  effect 
was  observed,  in  order  to  minimize  air  ionization  effects,  it 

is  important  that  circuits  be  designed  having  as  low  an  impedance 
as  possible  for  the  job  they  are  designed  to  do. 


SECTION  II 


TRANSIENT  RADIATION  EFFECTS 
ON  HALL  DEVICES 

By  W.  W.  Grannemann,  Harold  Cates  and  LeRoy  Meyer 


1.  Introduction 

A  complete  discussion  of  the  transient  radiation  effects 
on  Hall  devices  is  given  in  AFWL  TDR-64-38,  June  1964.  Any 
application  of  the  Hall  device  in  a  nacletar  environment  would 
require  a  knowledge  of  the  transient  radiation  effect  on  them. 
Therefore,  experiments  and  calculations  were  performed  to 
determine  the  gamma  radiation  effects  on  both  thin  film  and  bulk 
type  devices.  The  results  of  this  work  will  make  it  possible  to 
separate  the  gamma  radiation  effects  from  the  EM  and  circuit 
effects.  A  few  applications  of  Hall  devices  include  the  following: 
magnetic  flux  meters,  compasses,  magnetometers,  clip-on  DC 
ammeters,  computor  elements,  modulators,  amplifiers,  variable 
attenuators,  frequency  multiplexing,  and  DC  to  AC  conversion. 

The  commercial  devices  used  were  the  Ohio  Semiconductor 
Halltrons,  and  the  Beckman  th in-film  Hallefex  voltage  generators. 
These  are  the  indium  antimonide  (Halltron  HS-51  and  Hallefex  551) 
and/or  indium  arsenide  (Halltron  BS-51  and  Hallefex  350)  devices. 
Also  some  inSb  Hall  devices  were  made  in  the  University  of  New 
Mexico  laboratory  with  special  non-inductive  leads. 

The  gamma  radiation  sources  that  were  used  are  the  AFWL 
600  KV  flash  X-ray  and  the  White  Sands  gamma  Linac. 

2.  Theory 

a.  Gamma  and  Electron  Irradiation  of  Hall  Device  Materials 
The  effects  discussed  here  apply  primarily  to  InSb,  but 
much  of  the  same  phenomena  apply  to  other  semiconductor  materials. 

The  effects  studied  are  primarily  those  caused  by  gamma 
radiation,  although  electrons  have  been  used  with  Linac  experi¬ 
ments  . 
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Uhder  the  heading  of  TREES  (transient  radiation  effects 
on  electronic  systems),  there  are  several  effects  that  may  be 
important,  depending  on  the  environment  and  instrumentation. 

These  are  ionizing  effects,  sudden  creation  of  permanent  defects, 
and  the  sudden  rise  of  temperature.  In  many  cases  these  effects 
are  small  and  there  is  difficulty  in  separating  them  from  cable 
effects  and  noises.  While  these  small  effects  are  unimportant  as 
far  as  transient  radiation  effects  on  electronic  circuits  are  con¬ 
cerned,  it  is  desirable  to  understand  them  and  see  them  clearly 
enough  in  experiments  so  that  predictions  can  be  made  as  to  their 
effects  in  higher  dose  rate  environments. 

The  excess  electrons  and  holes  temporarily  created  in  a 
semiconductor  by  exposure  to  ionizing  radiation  enhance  the  con¬ 
ductivity  and  Hall  constant  of  the  media.  The  conductivity  and 
Hall  constant  are  a  function  of  the  majority  carriers,  in  addi¬ 
tion  to  this  current  in  the  volume  of  the  device,  there  may  be 
an  additional  current  component  from  surface  conduction,  or  con¬ 
duction  to  some  othetr  elements  of  an  encapsulating  or  supporting 
structure.  The  degree  to  which  these  excess  carriers  affect  the 
operation  of  the  device  will  depend  on  the  function  of  the  device. 
An  estimate  of  the  magnitude  of  the  radiation-grnerated  current 
can  be  obtained  for  a  device  by  knowing  the  geometry  and  electri¬ 
cal  characteristics. 

At  normal  temperatures  the  generation  and  recombination 
processes  constantly  occur  due  to  the  influence  of  thermal  energy. 
A  state  of  equilibrium  is  maintained  by  these  processes  as  deter¬ 
mined  by  Fermi -Dirac  statistics.  This  distribution  is  independent 
of  the  manner  in  which  carriers  are  either  generated  or  recombined. 

The  change  in  carrier  density  in  a  semiconductor,  6n,  can  be 
estimated  from  the  known  radiation  dose  rate,  the  efficiency  of 
the  material  in  producing  carriers  (.56  to  .72  ev/pair  for  Insb), 
and  the  carrier  recombination  rate  for  the  material.  If  the 
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radiation  pulse  is  short  compared  to  carrier  lifetime*  r,  the 
total  number  of  excess  carriers  at  the  end  is  proportional  to 
the  radiation  dose.  If  the  carrier  lifetime  is  much  shorter 
than  the  radiation  pulse,  the  excess  carrier  density  (fin)  at  any 
time  during  the  pulse  is 

fin  *  finfct  (60) 


where  t  is  the  time  duration  of  the  radiation  pulse*  and  fin^  is 
obtained  from  the  continuity  equation  as  follows: 


r 


time  rate  of 
increase  of 
conduction 
electrons 

fin. 


rate  of 
generation 
\  of  conduction 
electrons 


irate  of  re¬ 
combination 
of  electron- 
lion  pairs 


(61) 


The  generation  rate  g  can  be  estimated  if  the  energy  of  the 
radiation  beam  and  the  absorption  coefficient  for  the  material 
are  known.  The  energy  absorbed  is  given  by  the  expression 


Eab.  (8V) 


cm 


a  gram 


cm 


ram  .  10 'ev 

57155  1.9xl0-$.rg. 


(62) 


where 

y  ■  absorption  coefficient 

ol 

IQ  ■  energy  of  gamma  beam 

The  energy  lost  per  electron  ionized  in  any  material  can  be 
approximated  as  2  to  4  times  the  ionization  potential  of  the 
material  ( Reference  9 ) .  This  ionization  potential  in  solid-state 
materials  is  interpreted  to  be  the  forbidden  energy  gap*  which 
is  0.18  ev  for  inSb.  The  ionizing  efficiency  for  InSb  is  then 
about  G.36  to  0.72  ev  per  electron-ion  pair.  The  processes 
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involved  are  the  Compton  scattering  events  and  photoelectric 
events.  All  the  calculations  were  made  using  the  numbers  which 
gave  the  maximum  radiation  effect.  In  this  case  .36  is  used 
and  the  number  of  carriers  produced  is  then 


n_ 


energy  absorded  (ev) 
ev/electron-ion  pair 


(63) 


If  a  square  wave  radiation  pulse  of  duration  t,  is  assumed, 
the  generation  rate  g  of  conduction  electrons  is  given  by 


*  -  tr 


(64) 


An  equation  similar  to  equation  (6l)  can  be  applied  to 
holes.  The  excess  conductivity  can  be  shown  to  be 


6cr  =  q(nn6n  +  Mp6p) 


(65) 


where 


Un  *  electron  mobility 
M  =  hole  mobility 

tr 

6  =  excess  number  of  holes 

P 


Equation  (65)  reduces  to 


6 a  =  6nqu 


(66) 


M  »  M 

n  p 


for  n-type  material  where  n  »  p  and 

For  steady  state,  R  =  g,  and  the  steady-state  concentra¬ 
tion  increment  is 


6n  =  gr 


(67) 


where 

r  =  mean  lifetime  of  the  carriers. 
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Equation  (67)  correspbnds  to  equation  (3)#  page  F-2,  of 
the  TREES  Handbook  where  it  is  given  in  the  form 

n  -  nQ  =  gr  (68) 

where 

n  -  n  =  6_  in  this  report, 
on 

The  equation  given  in  the  TREES  Handbook  for  g  is 

g  =  Kf  (69) 


where 

K  =  the  energy-dependent-,generation  constant  ^electron-hole 
pairs  per  erg  g"l(c)J 

?  -  the  gamma  exposure  rate  (ergs  g~*(c)  sec”*) 

Since  the  constant  K  is  not  known  for  this  material  and 
energy,  the  above  derivation  was  needed  to  arrive  at  a  genera¬ 
tion  rate  g. 

These  empirical  equations  look  nice  mathematically,  but 
are  not  always  applicable  unless  all  the  constants  are  known. 

r  *  50  x  10"9  sec  for  InSb,  and  thus  equation  (67)  would 
apply  for  InSb  when  calculating  the  change  in  carrier  density. 

For  build-up, 

fin  *  gr[l  -  exp  ( -  ^)]  (70) 

And  for  decay, 

fin  «  gt  exp( -  |)  (71) 


d 

It 


i 


where 

t  ■  time  in  seconds 
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b„  Effects  on  Hall  Voltage  Output 

Since  current  is  related  to  voltage,  the  increase  in 
the  number  of  carriers  calculated  by  equation  (6?)  will  cause 
a  change  in  conductivity  of  the  material  6a  by  the  following 
relation 


6a  =  q6nMn  (72) 

Since  current  is  related  to  the  voltage,  V,  across  the  device 
by  the  relation 

i  *  crvwd/t  (75) 


where 

w  <=  width  of  device 
d  =  depth 
l  =  length 

the  change  in  current  through  the  device  will  then  be 

6i  =  6aVwd/-t  (7*0 


The  change  in  Hall  voltage  output,  may  be  calculated  for 

three  different  modes  of  operation,  if  it  is  assumed  that  the 
conventional  Hall  voltage  equation  holds  under  the  radiation 
conditions.  This  relationship  is  given  by 


(75) 


where 

*  Hall  coefficient  (cnrVcoul) 
B  -  Magnetic  field  (grams) 

I  «  Current  lines  (amperes) 
d  *=>  Material  thickness  (cm) 


The  three  cases  considered  are  for  the  Hall  device  operated 
with  the  bias  circuit  as 

1.  a  constant  voltage  source 

2.  a  constant  current  source 

5.  a  constant  power  source 

(l)  Constant  Voltage  Source 

This  circuit  arrangement  is  shown  in  Figure  14. 


Figure  14.  Constant  Voltage  Source  Bias  Circuit 


By  using  the  relations  for  the  Hall  coefficient  and  current 
density  in  the  Hall  voltage  equation,  one  gets  the  following 
results : 

»h  -  k  (76) 

where 

n  ■  number  of  conduction  electrons 
q  *  charge  of  an  electron 
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and 

I  _  qn  MVS 
wcf  “  Z 

I  *  =  qn  n  ^  Vs  (77) 

where 

M  =  mobility 

Vg  =  source  voltage 

d,w,  and  l  =  dimensions  of  sample 


Then 


V1 

”3 


(78) 
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s 


(79) 


Equation  (79)  indicates  that  for  a  constant  voltage  source 

the  Hall  voltage  should  remain  constant  during  a  radiation 

pulse,  since  B,  m,  V  and  the  device's  dimensions  are  not  a 

s 

function  of  the  radiation  energy  absorbed. 

(2)  Constant  Current  Source 

The  constant  current  source  has  a  large  output 
impedance  compared  to  the  Hall  device  input  impedance.  This 
circuit  arrangement  (Figure  15)  keeps  the  current  through  the 
device  constant.  Then  using  equations  (75)  and  (76)  one  may 
obtain 


(80) 


The  current  in  the  Hall  voltage  equation  will  be  a  constant. 
Then  if  n  changes  by  an  amount  6 n,  a  quantity  6V^,  the  change 
in  Hall  voltage  is  defined  by 


vh +  svh 


nTln 


(81) 
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where 


K  =  BI 
K1  qd 


so  that 


6V, 


K. 


K, 


n 


(82) 


'h  ”  n  +  6n 

(3)  Constant  Power  Source 

If  energy  is  supplied  to  the  Hall  device  at  a  constant 
rate,  the  result  that  Hall  voltage  is  a  function  of  is  obtained 


in  the  following  manner.  The  power  equation  is 


Jn 


P  =  I  R 


(83) 


The  resistance,  R,  can  be  expressed  in  terms  of  the  resistivity 
and  dimensions  of  the  device  by 

>1 


Then 


wd 


x  ■  [£*] 


1/2 


(84) 


(85) 


Next,  using  equation  (84)  and  the  following  relations  for  R^ 
and  p  in  the  Hall  voltage  equation  one  may  obtain 


*h  -  UP. 


and  p  = 


nqM 


V. 


= '  £)1/2 
vn' 

k2 
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when  K2  *  E(|^)1/2 
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Figure  15 .  Constant  Current  Source  Bias  Circuit 
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If  n  changes  by  an  amount  5n,  the  change  in  Hall  voltage,  6v^ , 
will  be  given  by 


+  5vh 


so  that 


(89) 


Equations  (82)  and  (90)  indicate  that  the  Hall  voltage  should 
decrease  when  n-type  material  is  irradiated  by  gamma  rays. 

c.  Theoretical  Results 

In  the  AFWL  TDR-64-J8  the  constant  current  case  is 
taken  and  the  change  in  Hall  voltage  calculated,  for  a  typical 
operating  condition,  always  assuming  the  worst  case,  for  various 
radiation  dose  rates.  The  results  were  that  there  was  no  signi¬ 
ficant  change  in  Hall  voltage  for  dose  rates  up  to  10^  r/sec. 

An  n-type  Indium  Antimonide  bulk  type  device  was  used  in  the 
calculations . 

At  a  dose  rate  of  greater  them  10^  r/sec  both  holes  and 
electrons  had  to  be  taken  into  account  (which  involves  a  much 
more  complicated  expression),  and  a  transient  radiation  effect 
was  predicted. 

The  thin  film  devices  have  even  less  calculated  transient 
radiation  effect  than  the  bulk  type  device  because  they  have 
less  material  in  which  to  abosrb  energy.  This  does  not  eliminate 
the  possibility  of  having  a  transient  radiation  effect  due  to 
lead  arrangements  or  other  inherent  error  effects  in  Hall  devices. 
There  are  a  number  of  causes  of  Hall  voltage  errors  which  would 
have  to  be  investigated  for  a  particular  application.  These 
are:  1)  loops  id  lead  arrangements  which  are  necessary  in 
making  contact  to  the  material,  but  could  cause  errors  when 
measuring  changing  magnetic  fields;  2)  displacement  voltage 


which  is  due  to  misalignment  of  Hall  voltage  leads?  j)  r-f 
pick-up?  4)  air  ionization?  and/or  5)  cable  effects. 

3.  Experiment 

In  the  60C  KV  flash  X-ray  experiments  no  transient  radiation 
effects  were  observed  when  the  constant  current  circuit  configura¬ 
tion  was  ue^d.  Also  there  was  no  change  in  Hall  voltage  when  the 
constant  '  »ge  circuit  was  used.  The  dose  rates  were  10^  r/sec. 

It  wa.  und  that  the  bulk  type  device  had  an  upper  cut-off 
frequency  w*.  round  150  KC.  To  be  sure  that  the  transient  radia¬ 
tion  effect  is  not  being  lost  due  to  band  width  limitations,  the 
White  Sands  <  imma  Linac  was  used.  It  had  a  gamma  pulse  width  of 
10  microseconds.  No  change  in  Hall  voltage  was  observed  for  dose 
rates  up  to  10^  r/sec  using  either  constant  current  or  constant 
voltage  circuits. 

4.  Conclusions 

The  Hall  devices  were  found,  by  experiments  with  the  flash 
X-ray  and  gamma  Linac,  to  be  radiation  resistant  for  dose  rates 
up  to  10^  r/sec.  Theoretical  calculations  indicate  a  radiation 

Q 

resistance  up  to  at  least  l(r  r/sec.  These  results  were  obtained 

using  the  constant  current  circuit  configuration.  Theoretical 

calculations  for  the  constant  voltage  circuits  indicate  no 

change  in  Hall  voltage  for  n-type  material.  At  dose  rates  above 
o 

10^  both  holes  and  electrons  become  involved  and  the  Hall  voltage 
equation  becomes  more  complicated.  A  transient  radiation  effect 
is  predicted  at  the  higher  dose  rates .  Higher  dose  rate  machines 
are  needed  to  check  the  theory  at  the  higher  dose  rates. 
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SECTION  III 


THE  FABRICATION  AND  IRRADIATION  OF 
TITANIUM  DIOXIDE  DIODES 
By  James  Woodward  and  W.  W.  Grannemann 

1 .  Introduction 

Probably  the  biggest  need  in  the  transient  radiation  effects 
program  is  for  a  highly  radiation  resistant  diode  and  transistor 
which  exhibit  both  desirable  electrical  characteristics  and  resis¬ 
tance  to  transient  radiation  and  steady  state  radiation.  Titanium 
dioxide  diodes  have  some  interesting  possibilities  in  this  respect. 
Recent  research  at  the  University  of  New  Mexico  indicates  that 
titanium  dioxide  diodes  are  highly  radiation  resistant  to  perma¬ 
nent  type  damage.  The  study  in  this  section  was  instigated  to 
provide  some  information  on  the  transient  radiation  properties 
of  titanium  dioxide  diodes. 

Titanium  dioxide  diodes  are  not  presently  commercially  avail¬ 
able.  So  this  section  covers  the  fabrication  of  the  devices,  a 
mathematical  approximation  of  the  transient  radiation  effects,  and 
the  results  of  irradiating  the  devices  with  the  flash  X-ray. 

2 .  Fabrication  of  Titanium  Dioxide  Diodes 

Titanium  dioxide  diodes  can  be  fabricated  by  several  dif¬ 
ferent  methods  in  the  laboratory.  Three  of  the  most  practical 
and  successful  methods  are  listed  below: 

1)  Titanium  Dioxide  on  Titanium  Base  Using  Forced  Steam 

Dne-Step  Process  (Reference  10) . 

Heated  titanium  will  readily  coat  with  a  semiconducting 
titanium  dioxide  film  when  steam  is  passed  over  the  sample.  As 
the  oxide  is  formed  some  reduction  takes  place  by  the  surplus 
hydrogen,  causing  oxygen  vacancies  in  the  oxide  film.  The  final 
result  is  an  n-type  semiconductor  film  which  gives  diode  action 
to  impressed  voltages. 

2)  Titanium  Dioxide  on  Titanium  Base  Using  a  Two-Step 
Heating  Process  (Reference  11) . 

In  this  method  an  oxide  film  is  formed  by  heating  the 
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titanium  sample  in  oxygen.  This  is  followed  by  reduction  in  a 
hydrogen  atmosphere  creating  the  n-type  titanium  dioxide  film 
on  the  titanium  sample. 

3)  Rutile  Point  Contact  Diode  (Reference  12). 

In  this  method  of  manufacture  crystal  rutile  is  heated  in 
a  hydrogen  atmosphere  causing  oxygen  vacancies  in  a  surface  layer 
of  n-type  semiconductor.  The  sample  is  cleaved  to  expose  a  clean 
surface  which  is  secured  to  a  copper  stud,  and  a  point  contact  is 
applied  to  the  titanium  dioxide  semiconductor, 
a .  Diode  Manufacture 

It  has  been  determined  that  the  one-step,  steam  reduction 
process  gives  the  most  reliable  and  repeatable  results.  This  is 
most  probably  due  to  the  fact  that  reduction  of  the  titanium 
dioxide  film  takes  place  continuously  during  the  film  coating 
process,  thereby  providing  a  nearly  homogeneous  distribution  of 
oxygen  vacancies.  The  other  two  processes  are  basically  surface 
processes;  hence  the  reduction  can  be  considered  to  be  more  con¬ 
centrated  near  the  surface.  Diodes  made  by  these  methods  are 
not  consistent  in  characteristics  and  are  more  prone  to  aging 
(Reference  13 )•  Because  of  the  aging  factor,  the  titanium  dioxide 
diodes  were  fabricated  using  the  one-step  steam  process. 

Chemically  pure  samples  of  titanium  .5  x  5  x  .016  inch  in 
size  were  prepared  from  sheet  Ti-55A  (this  is  chemically  pure 
titanium).  A  copper  lead  was  attached  to  one  .5  x  .3  inch  surface 
by  brazing  with  a  helium  shielded  arc.  Attaching  the  copper  ohmic 
lead  to  the  pure  metal  in  this  manner  prior  to  coating  the  sample 
with  the  oxide  film  was  considered  mors  advantageous  than  attempting 
this  step  later  on  in  the  manufacturing  process.  The  alternative 
was  to  completely  coat  the  sample  with  an  oxide  film  and  then 
scrape  away  part  of  the  film  so  that  an  ohmic  contact  could  be 
made  to  the  pure  titanium.  This  alternative  procedure  risks  heat 
and  mechanical  damage  to  the  semiconductor  film;  therefore,  it 
was  not  attempted. 
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The  pure  titanium  samples  with  copper  leads  attached  were 
then  polished  to  provide  a  smooth  surface  on  which  to  coat  the 
oxide  film.  This  is  considered  an  important  step  as  rough  sur¬ 
faces  were  found  to  readily  cause  shorting  through  the  semicon¬ 
ductor  film. 

Oxide  films  were  coated  on  the  samples  by  placing  them  within 
a  1  1/2  inch  quartz  tube  which  had  been  inserted  in  an  electric 
furnace.  A  constant  flow  of  superheated  steam  was  passed  through 
the  quartz  tube  and  over  the  samples  at  all  times  while  the 
samples  were  at  elevated  temperatures.  The  equipment  used  required 
approximately  8  to  10  minutes  to  bring  the  samples  to  temperatures 
in  the  vicinity  of  700 °C. 

A  variety  of  temperatures  and  time  periods  for  depositing 
the  film  was  investigated.  Temperatures  below  650°C  for  periods 
up  to  1  hour,  in  general,  produced  oxide  films  which  readily 
shorted  out  under  applied  voltages.  Satisfactory  results  were 
obtained  by  keeping  the  samples  in  an  atmosphere  of  750°C  for 
approximately  30  minutes  while  passing  superheated  steam  over 
them.  It  is,  however,  estimated  that  longer  exposures  of  about 
3  hours  would  have  produced  satisfactory  results  at  temperatures 
in  the  650°C  range . 

Oxide  films  obtained  in  this  manner  had  a  black -blue -gray 
coloring  which  could  not  be  readily  scraped  off.  The  film  is 
estimated  to  be  10  cm  thick. 

Leads  to  be  applied  to  the  semiconductor  surface  were  pre¬ 
pared  from  thin  copper  sheet.  The  strips  of  sheet  copper  were 
bent  so  that  an  area  of  about  1/8  x  1/1 6  inch  was  placed  flat 
against  the  oxide  film;  the  strips  were  fastened  with  conducting 
epoxy.  The  conducting  cement  required  a  20  minute  temperature 
cure  at  200°F.  This  did  not  cause  any  nbticeable  change  in  the 
diode  characteristics. 

b.  Typical  Characteristics  of  the  Manufactured  Diodes 

Titanium  dioxide  diodes  fabricated  by  the  above  method  exhibit 
characteristics  which  compare  favorably  with  the  older  types  of 
germanium  diodes  (Reference  14) .  The  manufactured  diodes  had 
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consistently  similar  characteristic  curves.  The  voltage  versus 
current  for  a  typical  diode  is  shown  in  Figure  16.  When  voltages 
were  first  applied  across  the  diodes  most  of  them  showed  an  ini¬ 
tial  characteristic  curve  of  less  slope  than  that  shown  in 
Figure  16.  After  5  or  10  seconds  the  slope  tended  to  increase 
and  stabilize,  never  again  returning  to  the  curve  initially 
observed.  This  is  believed  to  be  the  "healing  effect"  frequently 
discussed  in  literature  about  diode  characteristics.  Also  obser¬ 
vable  was  a  slight  hysteresis  effect  when  voltages  were  rapidly 
swept  across  the  diode  (Reference  15)-  The  reason  for  this  phe¬ 
nomenon  was  not  immediately  clear,  but  further  investigation  of 
this  characteristic  was  not  attempted, 
c.  Additional  Preparation  for  Testing  of  the  Samples 

The  completed  diodes  were  cast  in  a  resin  sealant  to  pro¬ 
tect  the  diode  from  mechanical  deunage  during  subsequent  handling. 
Without  potting  the  diode  in  this  manner  it  was  difficult  to 
attach  cable  leads  for  further  testing  without  twisting  off  the 
semiconductor  surface  contact  lead  as  the  sample  was  moved  about 
or  stored.  The  resin  sealant  also  served  to  isolate  the  semi¬ 
conductor  film  from  the  atmosphere  and  reduced  the  possibility 
of  aging  or  contamination  from  atmospheric  impurities  (Reference 
16). 

Additional  resin  sealant  was  poured  around  the  diode  after 
test  cable  leads  were  attached.  The  quantity  of  resin  sealant 
poured  was  large  enough  to  effectively  reduce  or  eliminate  ioni¬ 
zation  effects  across  the  diode  leads  when  the  diode  was  irra¬ 
diated.  The  sealant  completely  enclosed  the  diode  as  well  as 
the  soldered  connections  to  standard  shielded  cable  used  at 
irradiation  testing  facilities. 

3.  Computation  of  Radiation  Effects  Upon  Titanium  Dioxide  Diodes 

It  must  from  the  beginning  be  pointed  out  than  any  computa¬ 
tion  of  a  transient  current  in  a  semiconductor  material  can  only 
be  as  accurate  as  the  physical  constants  and  crystal  parameters 
used.  It  is  beyond  the  scope  of  this  report  to  determine  the 
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Figure  16.  Typical  Characteristic  (Currant 
Versus  Voltage)  Curve  of  Fabricated 
Titanium  Dioxide  Diodes 


51 


i*  -4iiR' 


exact  values  of  the  ionization  energy,  electron  mobility,  etc. 
of  the  fabricated  crystal.  Instead,  reasonable  values  for  these 
and  other  constants  have  been  selected  from  the  literature  to 
give  as  near  as  possible  the  correct  order  of  magnitude.  Other 
assumptions  as  to  the  nature  of  the  crystal  doping,  the  hole 
current,  the  diffusion  current,  etc.  also  detract  from  the  accu¬ 
racy  of  these  calculations. 

At  best,  the  computations  which  follow  can  be  considered  a 
crude  approximation  of  what  actually  takes  place  in  the  Tic^ 
semiconductor  crystal  when  exposed  to  a  pulse  of  radiation.  The 
general  method  as  well  as  the  computations  may,  however,  serve 
as  a  starting  point  for  the  laboratory  investigation  of  radiation 
effects . 

The  Air  Force  Weapons  Laboratory  flash  X-ray  facility  at 

• 

Kirtland  Air  Force  Base  was  used  to  radiate  the  manufactured 

titanium  dioxide  diodes.  Computation  of  radiation  effects  upon 

the  diode  must  of  necessity  start  with  consideration  of  the  range 

and  limitation  of  the  flash  X-ray  system.  X-ray  tube  location, 

shielding,  and  instrumentation  considerations  dictated  that  the 

sample  to  be  radiated  should  be  placed  6  inches  from  the  X-ray 

tube.  Data  available  at  the  facility  for  dose  and  dose  rate  at 

the  6  inch  location  are  shown  in  Table  II  for  the  300  KV  to 

2 

500  KV  range.  The  dose,  designated  I  ,  is  given  in  erg/cm  ,  and 

°  2 

the  dose  rate,  Io/sec,  is  given  in  erg/cm  sec. 

After  consideration  of  tube  life  and  waiting  periods  between 
pulses  (7-10  minutes  for  cooling  after  a  400  KV  burst)  it  was 
decided  that  400  KV  operation  would  best  satisfy  the  test  require¬ 
ments.  Computations  to  follow,  therefore,  are  based  upon  a  dose 
(lQ)  of  .9  x  lO-1*  erg/cm^. 

Brewster  lists  values  for  the  mass  absorption  coefficient 
for  titanium  dioxide  crystal  (rutile)  (Reference  10) .  The  coeffi¬ 
cients,  designated  W,  associated  with  various  wavelengths  are 
shown  in  Table  III. 
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Table  II.  Dose  and  Dose  Rate  at 
6  Inches  From  X-Ray  Tube 


Zo  2 
erg/cm 

I  /sec 

°  2 
erg/ cm 

X  lO-1* 

A 

-  sec  x  10 

300 

KV 

0.60 

0.30 

400 

KV 

0.90 

0.45 

500 

KV 

1.80 

0.90 

Table  III.  Mass  Absorption  Coefficients  of  Rutile 


Wavelength  (X) 

in  Angstrom  Units  .01  0.1  1  2 

w  in  cm2/ gm  .0543  .1790  38.71  250.4 


X-ray  wavelength  (x)  is  given  by  the  expression  {p-  where 
h  is  Planck's  constant,  c  the  speed  of  light,  and  E  the  energy. 
For  400  KV, 


X  *  .031  A‘ 


Graphical  interpolation  on  log  paper  of  the  values  in  the  above 
chart  gives  W  *  .075  cm2/g  for  .031  A°. 

The  total  absorbed  energy,  in  the  titanium  dioxide  semi¬ 

conductor  is  given  by  the  expression 


E«b.  ■  WIo*  (9!) 

where  w  is  the  mass  absorption  coefficient,  IQ  the  dose,  and 
g  the  mass  of  the  semiconductor.  The  thickness  of  the  semi- 
conductor  is  taken  as  10  x  cm  and  the  total  area  of  the  serai- 

2  "x 

conductor  surface  as  3-22  cm  ,  giving  a  volume  of  3-22  x  10  citr 
and  a  mass  of  13.8  x  10  grams.  (The  density  of  rutile  is  4.28 
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(92) 


ll  ,11)  .V^a,W,»V..  ,  r:  - - ~:r.T-  -.::: 


g/cnrr^)  (Reference  17) . 

Babs  *  5.85  x  107  ev 

For  gamma  irradiations  of  this  frequency  range  almost  all 
of  the  absorbed  energy  of  the  incident  gamma  rays  upon  the  semi¬ 
conductor  material  is  converted  into  free  electron  production 
and  motion  by  photoelectric  effect  and  Compton  events.  The 
energy  lost  per  ionization  in  any  material  can  be  approximated 
as  2  to  4  times  the  ionization  potential  of  the  material.  This 
approximation  is  commonly  found  in  the  literature. 

This  ionization  potential  in  semiconductors  is  interpreted 
to  be  the  forbidden  energy  gap  which  is  3 .05  ev  for  rutile  TK^ 
(Reference  17).  This  value  can  be  taken  to  be  a  close  approxi¬ 
mation  for  the  semiconductor  material  in  a  back  biased  diode. 

The  energy  lost  per  electron  produced  is  therefore  taken  to  be 
approximately  6.10  ev. 

The  number  of  electrons  (n)  produced  is  computed  by  the 
expression 

E  _ 

N  *  =  .957  x  107  electrons  (95) 

If  a  square  wave  radiation  pulse  of  .2  microsecond  duration 
is  assumed,  the  rate  g(t)  of  generation  of  conduction  electrons 
is  given  by 

g(t)  *  *  4.78  x  101*  electrons/sec  (94) 

where 

t  ■  time 

The  mobility  of  electrons  (u  )  in  reduced  rutile  at  room 

n  2 

temperature  is  believed  to  be  less  than  10  cm  /volt-sec. 

Current  flow  in  the  reduced  titanium  dioxide  semiconductor 
material  is  known  to  be  purely  electronic  (Reference  18).  X-ray 
irradiation  produces  additional  electrons  for  conduction.  Hole 
current  is  non-existent  or  at  best  negligible. 
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The  electron  mobility  (n  )  is  related  to  the  diffusion 
constant  by  the  Einstein  relationship  (Reference  19) . 

Mn  ”  StT  Dn 

where  at  500°K#  ~  *  .026  volt 
Therefore, 

o 

D  £  =  0.26  cm  /sec 

n 

The  transient  current  flow  through  the  back  biased  diode 
will  be  caused  by  the  excess  electrons  generated  by  the  inci¬ 
dent  radiation.  The  diffusion  current  occurs  whenever  there 
is  a  density  gradient  of  carriers  (^n).  This  current  is  related 
to  the  density  gradient  of  electrons  by  the  expression. 


\  "  q  Dn  7n 


(96) 


In  the  case  of  a  thin  film,  the  density  gradient  can  be 
considered  equal  to  n/l,  where  n  is  the  concentration  of  charge 
carriers  and  t  is  the  thickness  of  the  semiconductor  film.  Then, 

in  “  q  Dn  T  (97) 


The 

function 


concentration  of  charge  carriers  at  any  time  t  is 
of  both  the  generation  g(t)  and  the  current  flow. 


a 

so 


that 


n 


g(t)  -  \-S  dt 


number  of  electrons 
which  recombine 


(98) 


If  it  is  assumed  that  the  electrons  recombining  with  ions  are 
negligible,  the  in  computed  is  a  "maximized”  value  greater  than 
what  should  be  found  experimentally.  For  the  time  period  of  the 
radiation  pulse. 


qDn 

TTV5T) 


(99) 
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and  for  t  >  .2  x  10  sec,  when  electron  generation  has  ceased 

q  d _ (  v 

-  -  q 


[-  ^  at  ] 
t  (voi)  L  J  q 


(100) 


By  substituting  values. 


in  =  .618  t  -  8.07  •  102*  ^  in  at  (t  <  .2  x  10"6  sec) 

(101) 

and  solving  for  i^, 

in  *  7.65  •  10 "5  [l  -  exp(-8.07  X  105t)]  (t  <  .2  x  10"6  sec) 

(102) 

at  t  »  .2  x  10'  when  the  radiation  pulse  ends, 

in  *  7-65  *  10'5  [l  -  exp( - . 1614) ]  -  .000115  ma  (105) 


The  current  as  a  function  of  time  after  the  radiation  pulse  has 
passed  is 

in  s  1.13  •  10-7  exp(-8.07  x  105  t)  (t  >  .2  •  10"6  sec) 

(104) 

A  graph  of  the  transient  current  rise  and  fall  as  expressed 
by  the  two  equations  above  is  plotted  graphically  in  Figure  17 
along  with  the  radiation  square  pulse. 

The  computations  above  do  not  take  into  consideration  elec¬ 
tron-ion  recombination,  surface  effects,  and  trapping  effects.  A 
detailed  numerical  analysis  of  these  three  effects  must  start 
with  a  thorough  study  of  the  doped  semiconductor  crystal,  measure 
ment  and  calculation  of  the  various  crystal  parameters,  and 
selection  of  a  suitable  semiconductor  model. 

Ths  Schottky  model  (Reference  20)  appears  best  for  our  par¬ 
ticular  diode.  This  model  pictures  a  metal  in  contact  with  a 
uniform  semiconductor. 

The  fabrication  process  for  our  titanium  dioxide  diode, 
whereby  a  semiconductor  film  is  deposited  upon  pure  titanium. 
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produces  an  n-type  semiconductor.  Reduction  by  steam,  causing 
oxygen  vacancies,  takes  place  continuously  throughout  the  period 
that  the  film  is  deposited.  The  vacancies  are  assumed  to  be  uni¬ 
formly  distributed  throughout  the  crystal.  This  type  crystal 
lends  itself  to  the  Schottky  model. 

If  further  analysis  reveals  the  existence  of  some  sort  of 
surface  layer  between  the  titanium  metal  and  the  semiconductor, 
consideration  should  be  given  to  the  Bethe  model  (Reference  21) . 
This  model  takes  into  consideration  an  insulating  layer  of  a 
specified  thickness  between  the  metal  and  a  uniform  semiconductor. 

Consideration  of  the  electron-ion  recombination  phenomenon 
surface  effects,  and  the  trapping  phenomenon  as.  discussed  above 
would  modify  the  transient  pulse  so  far  computed  and  shown  in 
Figure  17 . 

Electron-ion  recombination  will  reduce  the  number  of  elec¬ 
trons  available  for  the  diffusion  current  and,  therefore,  will 
reduce  to  some  extent  the  peak  value  of  the  transient  pulse. 
Thermal  release  of  charge  carriers  from  the  trapping  centers 
will  cause  a  slower  decay  of  the  transient  current  peak. 

The  transient  current  pulse  shown  in  Figure  17/  therefore, 
displays  two  important  limits  for  the  expected  transient  pulse. 

The  peak  value  expected  to  be  found  experimentally  should  not 
exceed  the  peak  value  shown,  and  the  time  constant  for  the  rise 
and  fall  of  the  transient  pulse  should  not  be  less  than  that 
computed  and  shown  in  Figure  17 . 

4.  Laboratory  Investigation  of  Transient  Effects 
a .  Instrumentation 

Within  the  screen  room  are  Tektronix  type-551  oscillo¬ 
scopes  with  type  L  and  type  CA  plug-in  units  for  recording  test 
results.  These  also  have  Polaroid  camera  attachments.  Polaroid 
film  with  a  film  speed  of  10,000  ASA  is  normally  used. 

A  bank  of  power  supplies  is  available  for  instrumentation. 

It  includes  a  Lambda  regulated  power  supply  transistorized 
small  voltage/current  supply  and  various  dry  and  wet  cell  bat¬ 
teries  . 
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RG-58  C/U  cable  with  appropriate  connecting  and  matching 
hardware  is  normally  used  for  device  circuitry  and  instrumen¬ 
tation  . 

b.  Radiation  Test  Set  Up 

It  was  considered  desirable  to  establish  test  circuitry 
which  would  provide  information  about  both  the  voltage  across 
the  diode  and  the  current  flow  through  the  diode  during  the 
radiation  pulse  and  afterward  until  steady  state  conditions 
were  again  obtained.  This  meant  displaying  current  versus 
time  and  voltage  versus  time  on  the  scopes  where  the  pulse 
width  was  .2  x  10 second.  The  Tektronix  scopes ,  type  55 1# 
have  this  capability  with  the  results  permanently  obtained  for 
permanent  record  by  Polaroid  camera.  Carey  however,  must  be  taken 
to  provide  a  circuitry  which  is  accurately  responsive  during  time 
periods  of  such  short  duration. 

The  first  circuitry  investigated  included  a  circuit  loop 
with  diode,  power  supply,  and  current  sensing  resistor  in  series. 
Display  probes  from  the  scopes  were  connected  across  the  current 
viewing  resistor  and  the  diode  itself  (see  Figure  18) .  All  leads 
were  coaxial  cable  with  the  copper  shielding  grounded. 

It  was  soon  found  that  the  extensive  footage  of  cable  and 
the  interaction  of  two  display  connections  were  causing  oscilla¬ 
tions  or  "ringing"  which  obscured  any  valid  test  results. 

Efforts  then  were  made  to  reduce  the  cable  footage  and 
examine  only  the  transient  current  through  the  viewing  resistor. 
It  was  also  noted  that  excessive  noise  and  oscillation  occurred 
when  the  diodes  were  forward  biased  and  exposed  tc  radiation 
pulses.  For  this  reason,  and  the  fact  that  the  back  biased  con¬ 
dition  is  of  more  interest  in  radiation  effects  studies,  it  was 
decided  to  conduct  all  tests  at  a  back  bias  value  of  6  volts. 

The  final  configuration  used  consisted  of  the  basic  circuit 
loop  but  with  a  Tektronix  1121  amplifier  close  to  the  test  sample 
(see  Figure  IQ) .  This  arrangement  eliminated  the  long  cables, 
and  the  system  proved  to  have  adequate  time  response.  A  Ruther¬ 
ford  pulse  generator  was  used  to  thoroughly  check  out  the  system 
and  to  determine  that  the  frequency  response  of  the  circuit  was 
within  the  time  limits  of  our  measurements. 
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Figure  18.  Trial  Circuitry  for  Examination 
of  Transient  Current  and  Voltage 


Lead  Shielding 


Radiation 


Figure  19.  Final  Circuitry  for  Testing  of  Diodes 
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Additional  shielding  of  the  power  supply  and  the  Tektronix 
1121  amplifier  was  accomplished  by  using  lead  blocks  and  a  lead 
sheet.  Only  the  diode  itself  and  less  than  1  inch  of  the  con¬ 
necting  leads  were  exposed  to  the  direct  radiation. 

An  additional  check  of  the  test  circuitry  was  made  under 
radiation  conditions  by  comparing  the  results  obtained  from  the 
irradiated  diode  with  the  results  obtained  from  the  same  diode 
after  it  had  been  lowered  1  inch  out  of  the  direct  radiation 
beam.  In  all  cases  the  transient  radiation  effects  upon  the 
lowered  diode  were  either  undetectable  or  negligible  in  compari¬ 
son  with  the  effects  observed  on  the  diode  in  the  radiation  beam. 
It  was  concluded  that  all  results  by  the  irradiated  diodes  can 
be  attributed  only  to  the  diode  and  its  short  leads  exposed  to 
radiation.  Pickup  elsewhere  in  the  circuit  directly  or  indi¬ 
rectly  caused  by  the  radiation  pulse  is  therefore  considered 
negligible, 
c.  Potting  of  Diodes 

A  commercial  diode  with  a  characteristic  curve  similar  to 
the  fabricated  titanium  dioxide  diodes  was  irradiated  under 
various  environments.  When  it  was  irradiated  at  atmospheric 
pressure  without  any  potting,  a  maximum  transient  current  excur¬ 
sion  of  .010  ma  was  observed.  The  same  diode  in  a  partial  vacuum 
(10  microns)  had  a  maximum  excursion  of  .008  ma.  When  encased 
in  resin  epoxy  (3/4  inch  diameter  cylinder)  the  maximum  excursion 
was  reduced  to  .004  ma  (see  Figure  20).  It  is  believed  that  this 
variation  in  results  is  caused  by  ionization  around  the  diode 
leads  providing  a  secondary  path  for  current  around  the  diode 
during  the  radiation  burst. 

Some  potting  material  was  used  in  fabricating  the  titanium 
dioxide  diodes  and  could  not  be  removed  without  structurally 
damaging  the  diode.  It,  therefore,  appeared  invalid  to  compare 
a  partially  potted  titanium  dioxide  diode  with  commercially 
unpotted  diodes.  It  was  determined  best  to  "pot"  all  diodes  in 
an  identical  3/4  inch  cylinder  of  potting  epoxy.  In  this  manner 
a  valid  comparison  could  be  achieved  and  at  the  same  time  ioniza¬ 
tion  effects  of  all  samples  would  be  reduced  under  similar  environ 
ments. 

62 


mmmm m 


d .  Noise  Levels 

The  Tektronix  oscilloscopes  used  to  display  the  transient 
current  have  a  variable  gain  control  for  the  vertical  or  current 
axis.  By  changing  the  gain  it  is  possible  to  determine  a  noise 
level  for  each  sample  and  circuit  combination.  The  noise  level 
observed  on  the  oscilloscopes  is  a  subject  presently  being 
investigated  at  the  TREES  laboratory;  however,  investigation  and 
analysis  so  far  have  not  revealed  any  way  of  reducing  the  noise 
level  with  the  present  electronic  equipment.,  screen  room,  and 
flash  X-ray  machine. 

The  noise  level,  therefore,  is  an  inherent  limitation  of 
the  flash  X-ray  system  and  its  instrumentation.  Experimental 
results  otherwise  obtainable  are  obscured  by  noise  if  the  magni¬ 
tude  of  the  results  is  at  or  below  the  noise  level. 

e.  The  Radiation  Pulse 

The  shape  of  the  radiation  pulse  can  be  measured  and  recorded 
on  an  oscilloscope  having  a  Polaroid  camera  attachment.  The 
shape  may  vary  from  pulse  to  pulse  and  may  be  adjusted  to  some 
extent  by  changing  various  machine  parameters  prior  to  pulsing 
the  X-ray  machine.  Although  the  radiation  pulse  is  not  a  square 
wave  in  form,  its  concentrated  form  and  short  time  duration  make 
the  approximation  of  a  square  wave  for  computation  purposes  a 
reasonable  one. 

f .  Irradiation  of  Diode  Samples 

Irradiation  of  a  commercial  diode  (1N459)  produced  the  tran¬ 
sient  current  pulse  shown  in  Figure  20-c.  The  maximum  transient 
current  value  is  .004  ma.  The  time  required  for  the  diode  to 
return  to  steady  state  conditions  is  of  the  order  of  5  micro¬ 
seconds.  Maximum  amplitude  was  achieved  at  1  microsecond.  The 
reason  for  the  continued  rise  of  the  transient  current  after  the 
pulse  duration  is  not  clear. 

Since  the  above  diode  has  a  characteristic  curve  which  matches 
that  of  the  fabricated  titanium  dioxide  diode,  it  was  expected 
that  results  obtained  by  irradiating  the  titanium  dioxide  would 
be  comparable  to  those  obtained  above. 

With  the  same  scope  settings,  circuitry,  bias,  etc.  the 
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Figure  20. 
(a -unpotted . 


Commercial  Diode  (1N459)  Transient  Current 
b-in  partial  vacuum,  c -potted  in  resin  epoxy) 
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titanium  dioxide  diode  was  irradiated.  No  noticeable  transient 
current  was  recorded.  The  gain  was  increased  to  the  maximum 
allowable  by  the  existing  noise  level.  Irradiation  at  this  set¬ 
ting  again  produced  no  detectable  transient  current.  At  no  time 
did  any  of  the  noise  oscillations  exceed  a  value  of  .001  ma. 
Repeated  irradiations  had  no  apparent  effect  on  the  diodes,  and 
results  were  consistently  repeated. 

To  further  substantiate  the  test  circuitry  and  test  com¬ 
parison,  it  was  decided  to  irradiate  a  variety  of  commonly  used 
commercial  diodes.  The  results  of  these  tests  are  shown  in 
Figures  21  through  23- 
5.  Conclusions  and  Recommendations 

Irradiation  of  titanium  dioxide  diodes  and  other  commercial 
diodes  at  the  flash  X-ray  facility.  Air  Force  Weapons  Laboratory, 
revealed  the  following: 

1)  No  permanent  or  semi-permanent  effects  were  observed  in 
any  of  the  diodes  tested.  The  components  responded  identically 
to  successive  bursts,  and  no  change  in  the  electrical  character¬ 
istics  was  found  subsequent  to  the  irradiation. 

2)  Titanium  dioxide  diodes,  when  exposed  to  pulsed  gamma 

radiation  of  .2  microsecond  duration  and  an  average  dose  rate 
-8  2 

of  .45  x  10~  erg/cm  -sec.,  produced  a  transient  current  which 
was  less  than  the  noise  level  of  the  test  circuitry  and  not 
detectable  by  the  equipment  used.  At  no  time  did  the  transient 
current  and  noise  current  combination  exceed  .001  ma.  Other 
commercial  diodes  showed  current  excursions  between  .0025  and 
.625  ma  for  the  same  400  KV  radiation  pulse. 

3)  If  the  electrical  properties  of  titanium  dioxide  diodes 
can  be  improved  through  further  development  and  the  use  of  thin¬ 
ner  films,  it  appears  that  these  diodes  will  give  a  highly  radia¬ 
tion  resistant  diode.  It  has  been  proven  that  the  titanium  dioxide 
diodes  are  much  more  resistant  to  radiation  both  of  the  transient 
and  permanent  damage  type  than  similar  germanium  or  silicon 
diodes. 

4)  It  is  recommended  that  the  titanium  dioxide  diodes  be 
further  investigated  under  higher  radiation  intensities  to 
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Figure  22.  Transient  Currei 
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determine  accurately  the  shape  and  extent  of  the  transient 
current  resulting  from  pulsed  radiation.  Further  development 
for  the  improvement  of  the  electrical  characteristics  of  titan¬ 
ium  dioxide  diodes  is  also  recommended. 
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SECTION  IV 


MEASUREMENT  OF  THE  CONTINUOUS  SPECTRUM 
OF  A  FLASH  X-RAY  MACHINE 
By  Harold  Southward 

1.  Introduction 

A  complete  discussion  of  the  technique  used  to  measure  the 
continuous  spectrum  of  a  flash  X-ray  machine  is  given  in  AFWL 
TDR-64-36,  June,  1964,  by  H.  D.  Southward.  The  use  of  a  curve 
of  transmitted  X-ray  intensity  versus  attenuator  thickness  to 
calculate  the  spectrum  of  continuous  X-rays  was  reported  by 
L.  Silberstein.  Silberstein 1 s  .experimental  method  consisted 
of  measuring  the  intensity  of  a  well-collimated  beam  of  X-rays 
transmitted  through  various  thicknesses  of  metal  scattering 
material.  Photographic  film  was  used  as  the  radiation  detector. 
The  transmission  curve  was  fitted  to  an  empirical  formula,  and 
the  spectrum  was  calculated  from  a  process  which  basically 
involved  the  inversion  of  a  Laplace  transformation.  (Reference 
22). 

It  is  important  to  consider  the  following  points  when 
Silberstein' s  method  is  used:  (l)  The  response  of  the  radia¬ 
tion  detector  as  a  function  of  photon  energy  should  be  known; 

(2)  The  response  of  the  radiation  detector  should  be  linear 
with  the  incident  X-ray  intensity  for  a  given  spectral  distri¬ 
bution  (i.e.,  it  should  not  be  dose  rate  dependent);  (5)  The 
absorption  coefficient  y  of  the  scattering  material,  which  is 
a  function  of  photon  energy  E,  should  be  accurately  known 
along  with  its  derivative  dy/dE;  (4)  An  idealized  geometry 
should  be  approximated  so  that  only  single  scattering  occurs 
within  the  scattering  material;  (5)  Background  radiation  which 
arises  principally  from  scattering  of  X-rays  off  the  floor, 
ceiling,  walls  and  surrounding  objects  should  be  avoided;  (6) 

And  finally,  an  empirical  formula  must  be  found  which  will  fit 
the  experimental  data. 
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The  method  described  herein  uses  a  scintillation  detector 
to  measure  the  intensity  of  the  X-ray  beam  transmitted  through 
the  scattering  material,  thereby  eliminating  the  objections 
which  tend  to  arise  when  film  is  used. 

2 .  The  Apparatus  and  Basic  Data 

Figure  24  shows  a  block  diagram  of  the  equipment  used  to 
obtain  the  transmission  data.  The  detector  consisted  of  a 
2"  x  2"  Nal(T£)  Harshaw  crystal  mounted  on  a  2"  type  114  ITT 
photodiode.  The  photodiode  was  chosen  because  of  its  relative 
insensitivity  to  variations  in  high  voltage  and  because  of  its 
wide  range  of  linearity.  The  signal  and  its  time  integral  were 
recorded  from  an  oscilloscope  with  a  Polaroid  camera.  Figure  25 
shows  some  typical  data.  The  Nal(T£)  crystal  did  not  follow  the 
shape  of  the  burst  of  X-rays  because  its  0.25  microsecond  time 
constant  was  nearly  the  same  as  the  0.2  microsecond  duration  of 
the  radiation.  However,  the  Nal(Tt)  crystal  has  the  advantage 
of  a  higher  sensitivity. 

The  X-ray  machine  output  was  monitored  with  detection 
apparatus  similar  to  that  shown  in  Figure  24  except  that  a 
Pilot  "B"  scintillator  was  used.  A  signal  and  its  time  inte¬ 
gral  are  shown  in  Figure  26.  The  transmission  data  were  nor¬ 
malized  to  the  output  of  the  monitor  in  order  to  correct  for 
small  variations  of  the  output  of  the  flash  X-ray  machine. 

The  shape  of  the  monitor  signal  was  also  helpful  in  determining 
whether  the  flash  X-ray  machine  malfunctioned  when  it  was  fired, 
and  formed  a  good  basis  for  rejection  of  data.  Typical  trans¬ 
mission  data  are  shown  in  Figures  27  and  28  along  with  their 
empirical  curves  for  the  X-ray  machine  operating  at  400  and 
500  kilovolts  respectively.  The  data  are  represented  by  points 
and  the  empirical  curve  fitted  to  the  data  is  represented  by 
the  smooth  curve. 

5 .  System  Linearity 

A  simple  function  which  approximates  the  shape  of  the 
Flash  ¥ -ray  machine  output  is  shown  in  Figure  29.  If  it  is 
assumed  that  the  radiation  detecting  system  is  a  linear  system. 
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Figure  27*  Transmitted  Intensity  of  X-Rays  Produced  by  a 

Flash  X-Ray  Machine  at  a  Tube  Potential  of  400  kev 
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Transmitted  Intensity 


Figure  28.  Transmitted  Intensity  of  X-Rays  Produced  by  a 

Flash  X-Ray  Machine  at  a  Tube  Potential  of  500  kev 
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(105) 


the  response  of  this  system  to  an  idealized  pulse  is 

t 

R(t)  =  K  ^  F (r)e  "  T)  dT 

o 

Here  R(t)  is  the  system  response,  K  is  the  proportionality 
factor,  F(r)  is  the  idealized  input,  and  0  is  the  time  constant 
of  the  Nal(T-t)  crystal.  Figure  30  shows  a  comparison  of  the 
experimental  response  of  the  Nal  crystal  detector  system  with 
the  theoretical  response.  One  may  conclude  from  this  close 
agreement  that  the  radiation  detector  system  is  linear. 

Strictly  speaking,  this  conclusion  is  not  completely  true. 
However,  the  same  results  can  be  obtained  only  if  the  two 
measuring  systems  used  contained  the  same  nonlinearities. 

This  is  very  unlikely  due  to  the  vast  difference  in  the  char¬ 
acteristics  of  the  radiation  detectors. 

4.  Calculation  of  the  X-Ray  Spectrum 

Various  techniques  for  calculating  the  X-ray  spectrum 
from  transmission  data  have  been  published  (References  22,  23/ 
26  and  27) •  They  all  follow  the  general  procedure  given  by 
Silberstein.  The  basic  data  form  a  transmission  curve  I  (x), 
which  is  obtained  experimentally  by  measuring  the  integrated 
intensity  of  the  X-ray  beam  after  it  has  passed  through  x  cm. 
of  copper.  This  curve  was  fitted  to  an  empirical  formula  and 
the  spectrum  was  calculated  from  a  process  which  involved  the 
inversion  of  a  Laplace  transformation.  The  following  treatment 
is  due  to  Greening  (Reference  26): 

I  (x)  =  I  •-**-*<  Vx  +  d  -  4  d) 

Here  A,  B,  and  d  are  constants  and  I0  is  the  incident  intensity 
of  an  X-ray  beam  upon  the  absorber.  If  F(x)  is  the  intensity 
distribution  of  the  X-ray  beam  per  unit  wavelength  then  the 
transmitted  intensity  l(x)  through  a  material  having  absorption 
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coefficients  m(x)  is 

CP 

l(x)  =  j  F(X)e~M(^xdx. 
Amin 


Here 


j  F(X)dA 
Xmin 


I 

o 


(107) 


dee) 


The  minimum  wavelength  is  determined  by  the  relation 


^min 


he 

E - 

max 


(109) 


where  h  is  Planck's  constant,  c  is  the  velocity  of  light  and 
Emax  t*ie  max^mum  photon  energy  in  the  X-ray  spectrum.  It 
is  more  convenient  to  use  photon  energies  in  calculations. 

By  defining 


M  -  M 


(’'min) 


(j  -  A 


(110) 


it  can  be  shown  that  the  validity  of  equation  (110)  by  examining 

equation  equation  (106).  If  it  is  assumed  that  the  absorption 

coefficient  h(e)  is  a  monotone  decreasing  function  of  the  photon 

-Ax 

energy  E,  then  for  large  x,  l(x)oce“  .  Hence,  A  can  be  found 

from  a  table  of  m(e)  if  E  ^  is  known. 

max 

Now  one  considers 


H*  =  U'(x)  (111) 

It  is  possible  to  consider  the  intensity  distribution  as  a  dis¬ 
tribution  in  m',  9(4'). 

f(X)dX  -  q>Gi')  dn'  (112) 
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Equation  ( 107)  than  becomes 


00 


I(x)  -  j 


-li'x 


<p(m  '  )dM 


(115) 


which  is  a  Laplace  transformation.  By  using  equation  (lC5) 
for  l(x)  and  inverting  equation  (llj)  it  can  be  shown  that  the 
intensity  distribution  of  the  X-ray  beam  per  unit  photon  energy 
is  ~ 


c(E)l(E)  = 


Be 


Ad  +  BVd 


e 


B 


2V  7 r 


(M  -  a)V2 


-Md  dM  .  (114) 
*  HE 


Here  l(E)  is  taken  as  a  relative  distribution  and  e(E)  is  the 
energy  absorption  efficiency  of  the  detector  used  to  measure 
l(x) .  The  energy  absorption  efficiency  of  Nal  crystal  is  shown 
in  Figure  (31),  and  is  defined  as  the  ratio  of  the  photon  energy 
absorbed  in  the  crystal  to  the  incident  photon  energy.  This  cal¬ 
culation  was  performed  at  Argonne  National  Laboratory  using  a 
Monte  Carlo  technique  for  a  1/4"  diameter  collimated  beam  of 
X-rays  hitting  the  end  face  of  a  2"  x  2"  cylindrical  sodium 
iodide  crystal  (Reference  28). 

A  table  of  the  mass  absorption  coefficient  of  copper  and 
its  derivative  as  a  function  of  photon  energy  is  given  in 
Table  IV.  These  values  were  derived  from  those  given  by  Grodstein 
(Reference  29)  using  a  formula  of  the  form 


Loge  <r(E)  * 


A  + 
o 


8 

£ 

i*l 


A  .E' 
1 


(115) 


It  is  necessary  to  evaluate  the  constants  B  and  d  in 
equation  (ll4).  This  is  equivalent  to  fitting  the  transmission 
data  to  equation  (105)*  This  may  be  done  very  simply  by  the 
method  of  Greening  (Reference  26),  by  defining 


-  [lo*e  ^  ]  “  Ax 


(116) 
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Table  IV  Mass  Absorption  Coefficient,  a,  and  its  Derivative 
^jljr  as  a  Function  of  Photon  Energy,  E,  for.  Copper* 


E( raev) 

a(E)  SSL 
v  '  gm 

ff(E) 

'  '  gm 

(from  Ref.  29) 

2 

da  / cm  \ 
aE  vgm-lnev' 

.01500 

76.37 

76.8 

19030 

.02000 

34.86 

34.6 

5772 

.02500 

18.64 

2379 

.03000 

11.07 

11.1 

1154 

.03500 

7.101 

623.5 

.04000 

4.837 

4.83 

364.4 

.04500 

3.457 

226.5 

.05000 

2.572 

2.56 

147-9 

.05500 

1.978 

100.7 

.06000 

1.565 

1.58 

70.93 

.06500 

1.269 

51.49 

.07000 

1.051 

38.35 

.07500 

0.8856 

29.21 

.08000  r 

0.7585 

0.762 

22.69 

.08500  ) 

0.6588 

17.93 

.09000  i 

0.5792 

■  i  *  1 

14.39 

.09500 

Y  0.5149 

11.71 

.10000 

0.4622 

0.461 

9-642 

. 10500 

0.4185 

8.031 

.11000 

0.3819 

6.756 

.11500 

0.3509 

5.735 

.12000 

0.3245 

4.910 

.  12500 

0.3018 

4.234 

. 13000 

0.2822 

(continued) 

3.678 

2 

"  These  data  have  been  calculated  from  a  close  fit  of  the  data  in 
Reference  29.  Note  that  the  data  are  questionable  in  the  third 
significant  figure,  da  is  less  accurate  than  a(E) . 

cH? 


a.?. 
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Table  IV  (cont'd) 


E  (mev) 

2 

n(E)  Cm  . 

0 

cr(E)  cm * 

2 

drr  /  cm  \ 

'  '  gm 

'  '  gm 

(from  Ref.  29) 

dE  ' gm-mev  • 

.1350 

0.2650 

3.215 

.1400 

0.2500 

2.826 

.1450 

0.2368 

2.498 

.1500 

0.2250 

0.2  22 

2.220 

.1550 

0.2146 

1.981 

.1600 

0.2052 

1.776 

.1650 

0.1968 

1.599 

.1700 

0.1892 

L.445 

.1750 

0.1824 

1.310 

.1800 

0.1761 

1.192 

.1850 

0.1704 

1.088 

.1900 

0.1652 

0.9964 

.1950 

0.1605 

0.9148 

.2000 

0.1561 

0.156 

0.8422 

.2050 

0.1521 

0.7775 

.2100 

0.1483 

0.7192 

.2150 

0.1449 

0,6669 

.2200 

0.1417 

0.6198 

.2250 

0.1387 

0.5773 

.2300 

0.1359 

0.3388 

.2350 

0.1333 

0.5038 

.2400 

0.1308 

0.4720 

.2450 

0.1286 

0.4430 

.2500 

0.1264 

0.4165 

.2550 

0.1244 

0.3923 

.2600 

0.1225 

0.3701 

.2650 

0.1207 

0.3496 

.2700 

0.1190 

0.3309 

.2750 

0.1174 

0.3136 

.2800 

0.1159 

(continued ) 

0.2976 
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Table  IV  (cont'd) 


E  (mev) 

2 

f.(E )  cm 

nr  ( E  )  Cm 

2 

da  /  cm  \ 

gn r 

(from  Ref.  29) 

dE  'gm-itev  ' 

.2850 

0.1144 

0.2829 

.2900 

0.1130 

0.2692 

.2950 

0.1117 

0.2566 

.3000 

0.1105 

0.112 

0.2448 

.3050 

0.1093 

0.2339 

.3100 

0.1081 

O.2238 

.3150 

0.1070 

0.2143 

.3200 

0.1060 

0.2055 

.3250 

0.1050 

0.1973 

.3300 

0.1040 

0.1897 

.3350 

0.1031 

0.1825 

.3400 

0.1022 

0.1758 

.3450 

0.1013 

0.1695 

•3500 

0.1005 

0.1636 

.3550 

0.09968 

0.1580 

.3600 

0.09890 

0.1528 

.3650 

0.09815 

0.1480 

•3700 

0.09742 

0.1433 

•3750 

0.09672 

0. 1390 

.3800 

0.09603 

0.1349 

.3850 

0.09537 

0.1311 

.3900 

0.09472 

0.1274 

.3950 

0.09410 

0.1240 

.4000 

0.09348 

0.094 

0.1207 

.4050 

0.09289 

0.1176 

.4100 

0.09231 

0.1147 

.4150 

0.09174 

0.1120 

.4200 

0.09119 

0.1093 

.4250 

0.09065 

0.1068 

.4300 

0.09012 

0.1045 

.4350 

0.08960 

(continued) 

0.1022 

Table  IV  (cont'd) 


E  (raev) 

__2 

rrCE}  cm 

nm2 

cm 

2 

da  /  cnr  \ 

W)  gS~ 

gm 

(from  Ref.  29) 

dfe  'gm-Mev  ' 

.4400 

0.08910 

0.1001 

.4450 

0.08860 

0.09802 

.4500 

0.08812 

0.09608 

.4550 

0.08764 

0.09422 

.4600 

0.08717 

0.09245 

.4650 

0.08672 

0.09075 

.4700 

0.08627 

0.08913 

.4750 

O.O8583 

0.08759 

.4800 

0.08539 

0.08610 

.4850 

0.08496 

0.08468 

.4900 

0.08454 

0.08532 

.4950 

0.08415 

0.08201 

.5000 

0.08572 

0.0834 

0.08075 

.5050 

0.08532 

0.07954 

.5100 

0.08293 

0.07858 

.5150 

0.08254 

0.07726 

.5200 

0.08216 

0.07618 

.5250 

0.08178 

0.07514 

.5500 

0.08141 

0.07404 

.5350 

0.08104 

0.07317 

.5400 

0.08067 

0.07223 

.5450 

0.08051 

0.07133 

.5500 

0.07996 

0.07045 

.5550 

0.07961 

0.06960 

.5600 

0.07926 

0.06877 

.5650 

0.07892 

0.06797 

.5700 

O.O7858 

0.06720 

•  5750 

0.07825 

0.06644 

.5800 

0.07792 

(continued) 

0.06571 
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Table  IV  (cont'd) 


E  (m  ev) 

„_2 

a(E)  HL 
'  '  gm 

2 

a(E)  £2- 
'  i  gm 

(from  Ref.  29) 

2 

dcr  /  cm  \ 

HIT  'gm-mev  ' 

.5850 

0.07759 

0.06500 

.5900 

0.07727 

0.06430 

.5950 

0.07695 

0.06362 

.6000 

0.07665 

0.0760 

0.07695 

.6050 

0. 

0.07663 
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A  plot  of  x/y  versus  y  should  be  a  straight  line  if  no  char¬ 
acteristic  radiation  is  present  in  the  X-ray  beam.  In  practice, 
flash  X-ray  irradiation  studies  are  conducted  in  a  room  shielded 
from  electromagnetic  radiation.  The  radiation  was  filtered  by  a 
50  mil  window  of  copper  which  removed  most  of  the  characteristic 
radiation  from  the  X-ray  beam.  It  can  be  shown  that 


B  = 


and 


d 


c 


2 


(117) 


(U8) 


Here  c  is  defined  as  the  intercept  of  the  straight  line  on  the 
x/y  axis. 

A  typical  x/y  plot  is  shown  in  Figure  32.  The  first  set 
of  values  determined  by  this  method  is  generally  satisfactory. 
The  time-integrated  intensity  distributions  of  a  flash  X-ray 
machine  operating  at  400  and  500  kilovolts  charging  voltage 
are  shown  in  Figures  33  and  34. 

5.  General  Discussion 

The  method  described  for  obtaining  flash  X-ray  spectra  is 
limited  to  situations  in  which  the  shape  of  the  X-ray  pulse  is 
closely  reproducible  from  shot  to  shot.  Experience  has  shown 
that  the  intensity  distribution  changes  with  the  shape  of  the 
X-ray  pulse  in  time  as  shown  in  Figure  26. 

This  method,  employing  essentially  the  same  apparatus,  was 
also  used  for  measurements  on  a  high  intensity  continuous  X-ray 
machine  and  found  to  be  satisfactory. 

In  the  inversion  of  equation  (113)  it  is  assumed  that  u' (x) 
is  a  single  valued  function  of  A,  the  wavelength,  so  that  the 
derived  spectrum  is  unique.  A  study  of  a  table  of  absorption 
coefficients  for  the  various  elements  shows  that  this  method  is 
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Integrated  Intensity  (ergs/cm  -kev) 


Figure  32.  Plot  of  y  Versus  y  for  400-kev  Flash 
X-Rav  Data  Used  to  Determine  the  Co- 


Figure  32.  Plot  of  ~  Versus  y  for  UOO-kev  Flash  X-Ray  Data  Used  to 

Determine  the  Coefficients  B  and  d. 


Figure  3?*  Time -Integrated  Intensity  Spectrum  of  400-kev 
X-Ray  Machine  at  20  in.  From  the  Target.  The 
Total  Energy  was  90  ergs/ cm*  with  the  Geometry 
Shown  in  Figure  24.  The  X-Ray  Beam  was 
Filtered  with  .030  in.  of  Copper. 
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Integrated  Intensity  (ergs/cm^ -kev) 


Photon  Energy  (kev) 


Figure  34.  Time -In teg rated  Intensity  Spectrum  of 
500-kv  X-Rays  Produced  by  a  Flash 
X-Ray  Machine  at  20  Inches  from  the 
Target. 2  The  Total  Energy  was  192 
ergs/cm  .  The  X-Ray  Beam  was 
Filtered  with  .050  Inch  of  Copper. 
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4 

limited  to  spectral  measurements  at  energies  below  the  point 

I  where  pair  production  begins  to  dominate  and  cause  i-i(e)  to 

.V: 

increase  with  further  increase  of  E.  According  to  equation  (ll4) 
the  accuracy  of  computation  is  dependent  on  d(i/dE.  This  deriva- 

,  tive  decreases  in  magnitude  near  the  minimum  in  m(e)  .  The  maxi- 

I 

|  mum  photon  energies  contained  in  a  distribution  therefore  determine 

the  type  of  material  used  as  an  attenuator  in  the  transmission 
method.  The  maximum  photon  energy  to  which  the  method  may  be 
applied  depends  upon  the  attenuator  material  and  desired  accu¬ 
racy  of  the  results. 

A  simplification  of  the  technique  to  include  radiation 
detectors  sandwiched  between  successive  absorbers  is  a  possi¬ 
bility#  but  this  does  not  conform  to  the  idealized  geometry 
described  above.  The  problem  of  inversion  of  equation  (115)  is 
formidable  if  one  tries  to  include  a  radiation  build-up  factor 
to  correct  for  the  deviation  from  the  exponential  attenuation 
law. 

There  may  be  occasion  to  calculate  an  intensity  distribu¬ 
tion  as  a  function  of  the  time  during  the  X-ray  flash.  This 
could  be  accomplished  by  using  a  fast  scintillator#  such  as 
plastic  scintillator,  as  a  primary  radiation  detector.  Again# 
it  will  be  necessary  to  know  the  energy  absorption  efficiency 
of  this  detector,  c(E),  which  may  not  be  readily  available.  The 
measurement  of  a  time -dependent  spectrum  will  require  a  higher 
degree  of  reproducibility  of  the  flash  X-ray  pulse  than  the 
measurement  of  a  time -in teg rated  spectrum. 

Finally#  if  radiation  intensity  is  required,  it  is  possible 
to  calibrate  the  Nal(T-C)  detector  against  a  radioactive  standard 
or  X-ray  machine. 
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SECTION  V 


GAMMA  RADIATION  EFFECTS  ON  DIELECTRIC  MATERIALS 
By  LeRoy  Meyer  and  W.  W.  Grannemann 

1.  Introduction 

Many  studies  have  been  made  of  the  conductivity  induced 
in  dielectric  materials  by  exposure  to  ionizing  radiation.  The 
problem  is  of  interest  because  of  the  information  obtained  on 
the  conductivity  processes  in  dielectrics  and  on  the  application 
of  dielectrics  in  a  radiation  environment.  The  first  investi¬ 
gations  were  made  by  F.  T.  Farmer  in  1942#  and  were  followed 
with  experiments  by  John  F.  Fowler,  A.  Rose,  and  many  others. 

In  these  works  the  effects  were  treated  in  terms  of  conventional 
band  theory  (which  uses  the  model  of  trapping  and  recombination 
levels  in  the  forbidden  gap  to  account  for  the  observed  effects) . 
More  recently  the  "hopping"  model  has  been  suggested  as  being 
more  appropriate,  the  magnitude  of  the  mobility  being  the  deter¬ 
mining  factor  as  to  which  model  is  used  (Reference  JO). 

The  dielectric  materials  may  be  divided  into  two  classes, 
organic  and  inorganic.  The  organics  may  be  further  broken 
down  into  crystalline  and  polymeric  materials.  The  present  work 
is  concerned  primarily  with  a  study  of  the  basic  mechanisms  of 
radiation  induced  conductivity  (which  will  be  referred  to  as  RIC 
hereafter)  in  the  polymer  materials.  The  polymers  were  chosen 
for  most  of  the  experiments  because  they  were  more  readily 
available,  and  their  mechanisms  of  conduction  are  not  yet  com¬ 
pletely  understood. 

Both  transient  and  steady  state  experiments  were  performed. 
The  basic  mechanisms  involved  are  the  same  when  proper  time  con¬ 
stants  and  lifetime  are  taken  into  account.  In  many  cases  the 
steady  state  experiments  are  easier  to  perform  and  have  fewer 
side  effects  to  take  into  account. 

The  materials  used  were  polyethylene,  teflon,  and  mylar. 

Of  these  three  mylar  shows  greater  resistance  to  permanent  deunage 
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and  has  the  advantage  over  the  others  of  being  available  in 
film  form  (Reference  3l)»  Also,  it  is  a  pure  chemical  compound 
of  constant  chemical  composition.  The  other  materials  gave 
more  variability  in  RIC.  Teflon  is  probably  the  least  radiation 
resistant,  with  properties  similar  to  polyethylene.  Polyethy¬ 
lene  is  the  most  extensively  used  insulator  of  the  three;  most 
cables  have  polyethylene  insulation. 

As  will  be  shown,  the  RIC  in  polymeric  materials  is  low, 
and  special  precautions  are  necessary  in  making  measurements  to 
eliminate  the  many  side  effects  such  as  noise,  cable  effects, 
secondary  emission,  air  effects,  and  circuit  time  constants. 

The  radiation  sources  used  in  the  experiments  were  the  Ce- 
137  gamma  calibration  source,  the  250  KV,  15  MA  DC  X-ray  machine, 
and  the  600  KV,  2000  amp.,  0.2  Msec  flash  X-ray  machine.  A 
range  of  dose  rates  was  available  from  these  three  sources  from 
about  50  r/hour  to  10^  r/sec. 

The  radiation  effects  observed  in  dielectrics  are  closely 
associated  with  the  dielectric  properties.  These  properties 
(thermal,  structural,  electrical,  etc.)  will  be  discussed  at  the 
beginning  or  this  section  wher<2  the  ones  most  sensitive  to  radia¬ 
tion  will  be  pointed  out. 

The  Kail  effect  and  magnetoresistance  were  investigated  as 
a  means  of  measuring  the  mobility  of  charge  carriers  in  a 
dielectric.  These  mobilities  had  previously  only  been  predicted. 
The  Hall  effect  method  was  successfully  used  for  measuring  the 
mobilities  in  polyethylene  when  the  material  was  under  the  irra¬ 
diation  of  gamma  rays.  In  connection  with  these  experiments 
some  of  the  theory  of  galvanomagnetic  effects  in  intrinsic  poly¬ 
mers  is  presented. 

Experiments  were  also  performed  on  various  polymer  materials 
to  determine  the  effects  of  surface  area,  surface  roughness, 
electrode  material,  and  the  direction  of  radiation  relative  to 
the  sample  on  the  RIC  in  the  material.  Photovoltaic  effects 
were  also  observed.  The  long  tails  of  RIC  decay  in  polyethylene 


94 


were  measured  and  related  to  the  relaxation  time  of  the  material. 
The  RIC  in  the  air  surrounding  a  dielectric  sample  was  also  inves¬ 
tigated  and  found  to  vary  with  pressure  and  temperature. 

2 .  The  Properties  of  Dielectrics 

a.  Electrical  Properties 

The  important  electrical  properties  are:  the  dielec¬ 
tric  constant,  volume  resistivity,  surface  resistivity,  the  dielec¬ 
tric  strength,  the  dissipation  factor,  and  the  loss  factor.  Many 
of  these  electrical  properties  are  difficult  to  measure  under 
radiation  environment  for  the  resistivity  may  not  be  linear  under 
all  radiation  conditions. 

The  dielectric  constant  (c)  can  be  shown  to  be  related  to 
the  polarizabilities  of  the  material  arising  from  electronic, 
ionic  and  dipolar  contributions . 

The  dielectric  strength  depends  on  the  breakdown  potential 
of  the  material,  which  varies  with  material  thickness,  tempera¬ 
ture,  surrounding  media,  gas  and  voids  in  the  dielectric,  chemi¬ 
cal  deterioration  and  radiation  exposure. 

Resistivity  or  conductivity  is  discussed  in  detail  in  a  later 
section.  The  loss  and  dissipation  factors  were  not  observed  in 
this  study. 

The  properties  of  some  of  the  dielectric  materials  used  in 
:T  study  are  shown  in  Table  V.  Electrical  breakdown  is  initiated 
by  some  action  of  electrons  or  ions  at  voids  in  the  electrode¬ 
polymer  interface  and  grows  by  an  intermittent  discharge  mechanism. 
The  breakdown  path  slowly  meanders  through  the  sample  over  a  period 
of  hours  with  little  reference  to  the  overall  voltage  gradient. 

It  then  usually  terminates  in  the  formation  of  a  highly  conductive 
path  through  the  material  (Reference  32).  Under  radiation  environ¬ 
ments  the  electrical  breakdown  may  originate  from  the  ionizing 
effect  of  radiation  on  the  surface,  chemical  changes  within  the 
material,  or  internal  heating  due  to  gamma  radiation.  These 
reactions  then  lead  to  intrinsic  electrical  breakdown  by  electron 
avalanche,  thermal  breakdown  from  gamma  heating,  or  electric  dis¬ 
charge  in  cavities  due  to  gas  evolution  from  chemical  impurities. 
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of  Some  Dielectric  Materials 


b.  Structural  Characteristics 

The  structural  parameters  are  difficult  to  interpret 
in  terms  of  dielectric  qualities.  It  appears  that  structures 
that  are  stable  to  heat  are  also  stable  to  radiation. 

The  important  properties  are  tension,  compression,  bending, 
shear,  torsion,  elastic  strength,  fracture  strength,  ductility, 
resilience,  elongation,  tensile  strength,  and  Young's  modulus. 
These  will  not  be  discussed  further  except  for  a  few  comments 
from  the  literature.  For  example,  the  failure  of  dielectrics 
is  usually  due  to  physical  deterioration  and  to  mechanical 
rather  than  electrical  properties.  They  often  provide  good 
electrical  insulation  even  though  their  mechanical  properties 
are  severely  damaged  by  radiation. 

c .  Chemical  Changes  in  Polymers  Due  to  Irradiation 
Chemical  deterioration  may  lead  to  thermal  decomposi¬ 
tion  and  dielectric  failure.  An  arc  produces  ozone,  and  its 
subsequent  reaction  with  the  dielectric  causes  chemical  deteri¬ 
orations.  These  are  usually  long  term  effects. 

Typical  chemical  reactions  are  classified  as  molecular  chain 
scission,  cross-linking,  and  excitation  (Reference  35).  These 
processes  are  given  below  for  polyethylene  (which  received  a  dose 
of  about  109  ergs/gm  c).  The  molecular  chain  scission  is  given  by 

H  H  H  H  H  H  H 

-C-C-C  +  y  n  -  C-C  +  -C  *  -C  (119) 

III  II 

H  H  H  H  H 

The  molecular  cross-linking  reaction  is  given  by 

H  H  H  H 

I  I  I  I 

-C-C-C-C-  +  y  n  - 

K  H  H  H 
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H  H  H  H 

-C-C-C-C  +  (120) 

i  «  I  (  4 

H  H  H 
H  H  H 

•  <  i  i 

-C-C-C-C 

rill 

H  H  H  H 


r 


Each  bond  indicated  by  a  dash  is  formed  by  a  pair  of  electrons. 

The  energy  required  to  remove  these  electrons  is  of  the  order 
of  3  to  5  ev.  However,  one  can  remove  adjacent  hydrogen  atoms 
and  form  a  series  of  double  bonds  such  as  equation  (119).  As 
the  number  of  double  bonds  increase  the  binding  energy  of 
these  electrons  decreases  due  to  resonance  phenomena. 

One  consequence  of  double  bonds  is  that  the  material  becomes 
colored.  For  a  sequence  of  3  double  bonds  absorption  is  in 
the  ultraviolet?  with  4,  the  material  is  yellow?  for  8  and  9  it 
is  red?  and  for  20  it  is  brown  or  black. 

d.  Thermal  Properties 

Although  the  resistivities  are  high,  they  do  change  with 
temperature  in  the  same  manner  as  in  semiconductors.  The  equa¬ 
tion  for  resistivity  is 

P  =  PQ  exp  (  +  |^)  (121) 

where  E  is  the  activation  energy. 

The  source  of  heat  in  dielectrics  is  the  conduction  currents 
from  the  free  electrons  present,  and  energy  absorbed  from  the  lag 
between  charge  and  applied  voltage  in  an  alternating  field.  Heat 
is  a  function  of  electric  field  in  both  cases,  as  well  as  frequency 
in  the  second  case.  In  some  semiconducting  polymers  the  ratio  of 
thermal  conductivity  to  electrical  conductivity  is  much  less  than 
for  ordinary  semiconductors. 

e.  Optical  Phenomena 

The  prominent  optical  phenomena  is  florescence,  which 
is  associated  with  electron-recombination  processes.  The  opti¬ 
cal  efficiences  reported  in  the  literature  are  of  the  order  of 
1#  to  10#  depending  on  the  material  used. 

f .  Observed  Effects  when  Insulators  are  Irradiated 

The  observed  effects  reported  in  the  literature  are:  gas 
evolution  from  organics,  discoloration,  change  in  conductivity, 
embrittlement  or  softening,  and  change  in  mechanical  properties. 

The  effects  are  both  transient  and  permanent. 
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The  effects  vary  widely  depending  on  composition,  type  of 
material  (organic  or  inorganic),  type  of  radiation,  dose  rate, 
total  dose,  environment,  amount  of  material,  and  amount  of  energy 
absorbed.  The  effects  caused  by  interaction  of  these  factors  can 
be  extremely  complex.  Many  of  these  are  discussed  in  detail  with 
experimental  results  in  Section  6. 

Two  important  applications  of  dielectric  materials  are  in 
cables  and  capacitors.  In  cables  the  secondary  emission  signals 
are  dominant  over  RIC  for  high  dose  rates .  RIC  is  dominant  at 
low  dose  rates  (Reference  34) .  For  capacitors  in  a  gamma  radia¬ 
tion  environment,  little  or  no  change  in  the  value  of  capacitance 
has  been  reported.  This  is  what  might  be  expected  if  one  looks 
at  the  equation  for  capacitance 

C  =  (122) 

when 

€  =  dielectric  constant 

A  =  area 

d  *  plate  separation  distance 

From  equation  (122)  it  is  clear  that  any  change  in  C  would  have 
to  be  due  to  a  change  in  e .  The  c  can  be  related  to  the  polari¬ 
zabilities  of  the  material  through  the  susceptance.  The  polari¬ 
zability  can  be  made  up  of  three  parts:  electronic,  ionic,  and 
dipolar  contributions.  Any  change  in  C  would  most  likely  come 
from  a  change  in  the  dipolar  polarization.  A  build-up  of 
polarization  has  been  reported  in  the  literature  and  also  observed 
in  experiments  performed  on  polyethylene.  This  build-up  is 
small.  Even  when  combined  with  the  other  polarizabilities  the 
overall  change  in  c  is  quite  small.  Thus  any  noticeable  change 
in  C  is  not  likely. 

A  more  significant  factor  in  the  application  of  capacitors 
in  a  gamma  radiation  environment  is  the  leakage  conductivity.  Th<- 
leakage  time  constant  (RC)  is  in  the  millisecond  range  for  most 
circuits.  R  is  the  short  circuit  leakage  resistance  of  the 
capacitor. 
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3.  Theory 

a.  Phenomenon  of  Radiation  Induced  Conductivity  in 
Dielectrics 

There  are  two  aspects  to  the  phenomenon  of  RIC  in 
dielectrics.  The  first  consists  of  determining  the  nature  of 
the  charge  carriers,  how  they  form,  how  they  recombine,  and 
what  their  mobilities  are.  The  second  aspect  is  that  once  the 
properties  are  known  one  should  be  able  to  predict  the  behavior 
of  the  currents  as  a  function  of  radiation  intensity,  applied 
voltage,  properties  of  the  material,  and  temperature.  A  sum¬ 
mary  of  a  few  relevant  experimental  facts  that  determine  the 
model  to  be  used  will  be  given  first. 

Experimental  evidence  has  been  obtained  on  molecular 
crystals  which  indicates  that  carriers  are  formed  only  on  the 
surface  of  the  crystal  or  at  electrodes  by  excitons  (Reference 
30) .  A  modified  electron-band  of  solids  has  been  used  to 
explain  the  experimental  transport  properties.  In  this  theory 
the  intermolecular  binding  is  weak,  and  the  overlap  is  small, 
which  leads  to  narrow  energy  bands  and  therefore  implies  small 
carrier  mobilities. 

The  properties  related  to  conduction  of  current  in  polymers 
is  not  nearly  so  well  established.  The  mobilities  are  always 
low,  and  these  materials  are  usually  P-type.  Also,  there  is  no 
evidence  that  a  band  model  is  applicable  to  polymers.  Any 
interpretation  of  experimental  conductivity  and  activation 
energy  of  band  model  energy  gaps  is  not  yet  justified,  although 
it  is  often  done  in  the  literature.  Only  hole  conduction  and 
trapping  effects  have  been  observed  in  the  pure  materials. 

The  atoms  in  a  polymer  structure  are  held  together  by  the 
weak  forces  of  the  covalent  bonds.  What  is  left  is  an  environ¬ 
ment  of  free  electrons,  ionized  molecules,  and  free  radicals. 
These  are  free  to  move  about  within  the  material  and  may  recom¬ 
bine  in  various  ways. 

The  phenomena  that  occur  ir.  organics  that  do  not  occur  in 
inorganics  arer  changes  in  dielectric  properties,  increased 
sensitivity  to  environmental  interactions,  and  gas  evolution. 
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Because  of  covalent  bonding  between  atoms  most  polymers 
exhibit  dielectric  properties.  The  conductivity  of  organic 

4  -21 

polymers  is  now  known  to  range  from  10  to  10  mho/cm,  a 
25 

range  of  10  .  They  can  thus  be  classed  as  semiconductors  as 

well  as  insulators.  Most  of  the  polymers  in  the  higher  part 
of  the  conductivity  range  have  filler  materials  added  during 
the  manufacturing  process.  They  have  future  possibilities  as 
semiconductors  since  one  could  make  film,  and  mold  shapes  as 
components  of  semiconductor  devices.  There  the  most  serious 
limitations  would  be  high  temperature  tolerance. 

There  is  remarkable  agreemsrt  between  different  experimenters 
who  have  worked  with  commercial  specimens  of  polymers  (plastics). 
The  currents  and  exponents  that  have  been  measured  are  found  to 
give  fair  agreement.  There  seems  to  be  some  fundamental  property 
of  the  basic  molecular  structure  that  is  more  important  than  the 
impurity  content  (Reference  35) • 

The  semicrystalline  polymer  may  be  considered  to  be  a  mix¬ 
ture  of  crystallites  and  intervening  amorphous  regions.  It  may 
be  regarded  as  a  system  of  rigid  particles  (crystallites)  dis¬ 
persed  in  a  rubbery  medium  with  some  type  of  periodicity.  They 
are  not  crystalline  in  the  sense  that  their  molecules  are 
periodically  arranged  in  space,  such  as  in  an  NaCl  crystal,  but 
rather  that  the  polymer  chains  must  be  folded  somewhat  as  shown 
in  Figure  35#  (Reference  3 b).  The  calculated  dimenions  of  the 


Figure  35*  Folding  of  Polymer  Chains 
in  Laminar  Crystals 
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crystallites  were  much  less  than  the  length  of  the  polymer 
molecules.  Thus,  the  f ringed-micelle  model  of  the  structure 
of  polymer  crystals  no  longer  holds  (Reference  32).  Also 
mentioned  in  this  reference  is  the  fact  that  electron  micro¬ 
scope  pictures  have  shown  a  screw  axes  of  spiral  growth;  and 
electron  diffraction  patterns  have  shown  that  the  chain  axes 
of  the  molecules  are  normal  to  the  large  dimensions  of  the 
platelets.  Therefore,  while  the  molecules  themselves  are  not 
periodic,  the  intramolecular  arrays  or  chains  in  the  material 
might  be  said  to  have  a  pseudo  periodicity.  A  comparison  can 
be  made  between  the  chainlike  molecules  and  the  elongated 
crystallites  in  a  wire.  They  both  try  to  align  themselves 
parallel  to  the  fiber  axes. 

The  properties  of  the  polymer  are  not  unlike  that  of  the 
complete  relaxation  of  lattice  periodicity:  that  is,  the 
melting  of  a  solid.  When  melting  occurs,  the  mechanism  of 
transport  is  not  clear.  One  possibility  is  that  the  band  model 
is  replaced  by  the  so-called  "hopping"  model.  The  term  "hop- 
pincf'  has  been  given  to  carrier  motion  observed  in  certain  ionic 
cryptals  whose  electron-phonon  coupling  is  determined  by  elec¬ 
trostatic  interaction  between  the  electrons  and  ions.  In  non¬ 
polar  crystals  the  electron-phonon  coupling  is  determined  by 
interactions  between  electrons  and  induced  molecular  dipoles 
and  is  much  smaller.  Neither  the  "hopping"  model  nor  the  band 
model  is  completely  satisfactory.  The  magnitude  of  the  mobility 
has  been  used  to  determine  the  applicability  of  the  appropriate 
model.  When  the  mobility  is  greater  than  1  cm  /v-sec,  the  band 

theory  is  more  appropriate,  and  when  the  mobility  is  less  than 
o 

1  cm  /v-sec  the  "hopping"  model  is  more  suitable.  In  general, 
when  the  band  model  scattering  free  path  is  equal  to  or  smaller  ' 
than  the  lattice  spacing,  a  "hopping"  process  is  favored 
(Reference  30). 

H.  A.  Pohl  has  shown  that  the  conductivity  in  polymers  is 
electronic  and  not  ionic.  He  performed  the  following  experiment: 
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Over  a  seven  day  period  a  sufficient  current  was  passed  through 
the  polymer  to  electrolyze  every  carrier,  hydrogen,  or  oxygen 
atom  in  it.  The  resistivity  of  the  sample  did  not  change  dur¬ 
ing  this  period.  It  was  concluded  that  electrolysis  does  not 
occur.  The  conductivity  is  not  ionic  and  must  he,  most  probably, 
electronic.  That  the  conductivity  is  electronic  in  nature  also 
follows  from  the  measured  ability  to  pass  much  more  charge 
through  a  sample  than  it  would  store  as  a  capacitor.  Charge 
transport  by  ions,  rotation  of  dipoles,  or  other  chemical  pro¬ 
cesses  was  ruled  out  for  the  same  reason  by  G.  T.  Cheney  et  al, 
and  they  also  observed  that  RIC  is  a  linear  function  of  voltage 
and  does  not  depend  on  the  electrode  material  (Reference  37)- 

It  has  thus  been  established  beyond  reasonable  doubt  that 
the  polymers  are  electronic  conductors  and  that  they  have  an 
exponential  temperature  dependence.  Mathematically  the  band 
theory  equations  describe  the  semiconductivities  of  polymers. 

Holes  have  greater  mobilities  than  the  electrons.  This  does  not 
mean  that  these  materials  are  intrinsic  semiconductors  in  the 
sense  of  band  theory,  but,  rather,  that  mathematically  a  good 
fit  can  be  obtained  using  such  equations.  The  application  of 
band  theory  to  materials  of  low  mobility  is  questionable.  In 
the  light  of  these  facts  the  problem  of  conductivity  in  dielec¬ 
trics  shall  be  examined.  A  model  will  be  used  in  which  it  is 
assumed  that  the  charge  carriers  behave  as  particles  possessing 
characteristic  diffusion  constants  and  mobilities. 

Experimentally  a  positive  Hall  coefficient  has  been  obtained 
for  polyethylene.  This  could  arise  in  two  ways:  l)  an  intrinsic 
material  in  which  the  mobility  of  the  positive  charge  carriers  is 
greater  than  the  mobility  of  the  negative  carriers,  and  2)  an 
extrinsic  material  in  which  the  number  of  positive  charge  carriers 
i.s  gi~eater  than  the  number  of  negative  carriers.  This  would 
require  that  the  impurity  concentration  be  greater  than  the  host 
atom  concentration,  and  is  not  practical.  Assuming  then  that  the 
material  is  intrinsic,  one  can  calculate  the  concentration  of 
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holes  and  electrons  from  equation  (123)  below,  where  the  rest 
mass  of  the  electron  is  used  for  m*  . 


where 

n  =  density  of  electrons 
p  =  density  of  holes 
m*  =  effective  mass 
K  =  Boltzmann  constant 
T  =  temperature  (degrees  Kelvin) 
h  -  Planck's  constant 


E  =  activation  energy 


(123) 


The  equation  for  conductivity  in  intrinsic  material  is  given 
by 

cr  =  |  e  |  n  uh  (1  +  b)  (124) 


where 

e  =  electronic  charge 
=  hole  mobility 

b  =  p  /p.  ratio  of  the  electron  mobility  to  hole 
mobility 

Experiments  Li/,  E.  t.  Con-ad  -1.  found  that  ti«-  conductivity 
of  a  dielectric  in  a  radiation  environment  may  be  described  by 
the  equation 


a  =  ao  +  ai  +  (Reference  31 )  (125) 

where 

a.  ~  photoconductivity 

<jQ  *  dark  conductivity 

KfA  *  gamma  induced  photoconductivity 

K  *  constant  or  proportionality  of  the  material 

A  *  an  exponent  which  is  characteristic  of  the 
material  and  varies  between  0.5  and  1. 
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The  simple  power  relationship  between  a  and  ?  seems  reliable 
as  a  means  of  predicting  the  behavior  of  dielectrics  over  a  wide 
range  of  gamma  flux.  The  relationship  is  empirical  and  does  not 
lend  itself  to  investigation  of  microscopic  material  properties, 
but  the  design  engineer  may  use  the  relationship  with  reasonable 
confidence. 

S.  E.  Harrison  et  aL  used  equations  similar  to  equation  (125) 
to  explain  the  three  time  intervals  observed  in  RIC  during  a 
transient  gamma  pulse  (Reference  58) . 

RIC  is  a  function  of  dose  rate  as  can  be  seen  in  Table  VI. 
Most  investigators  used  only  a  limited  range  of  dose  rates; 
therefore  the  following  table  is  given  for  data  collected  from 
several  sources  on  induced  conductivity  in  mylar. 


TABLE  VI 

INDUCED  CONDUCTIVITY  IN  MYLAR 


Source 


Flux 

r/sec 


Induced 

Conductivity 

(ohm-om) 


Fowler 

1.3  x 

10" 1 

6  x 

10'20 

Conrad 

1.1 

X 

102 

1.1 

X 

10-16 

Conrad 

3.5 

X 

10* 

5.7 

X 

1C-I5 

Day 

6.5 

X 

10^ 

1.5 

X 

10"15 

H 

1.5 

X 

105 

3.7 

X 

10-1’ 

II 

3.3 

X 

105 

9.7 

X 

10-1* 

II 

7.5 

X 

105 

1.9 

X 

10" 15 

M 

3.8 

X 

106 

1.1 

X 

10“ 11 

b.  Hall  Effect  in  Polymers 

The  Hall  effect  is  one  of  the  most  valuable  tools  in 
the  study  of  semiconducting  materials.  It  could  be  equally 
valuable  in  the  study  of  organic  materials,  except  for  the 


105 


difficulty  that  arises  in  organic  materials  because  of  .the  low 

mobilities.  Mobilities  have  been  reported  for  various  organics 

-6  2  2  2 

that  range  from  i0’°  cm  /volt-sec  to  10  cm4/volt-sec . 

For  intrinsic  material  the  following  relations  exist: 

n  -  p  *  ni  (126) 

where 

n  =  electron  density 
p  *  hole  density 
n^  *  carrier  density 

The  Hall  coefficient  can  be  shown  to  be 

«h  ■  hts{t?e}  cm3/coul  <12?) 

where 

.  m  mobility  of  electrons 
*  mobility  of  holes 

e  *  1.6  x  10-1^  coul 


The  values  of  the  mobility  ratio,  b,  have  been  reported  in 

the  literature  to  range  from  0.25  to  0.9#  the  mobility  of  holes 

being  slightly  higher  than  that  of  the  electrons  (Reference  50) . 

From  equations  (124)  and  (127)  the  relation  between  hole 

mobility,  conductivity  and  the  Hall  coefficient  cam  be  obtained 

as  „ 

OR. 

-128) 


The  Hall  voltage  for  the  conditions  of  constamt  voltage 
bias,  and  for  intrinsic  material,  cam  be  shown  to  be  proportional 
to  the  difference  in  mobility  as  follows. 

Consider  the  Hall  voltage  equation 


(129) 
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where 


I  *  bias  current 

B  *  magnetic  field 

d  »  thickness  of  material 


For  constant  voltage  bias  this  equation  can  be  rewritten  as 


V 

h 


“h  5 


B 

a 


(130) 


where  R  is  the  resistance  of  the  material  to  bias  current.  Remem¬ 
bering  that  R  can  be  written  in  terms  of  cr  by 


R 


L 

ah 


(131) 


where 

L  ■  length  of  the  sample 
A  =  cross  sectional  area 

and  when  equations  (124)  and  (12?)  are  substituted  for  a  and 
the  following  equation  is  obtained. 

vh  -  (1-b)  tih  gp  (132) 


or 


K-ue>K 


where  K  is  a  constant 


VAB 

LcT 


(133) 


Therefore  the  Hall  voltage  of  an  intrinsic  material  for  the  con¬ 
dition  of  constant  voltage  bias  depends  on  the  difference  of 
mobilities. 

The  Hall  voltage  is  very  small  and  must  be  separated  from 
the  displacement  voltage  and  magnetoresistance  effects.  For 
example,  when  a  sample  of  polyethylene  is  placed  in  a  magnetic 
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field  and  a  bias  current  applied,  the  resulting  small  Hall  voltage 
is  difficult  to  read  accurately  even  with  ...  microvolt  meter. 

However,  when  the  experiment  is  performed  under  irradiation  from 
a  constant  source  of  gamma  rays  the  conductivity  increases.  There¬ 
fore,  the  current  increases,  and  the  Hall  voltage,  being  propor¬ 
tional  to  current,  also  increases  and  can  usually  be  read  in  the 
M-volt  range.  If  the  assumption  of  having  intrinsic  material  is 
correct,  then  according  to  equation  (13?)  there  is  a  change  in 
the  mobilities  of  the  holes  and  electrons  under  radiation  conditions. 
One  might  ask  what  the  true  meaning  is  of  a  Hall  voltage  which 

is  measured  by  a  meter  having  1  megohm  impedance  connected  to  a 

12 

dielectric  sample  having  an  impedance  of  10  ohms.  The  meter 
would  appear  to  be  a  short-circuit  across  the  sample.  In  terms 
of  a  voltage  source,  the  equivalent  circuit  would  be  as  shown  in 
Figure  36.  If  10  w  V  are  read  across  a  10^  ohm  meter,  a  current 


u 

/  O  S\- 


T 

v, 

1 


tO  HV 


Figure  36.  Equivalent  Circuit  for 
the  Hall  Effect  in  Polymers 

-12 

of  10  x  10  ampere  would  have  to  be  flowing.  This  would  mean 

12 

that  a  drop  of  10  volts  would  occur  across  the  10  ohms  of  sample. 
Therefore  the  true  Hall  voltage  developed  is  10  volts.  Usually, 
in  circuit  theory  where  high  impedance  sources  are  used  a  constant 
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current  source  would  be  more  appropriate.  In  this  case  one  might 

-12 

say  that  he  has  a  constant  Hall  current  source  of  10  x  10  amperes. 

The  Hall  voltage  always  comes  out  to  be  positive  as  is  dis¬ 
cussed  in  part  6. a  of  this  section  which,  according  to  equation 
(133)#  means  that  the  hole  mobility  is  greater  than  the  electron 
mobility.  This  agrees  with  what  has  been  reported  in  the  litera¬ 
ture  on  dielectric  crystal  samples. 

c.  Magnetoresistance  Effect  in  Polymers 

When  a  sample  of  polyethylene  was  placed  in  a.  gamma  field 
of  190  r/hour  and  a  magnetic  field  was  suddenly  turned  on  in  a 
direction  perpendicular  to  the  direction  of  current  flow  through 
the  material,  a  change  in  resistance  was  noted.  This  is  called 
magnetoresistance.  It  is  caused  by  a  departure  from  the  free 
electron  model,  because  the  electrical  charge  carriers  do  not  all 
have  the  same  velocity.  The  relative  change  of  resistance,  Ap,  with 
the  sample  in  a  transverse  magnetic  field,  is  usually  referred  to 
the  resistance  of  the  sample  in  a  zero  magnetic  field.  This  ratio 
yields  the  magnetoresistant  coefficient  and  is  given  by  the 
expression 


(134) 

where 

PQ  *  resistivitv  without  the  .uaynetrc  field 
A  «  magnetoresistant  coefficient 
B  «=  magnetic  field 

A p  *  change  in  resistivity  when  the  magnetic  field 
is  turned  on. 


This  equation  holds  for  low  magnetic  fields.  A  similar 
formula,  called  Frank's  formula,  refers  the  change  in  resistivity 
to  the  resistivity  when  the  magnetic  field  is  on.  It  is  given  by 


AB 


7 


1  +  CB 


(135) 


where  C  ■  a  constant 
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Frank's  formula  holds  for  semiconductors  vdieLe  C  is  usually 

2 

taken  as  equal  to  m  »  which  agrees  with  values  calculated  from 
the  Hall  effect.  For  polymers  where  the  mobility  of  holes  is 
only  slightly  greater  than  the  mobility  of  electrons  a  more 
correct  value  of  C  might  be 

C  =  Q2b  Mh2  (136) 


where 

b  «  mobility  ratio 

Q  *  a  dimensionless  quantity  which  depends  on  the 

device  dimensions  and  other  factors  not  completely 
understood. 


The  relation  between  the  Hall  coefficient  and  the  magneto¬ 
resistance  can  be  shown  to  be 


f>o 


*h 

P 


2  2 
Q*  Ij 


(137) 


The  expression  for  the  magnetoresistance  coefficient  is  then 


A 


2  2 
v  ° 


(138) 


which  cam  be  expressed  in  terms  of  the  mobilities  by  using 
equations  (124)  and  (127)  and  the  relation  p  *  1/a  as 


A  «  Q2  Uh2  (1-b)2 


(139) 


Experimental  magnetoresistance  data  would  only  yield  infor- 

2 

mation  in  terms  of  the  mobility  product  and  Q  .  If  the  Hall 
voltage  data  are  combined  with  the  magnetoresistance  data  the 
value  of  Q  can  be  obtained.  It  was  about  2  for  the  samples 
used  in  this  series  of  experiments. 
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4 .  Experimental  Techniques 
a.  Steady  State 

Conductivity  measurements  on  dielectric  materials  present 
two  special  problems.  These  are:  l)  to  measure  low  currents  of 
order  10  amperes;  and  2)  to  prevent  electrical  leakage  in  the 
sample  and  measuring  apparatus.  When  the  sample  was  cleaned  with 
alcohol  and  not  touched  by  hand  before  measurements  were  made,  the 
leakage  currents  were  not  a  problem.  A  millivac  MV-07c  DC  micro¬ 
volt-ammeter  was  used  for  the  direct  measurement  of  currents  as 

-12  -6 
low  as  10  amperes,  and  voltages  as  low  as  10~°  volts. 

Room  temperature  measurements  of  static  conductivity  may  be 
misleading  because  of  the  low  currents  involved  and  the  long  time 
constant  of  change  in  the  conductivity  when  the  voltage  is  changed. 
Fowler  gives  the  time  constant  as 


T 


10-14  eE/KT 


(140) 


where  E  is  the  average  level  for  deep  traps.  If  E  is  about  1  ev, 
r  will  be  10,000  seconds  at  20°C. 

When  sufficient  time  is  allowed  between  readings  the  volt- 
ampere  relation  should  be  ohmic  in  both  static  and  induced  con¬ 
ductivity. 

b.  Transient  Conductivity  Measurements 

The  problem  of  observing  trams ient  effects  in  insulators 
is  also  a  difficult  one  because  the  currents  involved  are  of  the 
order  of  microamperes  or  less.  These  currents  are  easily  masked 
by  the  various  side  effects  of  the  flash  X-ray  system.  These 
side  effects  are:  cable  effects,  air  ionization,  electrodes- 
geometry  amd  material,  circuit  time  constant,  and  RF  interference 
of  the  flash  X-ray  machine. 

(l)  Cable  effects 

The  radiation  pulse  produced  in  cables  may  be  of  the 
same  order  of  magnitude  and  have  decay  characteristics  the  same 
as  those  produced  in  the  dielectrics  under  test.  All  jables 
were  shielded  with  2"  lead  bricks  up  to  the  feed-through  connectors 
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in  the  vacuum  system  (see  Figure  37  and  see  Appendix  XI  for  details 
of  the  vacuum  system) .  A  short  bare  wire  was  used  in  the  vacuum 
system  for  support  and  electrical  connection  to  the  sample.  With 
this  cable  arrangement  cable  effects  were  negligible  compared  to 
the  signal  level  of  the  dielectric. 

(2)  Air  ionization 

Air  ionization  around  the  sample  provides  a  shunt 
current  path  on  the  order  of  one  megohm.  A  vacuum  system  will 
reduce  the  air  effect  to  a  minimum  during  an  experiment.  Potting 
a  dielectric  sample  as  is  done  with  other  components  to  reduce  air 
effects,  is  not  practical.  The  air  effects  are  discussed  in  detail 
in  part  6.f. 

(5)  Electrodes 

The  electrodes  consisted  of  1/2"  square  copper  plates 
between  which  the  sample  was  held  by  pressure  contact  in  a  special 
holder  as  shown  in  Figure  38.  Other  types  of  electrodes  were  also 
used  to  determine  whether  they  had  any  effect  on  the  measurements. 
Some  were  painted  on  the  sample  with  copper  print  paint  and  silver 
paint.  Vapor  co«ced  indium  was  also  used.  They  all  gave  similar 
results. 

(4)  Circuit  time  constant 

The  circuit  used  is  shown  in  Figure  39*  The  measured 
rise  time  of  the  system  was  about  30  nanoseconds  and  the  decay 
time  was  about  120  nanoseconds.  Thus,  the  long  tails  which  were 
observed  can  not  be  attributed  to  the  circuit. 

(5)  RF  interference 

The  double-walled  screen  room  reduces  the  RF  interference 
by  about  120  db.  This  effectively  reduces  the  RF  noise  from  the 
flash  X-ray  system,  but  does  not  eliminate  it.  Another  screen 
room  around  the  flash  X-ray  system  would  be  required  to  eliminate 
this  source  of  error. 
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Figure  37.  Arrangement  of  Vacuum  System,  Lead 
Shielding  and  Instrumentation  in  the  Flash 
X-Ray  Experiments 


Figure  38.  Sample  Holder 


Amplifier 


Figure  39*  Circuit  for  Transient 
Conductivity  Experiments 


Millivac 


Figure  40.  Hall'  Voltage  Experimental  Set-Up 
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5*  Dosimetry 

For  practical  application  of  the  test  data,  the  energy 
absorbed  by  the  material  or  the  energy  to  which  it  was  exposed 
must  be  determined.  The  mechanisms  of  energy  deposition  are 
the  photoelectric  effect  and  Compton  scattering.  Most  of  the 
energy  deposited  is  accounted  for  by  the  Compton  scattering.  To 
calculate  the  absorbed  dose,  knowledge  of  elemental  structure  of 
the  sample,  environment,  sample  dimensions,  and  absorption  coef¬ 
ficients  for  each  element  present  are  required.  This  makes  cal¬ 
culation  of  the  absorbed  dose  in  insulating  materials  extremely 
difficult  because  of  the  variety  of  complex  molecular  structures 
present.  Thus,  the  use  of  an  exposed  dose  usually  is  more  prac¬ 
tical  than  the  mathematical  and  quantitative  chemical  analysis 
involved  in  calculating  the  amount  of  energy  absorbed. 

6 .  Experimental  Results 

a.  Hall  Effect  Experiments 

The  experimental  apparatus  used  for  measuring  the  Hall 
voltage  in  polyethylene  is  shown  in  Figure  40.  A  Millivac  micro¬ 
voltmeter  was  used  to  read  the  voltage  to  an  accuracy  of  +  2 
microvolts  on  the  10  microvolt  scale. 

The  1/8"  x  1/4"  x  1/2"  polyethylene  sample  was  mounted  between 
two  copper  plates  in  the  special  holder  as  shown  by  Figure  58.  The 
electromagnet  had  tapered  poles  with  a  hollow  bolt  holding  the 
poles  in  place.  The  1/2"  hole  through  the  poles  served  as  a  means 
of  irradiating  the  sample  with  a  collimated  gamma  beam  from  the 
Ce -137  source.  The  gamma  dose  rate  at  the  location  of  the  sample 
was  190  r/hour.  The  electromagnet  was  capable  of  being  varied 
from  zero  to  two  kilogauss.  The  bias  voltage  source  was  variable 
from  0  to  3000  volts. 

The  bias  current  measured,  with  3000  volts  applied  to  the 
sample,  was  2.8  upamps.  This  was  with  no  radiation  or  magnetic 
field  applied.  When  the  magnetic  field  was  applied  the  current 
dropped  to  2.4  ppamps  due  to  a  magnetoresistance  effect. 
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When  the  radiation  was  turned  on  the  bias  currents  increased 
to  the  10“9  ampere  range.  The  current  would  decrease  by  20 %  to 
40$  when  the  magnet  was  turned  on.  The  magnetoresistance  effect 
on  polymers  will  be  discussed  in  part  6.b. 

A  large  displacement  voltage  was  measured  at  the  Hall  voltage 
probes.  Cable  and  air  effects  also  had  to  be  accounted  for.  The 
voltage  read  at  the  Hall  voltage  terminal  was  positive  (or  negative, 
depending  on  the  polarity  of  the  biasing  current)  for  both  direc¬ 
tions  of  the  magnetic  field.  However  it  was  larger  in  one  direc¬ 
tion  than  the  other.  The  difference  between  these  two  readings 
always  gave  a  positive  Hall  voltage.  The  manner  in  which  the 
difference  is  obtained  is  important  here  because  of  the  many  side 
effects  involved.  If  a  series  of  Hall  voltage  readings  is  taken 
for  different  bias  voltages  with  the  magnetic  field  in  the  N-S 
direction  and  then  with  the  magnetic  field  reversed,  S-N  direction 
the  resulting  averages  will  not  be  consistent  because  of  meter 
drift  and  hysteresis  effects  in  the  sample.  The  magnetic  field 
must  be  reversed  for  each  bias  current  setting. 

With  a  3000  volt/inch  bias  on  a  polyethylene  sample  and  no 
radiation  applied,  a  Hall  voltage  of  4  MV  was  observed.  This  dis¬ 
placement  voltage  was  10  (-1 V  with  no  magnetic  field  and  no  radia¬ 
tion.  When  the  radiation  was  turned  on  the  displacement  voltage 
increased  to  235  MV  with  no  magnetic  field.  When  the  magnetic 
field  was  turned  on  the  displacement  voltage  dropped  to  127  MV, 
due  to  the  smaller  current  caused  by  the  magnetoresistance  effect. 
The  Hall  voltage  was  +  21  mV,  which  was  an  increase  of  a  factor 
of  5  over  the  Hall  voltage  without  radiation.  This  was  at  a  radia¬ 
tion  dose  rate  of  190  r/hour.  At  lower  bias  voltages  the  increase 
was  less.  This  indicates  that  the  mobilities  are  a  function  of 
the  voltage  applied  as  well  as  dose  rate. 

The  Hall  voltage  data  and  conductivity  data  were  used  to  cal¬ 
culate  the  mobilities  of  electrons  and  holes  in  polyethylene.  By 
means  of  equation  (123)  one  can  calculate  the  number  of  carriers 
as  a  function  of  excitation  energy.  Table  VII  gives  these  values. 
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Table  VII.  The  Number  of  Carriers  n  for 
Various  Excitation  Energies  as  Calculated 
Using  Equation  (5) 


E  (electron  volts)  n  (carrier  density) 


0.2 


0.4 

0.5 

0.6 


0.7 

0.8 

0.85 

0.9 

0.912 

1.0 

1.1 


1.5 


8.76  x 

3-3  X 

6.2  x 
12.1  x 

2.35  x 
4.39  X 
6.27  X 
8.46  x 
5.12  x 
1.53  X 

2.08  x 

0.39  X 


7T 

18 

16 

14 

13 

11 


10 

10 

10 

10 

10 

10 


10 

10 

10 


9 

9 

8 


The  Hall  coefficient  was  obtained  from  the  measured  value 
of  Hall  voltage  and  equation  (l29)»  The  mobility  ratio,  b,  can 
be  obtained  by  making  use  of  equation  (127)  and  a  value  of  n  from 
Table  VII.  It  was  found  that  only  values  of  n  for  activation 
energies  greater  than  0.912  ev  would  fit  the  data  for  Hall  vol¬ 
tage  measurements.  For  energies  greater  than  1.1  ev,  b  became 
equal  to  1,  which  indicates  that  electron  and  hole  mobilities 
would  be  equal.  Thus,  the  activation  energy  for  the  polyethylene 
used  in  the  experiment  was  in  the  neighborhood  of  1  ev.  This 
value  falls  within  the  range  of  E  given  by  Fowler  of  0.4  ev  to 
1.5  ev  for  polyethylene. 

By  using  the  trap  depths  of  1  ev,  the  value  for  b,  and  the 
values  for  the  mobilities  of  holes  and  electrons  were  calculated 
making  use  of  equation  (128).  The  results  of  measured  and  cal¬ 
culated  values  are  combined  in  Table  VIII  for  three  different 
bias  voltages. 
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b.  Magnetoresistance  Experiments 

(l)  Polyethylene 

(a)  Mobility  calculations 

From  the  magnetoresistance  data  on  a  .25"  x 

.35"  *  .175"  sample  of  polyethylene  and  equation  (139)  the  values 
2 

of  Q  MeM^  were  calculated  for  several  intensities  on  the  sample. 
They  are  shown  in  Table  IX. 

Table  VIII.  Hall  Voltage  Data  and  Calculated  Values 
of  Mobilities  for  Polyethylene  at  190  r/hour 


Bias 

Voltage 

Volts/in 

Bias 

Current 

-12 

x  10  amp 

Hall 

Voltage 

10~6  volts 

Resistivity  (p) 
13 

x  10  ohm  cm 

Hall  (R,) 
Coefficient 

x  10^ 

1500 

120 

+  10 

1.56 

1.22 

2000 

145 

+  15 

1.87 

1.52 

3000 

175 

+21 

2.02 

1 .62 

Bias 

Voltage 

Volts/in 


Mobility 

Product  (m  M,  ) 
h  e  n 

-^95 - J  x  10'6 

(v-sec) 


Mobility 

Ratio 

b 


Electron  (u  )  Hole  (|J.  ) 
Mobility  Mobility 

2 


cm 


v-sec 


x  1C 


-3 


cm 


v-sec 


x  10 


-3 


1500 

1.02 

0.923 

0.97 

1.05 

2000 

1.18 

0.932 

1.05 

1.13 

3000 

1.71 

0.94 

1.27 

1.35 

The  mobility  product,  MeM^,  calculated  from  the  Hall  voltage 
data  was  greater  than  the  Q2MeMh  product  by  a  factor  of  about  2 
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2 

(see  Table  Vlll).  This  would  indicate  the  Qc  constant  was  about 
2  for  polyethylene. 


i 

* 


Figure  4l.  Magnetoresistance  Experimental  Set-Up 

The  curve  of  current  versus  voltage  with  the  sample  out  of 
the  radiation  field  is  shown  in  Figure  42.  The  magnetoresis¬ 
tance  effect  is  noted  by  the  smaller  current  for  the  curve  plotted 
from  data  with  the  magnetic  field  turned  on. 
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Current  (mia) 


Figure  42.  Curve  of  Current  Versus  Voltage  for 
Polyethylene  with  no  Radiation  Applied 


» 
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Table  IX.  Values  of  Resistivity  and  Mobility 
Product  Obtained  from  Magnetoresistance  Data 
on  Polyethylene  in  a  Radiation  Environment 
of  190  r/hour  Gamma  Flux 


Volts/in 

p  xlO12 

0 

12 

V10 

AP 

Q2u  u.xl0‘6 
e  h 

8600 

6.58 

20.2 

2.07 

0.518 

7150 

6.41 

20.9 

2.27 

0.567 

5730 

5-69 

18.72 

2.29 

0.572 

4290 

4.7 

15-62 

2.32 

0.583 

2860 

3-32 

11-73 

2.53 

0.635 

1430 

1.85 

9-93 

4 .  y- 

1.095 

(b) .  Radiation  Effects  as  a  Function  of  Sample  Orientation 

A  .35"  x  •25"  x  •175"  sample  of  polyethylene  was 
arbitrarily  given  x  y  z  directions  as  shown  in  Figure  43*  The 
bias  current  was  flowing  in  the  x-direction#  and  the  sample  was 
irradiated  first  from  the  y-direction  and  then  from  the  z-direc- 
tion.  The  results  are  shown  in  Figure  44.  The  two  currents  are 
almost  identical  with  the  indirect  current  caused  by  irradiating 
from  the  y-direction  being  slightly  larger.  This  is  explained 
on  the  basis  of  the  energy  absorbed  in  the  y-direction  being 
slightly  greater  than  in  the  z-direction. 

The  sample  was  uniformly  irradiated  at  atmospheric  pressure. 
The  dose  rate  was  190  r/hour. 

(c)  Radiation  Effects  as  a  Function  of  Surface  Area 
Three  samples  were  prepared  with  the  same  cross  sec¬ 
tional  area  and  length,  but  different  surface  areas.  They  are 
shown  in  Figure  43  and  the  results  are  plotted  in  Figure  46. 

All  three  curves  are  within  experimental  error. 
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Figure  4jJ.  Sample  Orientation 

(d)  Radiation  Effects  as  a  Function  of  Surface  Roughness 
Two  samples  of  polyethylene  were  prepared  with  the 
same  dimensions,  but  one  had  grooves  cut  in  the  surface  in  such 
a  way  that  its  actual  surface  area  was  doubled.  They  are  desig¬ 
nated  smooth  and  rough  surfaces  and  are  shown  in  Figure  47?  the 
effects  were  almost  identical  as  can  be  seen  in  Figure  48. 
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Figure  44.  Curve  of  Current  Versus  Voltage  for 
Polyethylene  which  Compares  Results  of  Irra¬ 
diation  from  two  Different  Directions 
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.889  cm 


— .635cm  h- 


2 

Surface  Area  =  2.14  cm 
(not  including  ends) 


2 

Surface  Area  *  1-92  cm 
(not  including  ends) 


2 

Surface  Area  *  I.67  cm 
(not  including  ends) 


Figure  43.  Polyethylene  Samples  with  Three 
Different  Surface  Areas,  but  the  Same 
Volume  and  Cross  Sectional  Area 
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Current  (uua) 


Current  versus  Voltage 
2000  Gauss 
190  r/hour 


Figure  46.  Plot  of  Current  Versus  Voltage  for 
the  Three  Samples  of  Polyethylene  having 
Different  Surface  Areas,  but  the  same 
Volume 
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.12  cm 


2 

Surface  Area  =  4.2  cm 
(not  including  ends) 


(a)  Rough  Surface 


27  cm 


(b) 


2 

Surface  Area  *  2.42  cm 
(not  including  ends) 


Smooth  Surface 


Figure  47-  Smooth  and  Rough  Surface  Sample* 

of  Polyethylene 
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Current  versus  Voltage  f 

1/2"  x  1/4"  x  1/8"  rough**  Polyethylene  Sample  j 

2000  gauss  ? 

190  r/hr  | 

CE-137  source  f 

11  June,  1964  j 

O  -  Smooth  sample  -  Magnet  Off  I 

X  -  Smooth  sample  -  Magnet  On 
O  -  Rough  sample  -  Magnet  Off 
+  -  Rough  sample  -  Magnet  On 


_  -<3> 


1 _ i _ I _ I _ L 


500  1000  1500  2000  2500 

Voltage  (volts) 


Figure  48.  Curve  of  Current  Versus  Voltage  which 
Compares  the  Smooth  and  Rough  Surface  Samples 


(2)  A  1/4"  x  1/4"  x  1/2"  sample  of  teflon  was  irradiated 

using  the  Ce-137  source  at  a  dose  rate  of  190  r/hr.  The  results 

of  the  data  are  tabulated  in  Table  X.  The  resistivity  of  the 

sample  with  no  radiation  field  was  not  calculated  because  the 

-12 

currents  were  less  than  10  ampere.  The  published  value  of 

n  c  _  n 

p  is  greater  than  10  (ohm  cm) 

From  the  magnetoresistance  data  the  mobility  product  was 

obtained.  No  Hall  voltage  data  was  taken  for  teflon  so  b,  m.  , 

?  2  n 
and  Q  cannot  be  found;  but  if  Q  is  around  2,  as  was  found  for 

-4 

polyethylene,  and  b  «  1,  the  mobilities  would  be  around  9  x  10 
2 

cm  /v-sec.  A  curve  of  current  versus  voltage  for  teflon  is  shown 
in  Figure  49- 


Table  X.  Values  of  Resistivity  and  Mobility 
Product  for  Teflon  Obtained  from  Magnetoresistance 
Data  of  Sample  Being  Irradiated  at  190  r/hour 

Gamma  Flux 


V/inch 

Resistivity 

(ohm-cm) 

po(i9o5)x1°12 

Resistivity  in 
Magnetic  Field 
(ohm-cm) 

p  t190"  )xl012 
ra  2000G 

2  -7 

Q  MeMhxl0  ' 

6000 

10.53 

29.3 

4.45 

5000 

10.19 

28.4 

4.47 

4000 

8.93 

22.7 

3-86 

3000 

6.9 

15.87 

3.25 

2000 

4.32 

10.72 

3.70 

1000 

2.04 

5.89 

4.73 

(3)  A  3/64  x  7/32  x  3/8"  sample  of  Ti02  was  irradiated 
with  the  Ce-137  source  with  a  dose  rate  of  190  r/hr.  There  was 
no  observed  change  in  conductivity.  Next  a  magnetic  field  was 
applied  perpendicular  to  the  direction  of  current  flow.  There 
was  no  magnetoresistance  effect  observed  in  or  out  of  the  radia¬ 
tion  beam. 
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Teflon  Sample 


1/4"  x  1/4"  x  1/2" 

2000  Gauss 
190  r/nr 
CE1J7  Source 
12  June,  1964 

O 


Figure  49*  Curve  of  Current  Versus  Voltage  for 
Teflon  in  a  190  r/hr  Gamma  Ray  Environment 
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A  plot  of  current  versus  voltage  applied  to  the  sample  is 
shown  in  Figure  50*  The  calculated  conductivity  for  this  sample 
was  1.01  x  10-5  (ohm-cm)"^.  The  value  given  for  conductivity  (c) 
in  the  American  Institute  of  Physics  Handbook,  is  a  -  5  *  10~^ 
(ohm-cm)-  .  The  purity  of  the  material  was  not  known  in  either 
case.  The  conductivities  depend  somewhat  on  the  purity  of 
crystalline  materials. 

c.  Dose-Rate  Dependence  of  Polyethylene 

The  dose-rate  dependence  for  the  polyethylene  was  found 
to  be  in  the  range  given  by  other  experimenters .  The  current 
dependence  on  dose  rate  is  given  by 


i  oC  RA 

The  exponent  A  was  found  to  be  0.58  in  the  50  to  200  r/hour 
range . 

d.  Photovoltaic  Effects 

When  a  sample  of  polyethylene  is  irradiated  at  one  end 
with  a  collimated  beam  of  y-rays,  a  voltage  appears  across  the 
sample.  The  circuit  is  shown  in  Figure  51*  there  is  no  applied 
voltage  in  the  circuit.  The  polarity  of  the  voltage  pulse  is 
always  such  that  the  irradiated  part  of  the  dielectric  becomes 
negative.  This  effect  was  observed  on  three  different  sources: 
the  Ce-137,  the  200  KV  DC  X-ray  machine,  and  the  flash  X-ray 
machine.  This  effect  occurs  because  the  material  is  non-uniformly 
excited  with  y-rays.  It  is  a  well-known  effect  in  high  purity 
crystal  materials  using  light  excitation. 

The  photovoltaic  equation  in  general  is  given  by 


[5 

AnAp 


- 


undn 


hP  +  wnn 


1 


(141) 


where 

K  «  Boltzmann  constant 
e  *  Electronic  charge 
P  *  Hole  density 
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Current  (m^a) 


Voltage  ( Volta) 

Figure  50 •  Plot  of  Current  Versus  Voltage 
for  Ti©2 
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«**4tt*«5*J  1h*,> 


Since  the  distribution  of  carriers  inside  the  sample  is  not  known, 
this  formula  is  not  too  helpful. 

A  1/4"  x  1/4"  x  1/2"  sample  produced  about  20  microvolts  when 
one  half  of  the  sample  was  irradiated  at  the  Ce-137  source  with  a 
dose  rate  of  190  r/hr.  A  similar  experiment  was  performed  on  the 
250  KV,  15  ma  DC  X-ray  machine  at  a  dose  rate  of  5  r/sec .  The 
photovoltage  was  about  1  millivolt  in  this  case.  Corrections  were 
made  for  cable  pick-up  in  each  case. 

The  same  experiment  on  the  flash  X-ray  at  dose  rates  of  5  x  10 
gave  a  peak  pulse  of  about  1  mv.  The  flash  X-ray  beam  was  found 
to  be  nonuniform  over  the  cone  of  radiation.  At  20.35"  from  the 
X-ray  target  the  beam  was  found  to  vary  as  shown  in  Table  XI. 

Table  XI.  Dose  Measurements  of  the  400  KV 
Flash  X-Ray  Pulse  at  20.35  Inches 
from  Target  and  at  Various  Positions 
across  the  Beam 


Position 

across 

beam 

R 

Milli- 

roentgen 

%  of 

Center  Value 

Center 

33 

100 

1"  Off 

Center 

29 

88 

2"  Off 

Center 

28 

83 

3"  Off 

Center 

10 

28.7 

A  polyethylene  disk  3  11/16"  is  diameter  and  1"  thick  was 
irradiated  in  the  flash  X-ray.  There  was  a  7/8"  hole  in  the 
center  as  shown  in  Figure  52.  Due  to  the  nonuniformity  of  the 
X-ray  beam  there  should  be  a  photovoltage  built  up  with  the 
outer  electrode  becoming  positive.  This  was  observed  when  the 
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center  was  grounded  and  the  outer  electrode  connected  to  a  view¬ 
ing  resistor  as  shown  in  Figure  53(a)*  A  2  millivolt  peak  signal 
was  observed  on  the  first  shot,  with  the  peak  signal  getting 
larger  with  each  successive  shot.  On  shot  5  the  peak  pulse  was 
7.5  mv.  When  the  outer  electrode  was  grounded  (Figure  53(b))  a 
small  negative  pulse  was  observed  at  the  center  electrode. 

These  results  indicate  there  is  a  build-up  of  polarization 
in  the  sample.  It  had  not  been  irradiated  before.  This  build¬ 
up  has  been  observed  in  cables  by  other  investigators  who  have 
found  that  it  continues  until  about  a  1000  R  total  dose  has  been 
absorbed,  after  which  no  further  increase  is  observed. 

This  photovoltage  may  explain  the  pulse  observed  by  other 
investigators  in  experiments  where  they  obtained  a  pulse  with  no 
applied  voltage  and  the  device  received  nonuniform  irradiation . 
e •  Leakage  Conductivity  of  a  Parallel  Plate  Polyethylene 
Capacitor 

In  an  experiment  with  a  parallel  plate  capacitor  arrange¬ 
ment  (Figure  54)  the  capacitor  was  charged  to  200  volts  and  then 
the  battery  removed.  When  irradiated  uniformly  with  a  400  KV 
X-ray  shot,  the  observed  pulse  had  a  very  long  decay  time.  This 
suggests  that  charge  carriers  are  being  released  from  traps.  After 
several  milliseconds  the  monitoring  electronic  system  would  not 
follow  the  pulse;  the  decaying  current  pulse  is  shown  in  Figure  55* 
The  capacitor  size  was  .35"  x  .25”  *  .175"* 

f *  Air  Ionization  Effects  in  the  Flash  X-Ray  Experiments 
The  Hall  effect  measurements  were  not  successful  in  the 
flash  X-ray  experiments  because  of  the  many  side  effects  masking 
the  desired  pulse.  Also,  the  instrumentation  problems  were  much 
greater  with  the  transient  measurements  than  in  the  steady  state 
experiments.  For  the  size  of  samples  used  the  induced  currents 
were  calculated  to  be  in  the  nanoampere  range.  The  largest 
unwanted  side  effect  was  air  ionization.  Since  this  effect  is 
observed  in  circuit  experiments  as  well  as  experiments  on  com¬ 
ponents,  it  was  studied  in  some  detail.  A  high  vacuum  does  not 
entirely  eliminate  the  effect;  it  only  reduces  it. 
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(a) 


(b) 


Figure  53.  Photovoltaic  Experiment  Circuit 
with  Polyethylene  Disk 
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Current  (Ma) 


Irradiated 


A  parallel  plate  capacitor  was  used  with  air  as  the  dielectric. 
The  separation  of  the  plates  was  about  what  the  distance  between 
leads  would  be  for  circuit  components.  The  circuit  arrangement  is 
shown  in  Figure  5 6.  Here,  two  1  centimeter  square  copper  plates 
are  separated  a  distance  of  1  centimeter.  They  were  irradiated  by 
the  400  KV  flash  X-ray  0.2  usee  gamma  pulse  at  various  pressures 
inside  a  vacuum  chamber.  The  results  are  shown  in  Figure  57* 

These  pulses  were  traced  from  Polaroid  pictures  of  the  oxcillos- 
cope  traces  and  put  on  one  plot  for  easier  comparison  at  the  dif¬ 
ferent  pressures.  The  initial  time  constant  was  a  function  of 
the  measuring  circuit  and  is  approximately  2  microseconds.  The 
long  tails  are  observed  between  pressures  of  200  mm  Hg  and  about 
1  mm  Hg.  They  are  explained  by  the  "attachment  time"  of  electrons 
to  oxygen  molecules  in  air.  This  phenomenon  of  the  variation  of 
attachment  times  of  electrons  to  oxygen  molecules  in  air  as  a  func¬ 
tion  of  pressure,  and  in  a  gamma  radiation  environment,  is  covered 
in  detail  in  Section  I  of  this  report. 

The  creation  of  an  electron- ion  pair  in  air  requires  about  32 

2 

ev  of  energy.  The  absorption  coefficient  for  air  is  0.0157  cm  /g 
for  400  KVP  gamma  rays.  When  calculations  are  carried  out  for 
currents  resulting  from  the  shunt  leakage  path  due  to  air  ioniza¬ 
tion  across  components,  they  are  in  the  microampere  range  which 
agrees  with  the  experimental  results  obtained.  The  general  pro¬ 
cedure  for  calculating  the  RIC  in  air  is  as  follows:  The  energy 

2 

of  the  gamma  ray  beam  in  ev/cm  is  multiplied  by  the  absorption 
coefficient  of  air.  This  yields  the  energy  absorbed  in  air. 

Then  energy  absorbed  is  divided  by  32  ev  to  get  the  electron -ion 
pairs  created.  Another  approach  was  used  in  Section  1  of  this 
report. 

The  ions  created  by  the  X-ray  pulse  do  not  gain  much  energy 
for  the  ratio  E/p  <  10,  where  E  is  the  electric  field  across  the 
plates  and  p  is  the  pressure  in  mm  Hg.  For  values  of  the  ratio. 

E/p,  greater  than  10  the  ions  gain  energy,  which  explains  the 
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Figure  57-  Observed  Current  Pulses  for  Parallel 
Plates  with  Air  Dielectric  at  Various  Pressures 
(10KO  Viewing  Resistor,  20v  Bias  Voltage) 
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fact  that  the  tails  get  shorter  when  the  potential  gradient  is 
increased.  If  the  voltage  pulses  for  pressures  of  300  mm  Hg  to 
100  nun  Hg  are  observed  (Figure  57)*  one  sees  that  they  first  go 
positive  and  then  negative.  The  initial  positive  pulse  is  due  to 
injection  current  (ite.,  electrons  flowing  from  ground  through 
the  10  K  viewing  resistor  to  replace  those  knocked  off  as  secondary 
electrons  from  the  plate).  A  space  charge  is  set  up  between  the 
two  plates.  When  the  electrons  return  to  the  plate  due  to  the 
potential  gradient  and  the  fact  that  the  electrons  have  not  yet 
combined  with  oxygen  molecules  in  the  air,  the  pulse  goes  negative. 
At  atmospheric  pressure  and  pressures  below  10  microns  the  pulse 
stays  positive.  This  is  because  the  attachment  times  of  the 
electrons  to  oxygen  molecules  are  much  shorter  at  atmospheric 
pressure  and  the  mean  free  path  is  much  longer  at  low  pressures . 

When  the  polyethylene  sample  is  placed  between  the  plates 
similar  pulses  were  observed.  A  1/8"  x  1/4"  x  1/2"  polyethylene 
sample  was  used.  The  observed  induced  current  caused  by  the 
ionizing  X-ray  pulse  can  be  explained  as  being  entirely  due  to 
air  ionization*  secondary  emission,  circuit  matching  impedance 
and  time  constants,  and  space  charge  effects.  This  also  appears 
to  be  true  of  other  circuit  components  of  similar  size,  such  as 
resistors  and,  possibly,  diodes  and  transistors.  The  temperature 
also  has  an  effect  on  the  pulse  shape.  At  higher  than  room  tem¬ 
peratures  the  pulse  sometimes  went  positive  again  after  being 
negative. 

The  negative-going  pulse  was  not  observed  when  a  smaller 
viewing  resistor  (100  0)  was  used.  With  this  100  Q  viewing  resis¬ 
tor  the  circuit  time  constant  was  much  less  than  the  observation 
time  and  the  resulting  pulses  are  shown  in  Figure  58*  These 
curves  were  also  traced  from  Polaroid  pictures  of  the  pulses 
and  put  on  one  plot  for  easier  comparison.  The  long  tails  are 
again  explained  by  the  attachment  time  of  electrons  to  oxygen 
molecules.  The  initial  pulse  follows  the  X-ray  pulse  and  is  a 
combination  of  currents  from  secondary  emission  replacement  cur¬ 
rent,  RIC,  and  space  charge  effects. 
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urrent  (ua) 
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The  peak  current  as  a  function  of  pressure  is  shown  in 

Figure  59*  This  curve  is  for  the  polyethylene  sample.  The 

maximum  air  effect  can  be  seen  to  occur  around  150  mm  Hg.  This 

was  the  same  point  that  was  found  for  resistors  by  H.  T.  Cates. 

An  equivalent  circuit  of  the  parallel  plate  capacitor 

arrangement  is  shown  in  Figure  60. 

For  the  size  of  the  polyethylene  sample  used  here  the  value 

of  C  was  a  few  MM f  and  the  leakage  resistance  through  the  material 

12 

was  of  the  order  of  10  ohms.  The  calculated  induced  current  in 

the  material  was  -  10  ampere,  which  would  not  be  observed  on  the 

viewing  resistor  of  100  ohms.  The  shunt  air  resistance  is  of  the 
order  of  megohms. 

The  same  two  parallel  plates  were  placed  in  the  constant 
gamma  radiation  beam  of  the  Ce-137  source.  The  currents  observed 
due  to  air  ionization  as  a  function  of  voltage  are  shown  in  Figure 
6l.  The  lower  curve  was  in  a  transverse  magnetic  field  of  2000 
gauss.  Thus,  there  appears  to  be  some  magnetoresistance  effect 
in  air  at  these  currents  levels. 

g.  Correlation  of  Flash  X-Ray  Experiments  with  the  Ce-137 
Gamma  Source  Experiments 

Two  types  of  samples  were  prepared  for  irradiation  in 
both  the  flash  X-ray  and  the  Ce-137  gamma  source.  These  were 
thin  film  parallel  plate  capacitors  made  of  .008  inch  mylar  film 
and  0.005  inch  polyethylene  film.  However,  the  mylar  samples  had 
different  electrode  material.  The  samples'  dimensions  were  such 
that  the  predicted  currents  under  irradiation  in  the  flash  X-ray 
were  in  the  microampere  range .  Microamperes  could  be  read  out 
with  the  instrumentation  system  used,  but  any  current  less  than  a 
microampere  would  produce  signals  near  the  noise  level. 

(l)  Polyethylene  Film  Capacitor 

The  capacitor  was  made  with  a  2  inch  wide  strip  of 
.005  inch  polyethylene  16  inches  long.  A  guard  ring  was  made 
around  the  center  electrode  which  was  made  of  aluminum  foil.  The 
ground  electrode  was  made  of  silver  paint.  The  sample  was  bent 
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Figure  59.  Peak  Currents  as  a  Function  of  Pressure 


Figure  60.  Equivalent  Circuit  for  the 
Parallel  Plate  Capacitor  in  a  Radia¬ 
tion  Environment 


Current  (wia) 


i 
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Figure  61.  Air  Ionisation  Effects  -  Steady  State 
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in  the  form  of  an  accordion  bellows  to  reduce  the  exposure  area 
in  the  flash  X-ray  beam.  The  results  of  the  irradiation  in  the 
Ce-137  source  are  shown  in  Figure  62.  The  same  sample  was  irra¬ 
diated  in  the  flash  X-ray.  The  observed  current  pulse  followed 
the  shape  of  the  X-ray  pulse.  The  maximum  value  of  current  was 
60  microamperes  using  a  bias  voltage  of  22  1/2  volts.  The  guard 
ring  voltage  was  16  1/2  volts.  The  resistivity  was  calculated 
for  the  different  dose  rates  and  the  results  are  shown  in  Table 
XII. 

Table  XII.  Resistivities  at  Different 
Dose  Rates  for  Film  Polyethylene 


Dose 

1  Rate 

Resistivity 

(ohm-cm) 

1.59 

X 

10"2 

r/sec 

2.7  x  lO1^ 

2.78 

x 

10"  2 

r/sec 

1.72  x  1015 

6.52 

X 

10‘2 

r/sec 

1.33  x  1015 

1.5 

X 

106 

r/sec 

3.05  x  1010 

(2)  Mylar  Film  Capacitor 

Two  mylar  parallel  plate  capacitors  were  built.  One 
had  indium  evaporated  electrodes  and  the  other  had  silver  painted 
electrodes.  The  film  was  0.008  inch  thick.  The  surface  area 
was  24  square  inches.  The  capacitors  were  constructed  with  a 
guard  ring  and  positive  electrode  sandwiched  between  2  pieces 
oi  dielectric.  The  ground  electrodes  were  on  the  outside  of  the 
sandwich  as  shown  in  Figure  63* 

The  capacitance  of  this  type  of  construction  was  less  than 
that  for  the  polyethylene  capacitor.  It  was  about  1000  uu  farads. 
The  results  of  the  Ce-137  experiments  are  shown  in  Figures  64 
and  63*  The  sample  with  the  indium  evaporated  contacts  had  about 
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Figure  63.  Mylar  Film  Capacitor 

2  and  1/2  times  less  resistance  than  the  one  with  the  painted 
contacts.  Table  XIII  shows  the  values  of  resistivity  for  mylar 
as  measured  at  the  various  dose  rates.  These  values  are  found 
to  be  in  the  range  given  in  Table  VI  for  conductivity,  which  is 
the  reciprocal  of  resistivity. 


Table  XIII.  Resistivity  of  Film 
Mylar  at  Various  Dose  Rates 


Dose  Rate 
r/sec 

Resistivity 
ohm- cm 

6.4  x 

10"5 

4.87  x  1016 

1.4  x 

10~2 

2.36  x  1016 

2.7  x 

10~2 

9.26  x  1015 

4.3  x 

105 

2.28  x  1011 
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The  currants  were  measured  using  a  10  ohm  viewing  resistor 
in  the  circuit  of  the  flash  X-ray  experiments.  The  circuit  is 
the  same  as  that  shown  in  Figure  5 6  except  that  the  battery  has 
a  bypass  capacitor  for  the  AC  pulse.  The  total  charge  generated 
by  the  X-ray  pulse  was  also  measured  by  using  a  1  megohm  viewing 
resistor.  This  made  the  circuit  time  constant  long  compared  to 
the  pulse  time,  and  the  change  in  voltage  across  the  capacitor 
plates  could  be  observed.  The  charge  is  then  obtained  from  the 
equation 

q  -  CV 

This  value  of  charge  was  then  used  to  check  the  current  values 

measured  with  a  10  ohm  viewing  resistor.  They  checked  very  well. 

The  charge  transferred  across  the  capacitor  in  one  400  KV  flash 

c  _  \  a 

X-ray  pulse  at  a  dose  rate  of  4.3  x  10^  r/sec  was  11.9  x  10 
coulomb. 

*nie  Ce  >137  source  experiments  were  performed  at  atmospheric 
pressure  using  the  guard  ring  technique  to  eliminate  surface 
currents.  The  flash  X-ray  experiments  were  performed  both  at 
atmospheric  pressure  and  in  a  vacuum  of  about  10  microns.  The 
currents  were  almost  identical.  The  currents  obtained  in  the 
vacuum  were  just  slightly  less.  The  difference  in  current  readings 
was  within  experimental  error. 

The  long  tail  of  the  RIC  induced  pulse  was  also  observed  for 
the  silver  painted  mylar  sample  when  irradiated  with  the  flash 
X-ray.  This  tail  is  shown  in  Figure  66.  This  indicates  trapping 
is  involved  in  the  process  of  returning  to  equilibrium  after  the 
X-ray  pulse.  It  is  of  interest  that  this  tail  is  of  the  same 
order  of  magnitude  as  that  observed  for  the  leakage  conductivity 
of  polyethylene  in  part  6*e.  Figure  35* 

7*  Summary  and  Conclusions 

Some  of  the  characteristics  of  insulators  have  been  reviewed 
to  give  a  basic  understanding  of  their  relation  to  radiation 
effects.  The  largest  bulk  radiation  effect  observed  was  RIC 
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(radiation  induced  conductivity).  Air  ionization  is  a  significant 
effect  around  high  impedance  components  in  a  circuit  which  is 
exposed  to  the  air. 

The  Hall  effect  method  and  magnetoresistance  method  were 

developed  to  measure  the  mobility  of  polymeric  materials.  It 

2 

was  found  to  be  around  10  cm  /v-sec  for  polyethylene  and  slightly 
less  than  this  for  mylar.  The  positive  Hall  coefficient  observed 
indicates  the  charge  carriers  are  holes.  From  a  literature  survey 
the  conduction  process  was  found  to  be  electronic  and  not  ionic. 

The  number  of  carriers  produced  was  found  to  be  directly  pro¬ 
portional  to  the  amount  of  radiation  absorbed.  The  mechanism  of 
charge  transport  in  polymers  is  not  yet  completely  understood, 
and  a  completely  satisfactory  model  has  not  yet  been  established. 

It  appears  that  some  combination  of  the  band  model  and  ’hopping" 
model  might  apply. 

In  order  to  predict  the  amount  of  RIC,  the  number  of  charge 
carriers  must  be  known.  This  requires  a  knowledge  of  the  material 
composition,  absorption  coefficients,  density,  geometry,  and 
intensity,  as  well  as  the  type  of  radiation.  This  is  a  complicated 
problem  for  any  but  the  most  simple  geometry  and  controlled  environ¬ 
ment.  This  is  why  empirical  equations  have  been  resorted  to.  The 
difficulty  in  using  these  empirical  equations  is  that  all  the  con¬ 
stants  must  be  known  for  the  material  and  geometry. 

The  effect  of  different  contact  materials  on  dielectrics  under 
irradiation  was  found  to  give  some  variation  in  RIC.  Parallel  plate 
capacitors  made  with  evaporated  indium  electrodes  had  2  1/2  times 
greater  RIC  than  ones  made  with  painted  silver  electrodes.  There 
was  no  observed  RIC  in  TiC^  capacitors. 

A  photovoltaic  effect  is  observed  when  materials  are  not 
uniformly  irradiated. 

The  directional  study  indicated  only  a  slight  difference  in 
RIC  in  a  polyethylene  sample  when  irradiated  from  different  direc¬ 
tions.  These  differences  can  be  related  to  the  amount  of  energy 
absorbed  in  the  material. 


The  maximum  air  ionization  effect  was  found  to  be  at  a 
pressure  of  150  nun  Hg,  The  RXC  in  air  increased  with  decreasing 
pressure  to  about  150  mmHg  and  then  decreased  until  the  pressure 
was  around  10  microns.  It  remained  constant  at  pressures  below 
10  microns.  This  was  in  agreement  with  what  has  been  found  for 
resirtors.  Below  10  microns  the  RIC  was  about  1/5  that  at 
atmospheric  pressure. 

A  guard  ring  on  the  film  capacitors  was  found  to  reduce  the 
air  ionization  effect  at  atmospheric  pressure  to  approximately 
that  of  10  microns  Hg. 
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SECTION  VI 

ALTERNATING  CURRENT  BRIDGE  METHODS 
IN  TRANSIENT  RADIATION  EFFECTS  STUDIES 
By  L.  T.  Boatwright 


1.  Introduction 

Much  effort  and  many  ingenious  instrumentation  schemes  have 
been  pursued  in  transient  radiation  effects  studies,  the  principle 
objective  being  the  determination  from  various  component  responses 
exactly  which  of  the  several  basic  electrical  properties  are  being 
affected  (References  39  end  40) .  For  example,  in  a  capacitor,  the 
electrical  quantity  C  is  quite  dependent  on  the  dielectric  constant 
and  conductivity  of  the  insulating  materials.  Responses  of  active 
components,  such  as  semiconductors,  to  transient  radiation  are 
more  complex  because  charge  carrier  generation,  carrier  lifetime 
and  mobility  are  drastically  affected  (Reference  8).  It  is  well 
recognized  that  external  air  ionization  (and  internal,  depending 
on  the  component)  has  an  important  effect  on  the  operational  per- 
formance  of  the  part. 

2 .  Experimental  Work 

As  shown  in  Appendix  III,  a  Schering  bridge  network  for 
measuring  capacitance  is  relatively  insensitive  to  changes  in 
conductivity.  The  relationship  describing  the  magnitude  of 
changes  is 

-  tan  0  (142) 

where  9  is  defined  as  the  loss  angle  (very  small  for  high  quality 
capacitors),  and  C  and  R  are  the  effective  capacitance  and  series 
loss  resistance,  respectively,  of  the  test  item.  Thus,  a  change 
in  conductivity  as  seen  by  the  bridge  is  suppressed  by  the  tar.  9 
factor;  whereas  a  capacitance  change  should  be  evidenced  by  bridge 
unbalance,  the  magnitude  of  which  is  controlled  by  relationships 
between  the  other  bridge  arms. 
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Experiments  were  conducted  at  the  Air  Force  Weapons  Labora¬ 
tory  TREES  facility  on  a  limited  number  of  components  including 
mica  capacitors,  resistors,  general  purpose  diodes  and  miniature 
transistors.  The  facility  utilizes  a  Fexitron  X-ray  system  capa-j 
ble  of  delivering  a  0.2  microsecond  pulse  up  to  600  KVP .  Radia-  : 
tion  exposure  rate  utilized  in  the  experiments  described  was  \ 
approximately  0.65  x  10^  r/sec.  The  component  holding  fixture 
was  a  general  Radio  Type  674  02  adapter  which  was  connected  to 
an  older  model,  GR  Type  RF  bridge  by  3  feet  of  RG-58  cable. 

This  is  a  substitution  type  bridge  utilizing  a  modified  Schering 
circuit  and  is  rated  from  3  to  50  megacycles.  Signal  source  for 
the  bridge  was  a  GR  Type  1330-A  oscillator,  set  for  maximum  output 
at  20  megacycles  («9V)  .  The  bridge  and  oscillator  were  surroundec 
by  lead  shields  ranging  from  1/2  to  2  inches  thick.  Bridge  output 
was  recorded  on  a  Type  533  Tektronix  dual  beam  oscilloscope  con¬ 
nected  to  the  output  of  the  Tektronix  Type  1121  amplifier,  which 
was  in  turn  connected  to  the  bridge  by  15  feet  of  RG-8  coax  cable. 
Sensitivity  of  this  arrangement  is  limited  by  amplifier  noise, 
typically  about  50  microvolts.  Inclusion  of  3  feet  of  RG-58  cablf 
in  the  bridge  arm  causes  some  balancing  problems  because  of  the 
added  capacity. 

Certain  tests  were  made  to  determine  pickup  when  the  flash  , 
X-ray  unit  was  pulsed.  These  included  the  RG-58  cable  alone  in 
the  bridge  circuit,  with  a  shorted  and  open  test  fixture,  with 
2  inch  lead  blocks  between  test  fixture  and  radiation  source  and  i 
gamma  pulsing  with  and  without  the  test  items  in  place.  i 

3.  Teat  Results  j 

Qualitative  results  only  are  reported  since  the  bridge 
experiment  was  preliminary  and  further  consideration  must  be 
given  to  calibration  details,  transient  response  of  the  recording 
system,  and  data  analysis.  ft 

Figures  67  through  69  ate  responses  cf  a  20  unf  disk  caps-  j 
citor,  i>0  and  100  MMf  mica  capacitors  respectively  at  room  I 
temperature.  The  lower  trace  is  the  output  of  the  plastic  scin-  1 
tillator-photo-diode  X-ray  monitor  described  in  Reference  4l,  and! 
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closely  follows  the  actual  shape  of  the  gamma  pulse.  The  middle 
trace  is  bridge  output  at  null  taken  a  few  minutes  before  firing 
the  pulser.  The  top  trace  is  bridge  output  during  irradiation 
and  occurs  at  the  same  time  as  the  lower  pulse  on  the  dual  beam 
oscilloscope.  The  scope  is  triggered  by  the  lower  pulse. 

Displacement  of  the  upper  trace  during  irradiation  is  not 
bridge  unbalance,  but  can  be  attributed  to  a  form  of  replacement 
current  (Reference  8),  although  the  possibility  of  external 
pickup  cannot  be  completely  ignored. 

Figure  70  is  the  bridge  output  when  a  1N604  Silicon  Diode 
(rectifier)  is  pulsed.  Initial  displacement  of  the  bridge  out¬ 
put  indicates  the  creation  of  excess  charge  carriers.  The  bridge 
unbalance  after  X-ray  pulse  termination  indicates  a  change  in 

effective  capacitance  because  of  redistribution  of  charge  across 

-•* 

the  junction.  This  decays  exponentially  in  about  10  ^  second, 
which  roughly  corresponds  to  carrier  recombination  at  a  rate 
determined  by  the  minority  carrier  lifetime  (Reference  39). 

Figures  71  and  72  are  the  collector-base  junction  responses 
from  a  miniature  transistor.  Type  TO-51,  manufactured  by  TRW 
Semiconductors,  Inc.  The  response  of  Figure  71  is  with  the 
collector  grounded,  and  of  Figure  72  with  the  base  grounded. 

As  expected  the  response  is  similar  to  the  1N604  diode  and 
shows  excess  charge  carrier  generation  initially  and  a  change 
of  junction  capacitance  after  termination  of  the  radiation 
pulse.  The  exponential  decay  exists  for  approximately  5  milli¬ 
seconds  . 

Additional  experiments  are  planned  in  order  to  obtain 
qualitative  data  and  to  resolve  what  may  be  small  interference 
during  the  X-ray  pulse.  For  this  purpose  new  equipment  is  on 
order  and  should  enable  acquisition  of  data  up  to  a  50  mega¬ 
cycle  signal  input  to  the  bridge. 
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Figure  6? 

20  Disk  Capacitor  to  X-Ray  Pulse 
X-Ray  Monitor  Output  «  0.6  x  10?  r/sec 
Output  from  Bridge 
Bridge  Null  Before  X-Ray  Pulse 
0.2  tusec/division 
(Center  and  Upper  Trace) 

0.5  millivolt/division 
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Figure  68 

50  MUf  Mica  Capacitor  to  X-Ray  Pulse 
X-Ray  Monitor  Output  w  0.6  x  10?  r/sec 
Output  from  Bridge 
Bridge  Null  Before  X-Ray  Pulse 
0.2  Msec/division 
(Center  and  Upper  Trace) 

0.5  millivolt/division 
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100  pr.f  Mica  Capacitor  to  X-Ray  Pulse 
X-Ray  Monitor  Output  «  0.6  x  10'  r/sec 
Output  from  Bridge 
Bridge  Null  Before  X-Ray  Pulse 
0*2  usec/division 
(Center  and  Upper  Trace} 

CT5  millivolt/diviaion 


Figure  70 
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TZ  ZN604  Silicon  Diode  to  X-Ray  Pulse 
X-Ray  Monitor  Output  **  0.6  x  10?  r/sec 
Output  from  Bridge 
0.2  Msec/di virion 

0.5  millivolt/division 


161 


Figure  71 

Base -Collector  Junction  Response 
of  a  Type  TO-51  TRW  Transistor 

Lower  Trace:  X-Ray  Monitor  Output  pa  .65  x  10‘r/sec 
Upper  Trace:  Output  from  Bridge 
Center  Trace:  Bridge  Null  Before  X-Ray  Pulse 
Sweep:  0.2  m sec/division 

Vertical  (Center  and  Upper  Trace) 

Deflection:  O.5  millivolt/division 


Figure  72 


Collector-Base  Junction  Response 
of  a  Type  TO-51  TRW  Transistor 


Lower  Trace: 
Upper  1'race: 
Center  Trace: 
Sweep: 
Vertical 
Deflection: 


X-Ray  Monitor  Output  «  0.6  x  10^r/sec 

Output  from  Bridge 

Bridge  Null  Before  X-Ray  Pulse 

0.2  lisec/division 

(Center  and  Upper  Trace) 

0.5  millivolt/division 
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SECTION  VII 

A  MONTE  CARLO  CODE  FOR 
TRANSIENT  RADIATION  EFFECTS 
By  Win.  J.  Byatt  and  Dale  Sparks 


1.  Introduction 

The  prompt  gamma  pulse  from  a  nuclear  explosion  acts  as  a 
source  for  the  generation  of  Compton  and  photoelectric  electrons 
within  and  around  the  components  of  electronic  circuits.  Tie 
conductivity  of-  both  resistors  and  capacitors  suffers  transient 
changes  in  rated  values  because  of  the  release  of  electrons. 

It  is,  in  part,  our  purpose  to  study  the  time -dependent  changes 
in  conductivity. 

A  Monte  Carlo  code  has  therefore  been  constructed  which  will 

\  ■ . 

give  not  only  the  time-dependenc  absorption  and  scattering  of 
gamma  rays,  but  also  the  number  of  electrons  generated  as  a 
function  of  energy.  The  subsequent  time -dependent  behavior  of 
the  electrons  can  then  be  discussed  by  available  multiple  scat¬ 
tering  theories;  or  crude  estimates  of  the  time -dependent  behavior 
of  the  electrons  can  be  obtained  from  a  knowledge  of  recombination 
times  within  a  given  material.  If  the  time-dependent  behavior  of 
the  number  of  electrons  is  known,  then  the  conductivity  is  given 


by 


ff(t) 


ne(t)e2t(v) 

mv 


(143) 


where  n  (t)  is  the  number  of  free  electrons  per  unit  volume  of 
charge,  and  of  (effective)  mass,  m;  and  velocity,  v;  and  whose 
mean  free  path  is  t(v).  But  equation  (143)  can  also  be  approxi¬ 
mated  by 

a(t)  n#(t)  r  (144) 

where  r  is  a  recombination  time.  If  the  pulse  of  gamma  rays  has 
passed  through  the  material  so  that  there  is  no  longer  any  genera 
tion  of  electrons,  then  it  can  be  crudely  estimated  that  the  time 


163 


dependent  behavior  of  electron!  is  given  by 


»e(t)  -  n0«-t/T  (145) 

with  r  the  recombination  time.  But  by  equation  (144)  this  means 
that  after  the  pulse  has  passed,  rr( t)  should  decay  according  to 

o(t)  -  *0e~t/T  (146) 

if  there  is  single  recombination  time.  If  there  are  several 
competing  processes, 

o( t)  -  co  E  Aie-t/T  (147) 

Again,  it  is  to  be  stressed  that  equations  (146)  and  (147)  will 
hold  only  after  the  gamma-ray  pulse  has  passed  and  there  is  no 
longer  any  generation  of  electrons  by  Compton  scattering  or  photo¬ 
electric  absorption  of  photons.  It  will  be  observed  that  a  rela¬ 
tion  such  as  equation  (147)  has  been  used  by  workers  in  the  field 
of  transient  radiation  effects  to  fit  observed  conductivity -time 
behavior  (Reference  42) . 

The  above  introductory  remarks  serve  to  illustrate  the  impor¬ 
tance  of  having  an  accurate  knowledge  of  the  time -dependent  number 
density  of  free  electrons. 

2.  Summary 

A  Monte  Carlo  code  has  been  constructed  which  will  be  used 
to  calculate  quantities  of  importance  in  transient  radiation 
effects.  It  is  constructed  so  that  the  following  can  be  found: 

a.  The  time  spreading  of  the  prompt  gamma-ray  pulse  from  a 
weapon:  since  transient  radiation  effects  are  dose-rate  dependent, 
the  spreading  is  of  interest. 

b.  The  number  of  freo  electrons  generated  by  the  prompt 
pulse;  since  the  free  electrons  contribute  to  changes  in  conduc¬ 
tivity,  an  accurate  knowledge  of  the  free  electron  density  is  of 
importance . 

In  the  next  few  months,  the  code  will  be  used  to  accomplish 
the  above  program.  In  this  short  report,  the  results  of  two 
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problems  have  been  included:  one  in  plane  geometry,  and  the 
other  in  spherical  geometry,  for  which  some  knowledge  of  the 
solutions  is  available.  Compar ison  shows  that  the  code  is 
working  properly. 

3.  The  Transport  Equation 

A  code  will  be  described  which  has  been  constructed  to 

find  n  (t)  by  solving  the  equation 
6 

—  4.  X  .  V  I  +  9  I  m 

C  ot  *  T 

k  (148) 

S  +  \  k(k'-k)dk'  ^  6(l-k  +  V-  mx)  Id  ft* 

k-2  O 

where  I  is  the  intensity  of  photons  of  speed  c  in  the  wavelength 
range  (k,  k  +  dk)  at  the  space -time  point  (x,t)  traveling  in  a 
direction  whose  cosine  is  denoted  by  u?  i.e.,  I  *  l(x,  t,  |i,  k) . 
The  quantity  t  *  v  1  is  geometry  dependent,  being  given  by 


t  *  v  1  *  M  |I 

in  plane  geometry,  and  by 


(1*9) 


(150) 


in  spherical  geometry.  The  notation  9f  stands  for  the  total  of 
absorption  (<r  )  plus  scattering  (9  )  cross  sections;  that  is, 
9*9+9,  where  each  is  to  be  multiplied  by  an  appropriate 
density  if  9  is  given  in  units  of  cnr/gm  if  x  is  expressed  in 
centimeters.  The  quantity  k(k'  -*  k)  is  the  Klein-Nishina  cross 
section  with  k  expressed  in  units  of  h/mc.  It,  too,  has  dimen¬ 
sions  of  reciprocal  length  if  x  has  dimensions  of  length.  The 
presence  of  the  delta  function  in  the  integral  on  the  right-hand 
side  of  equation  (148)  insures  that  the  Compton  condition  is 
fulfilled  (i.e.,  that  energy  is  conserved  in  a  scattering  event). 
The  source  S  must  have  the  dimensions  of  intensity  per  unit 
length  if  ct  and  x  are  expressed  in  units  of  length. 
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Whether  problems  are  to  be  solved  in  plane  or  in  spherical 
geometry,  one  writes 


1 

c 


V  I 


dl 

Hi 


(151) 


The  introduction  o£  (151)  requires  that  the  independent  variables 
x,  t,  u  (and  hence  X)  are  functions  of  s.  The  parametric  curves 
x(s)t  m(s),  and  so  forth,  are  known  as  the  characteristics.  The 
characteristics,  which  change  direction  on  scattering,  or  terminate 
with  photon  absorption  are  the  paths  along  which  the  radiation 
flows. 

Since 

dtt  dl  at  t  dl  dx  .  31  dn 

ds  dt  ds  dm  ds  ^  ds  (^52) 

one  has,  in  plane  geometry 


as  .  i  ,  a* 

ds  c  ds 
with  the  solutions 


t  (s  ) 


x(s)  -  *(s0)  +  m(s  -  s0) 

and 

M  *  m0,  a  constant 


(155) 


us*) 


In  spherical  geometry,  the  differential  equations  for  the  charac¬ 
teristics  are 


dt 

HI 


I « 


dx 

Hi 


(155) 
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with  solutions 


(a  -  a.) 

t(.)  -  t(.G)  +  --■■-g- 

x2(»)  -  x2(s0)  +  (*  -  »0)2  +  2x(«0)  w(»0)(»  -  »0) 

and 

x(s)  m  (a)  *  x(*0)  +  (s  “  ®o^  (156) 

with  the  introduction  of  the  characteristics,  the  transport 
equation  can  now  be  written 


S  +  o  I 

cTs  UT- 


S  +  TI 


(157) 


where  A 
TI 


*  \  2  "*  X^dA'  \  6(l  -  A  +  A '  -  UX)  Id  fi' 


The  "solutiorf  can  then  be  written  as 


-aVs-s-O  ( 

l(s)  «  l(«0)e  °  *  \  • 


~<*r  (»“••) 


•[s(s')  +  Tl(s')J  ds' 


(158) 


The  Monte  Carlo  code  must  be  constructed  so  as  to  provide  a 
direct  simulation  of  equation  (158). 

4.  Description  of  the  Monte  Carlo  Code 

The  physical  simulation  will  be  described  by  fixing  attention 
on  a  problem  in  plane  geometry.  The  extension  to  spherical  geometry 
will  be  clear. 

Assume  the  situation  pictured  in  Figure  75*  A  source  of 
photons  is  incident  normally  on  the  face  x  ■»  o(s  ■  sQ)  of  a  slab 
of  material  of  thickness  t. 
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Figure  73.  The  Geometry  for  the  Scattering 
and  Absorption  of  Photons  in  a  Slab  of  Thickness  l . 


After  traveling  a  distance  s  into  the  material,  the  photon  will 
either  suffer  a  collision  or  be  absorbed.  The  first  task  is  to 
find  the  point  s  at  which  such  an  event  takes  place.  If  the  wave* 
length  (or  energy)  of  the  photon  is  known,  crT  can  be  found  from 
the  tables  of  G.  W.  Grodstein  (Reference  43).  On  writing 


and  solving  for  (s  -  Sq),  one  has 


(• 


(159) 


(160) 


If  a  random  number  R 


-si 


1 


-jy,  uniform  in  the  interval  o  *  R  *  1, 
is  now  cast,  (a  *  s^)  is  determined.  The  point  at  which  the  photon 
undergoes  its  first  interaction  with  the  material  is  thereby  deter¬ 
mined.  Knowing  (a  -  Sq)  one  can  then,  for  either  plane  or  spheri¬ 
cal  geometry,  calculate  the  time  at  which  the  event  has  occurred 
and  the  position  x(s)  at  which  it  has  occurred.  If  the  position 
is  outside  the  material,  the  energy  is  recorded  and  another  photon 
history  is  begun.  If  the  event  occurred  within  matter,  its  nature 


168 


must  new  be  determined.  It  is  already  assumed  that  one  has  either 
a  Compton  scattering  or  a  photoelectric  absorption.  The  relative 
probability  of  scattering  is  given  by  the  number 


(161) 


while  the  probability  of  absorption  is  1  -  R^,  Another  random 
number  R,  uniform  in  0  s  R  *  1,  can  then  be  cast  such  that  if 
R  >  R^,  an  absorption  has  occurred.  If  the  event  is  an  absorption, 
the  energy  of  the  ejected  electron  is  given  by 


E  -  Eb  (162) 

where  E  is  the  original  photon  energy  and  E^  is  the  binding  energy 
of  the  electron.  The  energy  (162)  is  recorded  as  an  electron  energy. 
The  initial  position  of  the  electron  is  known  and  the  directional 
distribution  as  a  function  of  energy  is  given  in  Davisson  and  Evans 
(Reference  44).  On  sampling  from  the  distribution,  one  has  the 
photoelectron  direction  relative  to  the  direction  of  the  incoming 
photon . 

If  the  event  is  a  Compton  scattering,  the  new  direction  of 
the  photon  and  its  new  wavelength  must  be  calculated.  The  two 
quantities  are  related  by  the  requifement  that 


X  -  X  -  1  (1 63) 

so  that  X  *  X  +  -1  will  be  between  X  -  2 (for  ^  -  -1)  and  X 

(for  no  deflection).  It  is  seen  then  that  if  a  random  number  R 
is  cast,  uniform  in  the  interval  0  *  R  *  1,  and  if 

X 

^  k(X'  -»X)dX' 

r  -  hzl -  (164) 

J  k(X*  -  X)dX' 

X-2 
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one  must  determine  A^,  the  upper  limit  of  the  integral  in  the 
numerator,  so  that  equation  (164)  is  satisfied.  Then  equation 
(l63)  can  be  employed  to  find  |a^.  Since 

1^  =  UU'  +  ,/l-M2  cos  (tp  -  <P')  (165) 


with  m 1  the  cosine  of  the  angle  before  the  scattering,  and  |i  its 
new  cosine,  one  can  solve  for  m •  It  is  necessary,  in  view  of 
azimuthal  symmetry,  to  cast  a  random  number  in  0  s  ?)  -  tp'.  s  2tt 
and  compute  its  cosine  to  find  (-1  uniquely.  Since  the  change  in 
energy  of  the  photon  is  known,  the  energy  given  to  the  (Compton) 
scattered  electron  is  also  known.  With  the  energy  of  the  electron 
known,  its  direction  can  be  found  (Reference  45) . 

With  the  new  energy  and  direction  of  the  photon  given,  the 
new  total  cross  section  can  be  found  and  the  cycle  given  above 
repeated  until  the  photon  either  escapes  from  the  material  or  is 
absorbed . 

5.  Flow  Charts  and  Numerical  Results 

The  Monte  Carlo  Code  has  been  applied  to  a  variety  of  prob¬ 
lems.  Two  are  discussed  here.  Since  the  code  can  be  developed 
in  steps  by  means  of  subroutines,  the  problems  have  been  chosen 
so  as  to  determine  whether  a  given  subroutine,  or  combination 
thereof,  was  operating  properly.  The  coding  was  done  in  Fortran 
language  and  a  separate  program  listing  is,  as  a  consequence,  given 
for  each  problem.  Table  XIV  is  a  listing  for  the  first  problem 
discussed  below.  The  flow  charts  given  in  this  section  are  for 
the  following  situation:  Photons  of  an  initial  energy  of  1  mev 
are  incident  normally  on  the  face  of  a  carbon  slab  of  two  mean 
free  paths'  thickness.  It  is  required  to  find  both  the  unscat¬ 
tered  and  the  scattered  photon  intensities  which  are  transmitted 
through  the  material.  A  printout  was  also  requested  of  the  number 
of  electrons  ejected  within  the  material  as  a  function  of  energy 
and  the  quantity  EN(E),  the  energy  distribution  of  the  electrons. 

No  record  of  the  spatial  distribution  of  electrons  was  called  for, 
nor  was  the  time  dependence  called  for. 
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Table  XIV.  Program  Listing  for  the  Monte  Carlo 
Calculations  of  the  Scattering  of  Photons  in  a  Slab 


*0808 

DIMENSION  LIST  (lOO) 

DIMENSION  EE (101),  N(lOl),  TIG(lQO),  TIGA(lOO) 

COMMON  EE,  N 

23  FORMAT  (//,32HANGULAR  DISTRIBUTION  OF  PHOTONS) 

26  FORMAT  (/,  3X,12HC0SING  THETA,  3X,  14HNO.  OF  PHOTONS) 

27  FORMAT  (1F10.7, 10X.115) 

3  FORMAT  (IX,  l4,  10X,  14  ,  10XF10.2 

4  FORMAT  (4F6.2) 

5  FORMAT  (2l4,3El6.8) 

6  FORMAT  (//f3H  J  ,  10X,  18HNUMBER  OF  ELECTRON,  iaX,6HENERGY) 

7  FORMAT  (214, 14H  PHOTOELECTRIC) 

18  FORMAT  (214,  15HPHOTON  ABSORBED) 

1=0 

24  CONTINUE 

READ  4,  (TIG(K),  TIGA(K),  K=l, 100 
DO  2  J=l, 100 
LIST(j)=0 
EE(J)  =  0.0 
2  N(J)  =  0 

20  READ  4,  E4,  COS  ,  D  ,  CONV 
J  =  E4* 100.0 
LI  *  J 

SIG  *  (TIG(J)  +  TIGA(J))  *  CONV 

SIGA  =  TIGA(J)  *  CONV 

WLE  =  0.50903/E 4 

DO  8  J=l,  100 

WLEN=WLE 

E1=E4 

COSP=COS 

K=0 

DX=0 . 0 

13  CALL  RAND  (l,R) 

1=1+1 

S»-(L0GF(R))/SIG 

DX=DX+COSP*S 

IF  (DX-D)  14,10,10 

14  IF  (DX)  10,10,9 

10  IF  (COSP)  30,30,31 
31  L=(COSP  +  O.O5)  *20.0 
LIST(L)  »  LIST(L)  +  1 

30  IF  (SENSE  switch  3)  25,8 

25  PUNCH  5,  J,K,  COSP,  WLEN 
GO  TO  8 

9  CALL  RAND  (l,R) 
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Table  XIV  (cont'd) 


IP  (R-SIGA/SIG)  11,11,12 

11  PUNCH  7,  JtK 
GO  TO  8 

12  CONTINUE 

CALL  COMP  (WLEN,  ADM,  R) 

K=K+1 

COST- 1 . 0-WLEN  -ADM 
WLEN=ADM 

CALL  SCORE  (  El,  WLEN,M) 

IP  (ElrO.Ol)  17, 17# 16 
17  ip  (sense  switch  3)  19,8 

19  PUNCH  l8,  J,K 
GO  TO  8 

16  sig  =  (tig(m)  +  TIGA(M)  *  CONV 
SIGA  «  (TIGA(M))  *  CONV 
CALL  RAND  (l,Rj 

r=r*6 . 2831055 

COSP=COST*COSP -COSF ( R) *SQRTF ( ( 1 . 0-COST**2) * ( 1 . 0-COSP**2) ) 
IP  (SENSE  SWITCH  l)  15/15 
15  PUNCH  5,  Y  K,  COSP 
GO  TO  13 
8  CONTINUE 
PUNCH  6 
DO  1  J=l,  LI 

1  PUNCH  3  ,  J  ,  N(J)  ,  EE(j) 

PUNCH  23 
PUNCH  26 
DO  28  M=l, 21 
THETA  -  M 

THETA  -  THETA*0.05  -  0.05 
28  PUNCH  27,  THETA, LIST (M) 

IF(  SENSE  SWITCH  9)  21,20 

21  PAUSE 

22  IF  (  SENSE  SWITCH  2)  24,20 
END 


i 

s 


! 


I 


J 
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If  one  has  the  above  in  mind,  the  flew  charts  as  given  in 
Figures  7^  -  77  are  appropriate.  Figure  78  gives  the  angular 
distribution  of  photons  which  emerge  at  the  back  face  of  the 
material  when  500  photons  were  incident.  Note  that  in  region  1 
on  Figure  78  the  number  of  unscattered  photons  is  also  given. 
There  are  65  of  these.  If  x  is  expressed  in  units  of  mean  free 
path,  then 

500  e~2  =  68 


should  be  the  number  of  photons  traversing  the  material  suffering 
neither  scattering  nor  -^sorption.  The  ramainder  of  Figure  78, 
the  angular  distribution  of  scattered  photons,  is  roughly  Gaussian. 
Such  a  result  is  to  be  expected  on  the  basis  of  the  analytical  work 
of  L.  L.  Foldy  or  I.  Ogievetskii  (Reference  ^6).  The  above  two 
checks  on  the  performance  of  the  code  are  comforting.  One  may 
inquire,  finally,  as  to  the  connection  between  Figure  78  and 
equation  (158).  The  latter  is 


1(b) 


-a(a-s') 
e  ds 1 


s 


o 


In  the  problem  under  consideration,  the  volume  source  of  photons 
which  could  be  generated  by  electron  absorption  was  not  calculated. 
Thus,  S  ss  0  in  the  present  case.  The  connection  between  equation 
(158)  and  Figure  78  is  now  clear.  For  (i  >  0,  Figure  78  is  the 
Monte  Carlo  solution  to  (158).  The  photons  in  region  1  are 
equivalent  to  the  term 


while  the  given  distribution  of  scattered  photons  is  the  physical 
simulation  of  the  second  term  on  the  right-hand  side  of  (158). 

It  will  be  observed  that  the  intensity  is  greater  than  one  would 
expect  on  the  basis  of  a  "good  geometry"  experiment;  and  this,  in 
turn,  leads  to  the  concept  of  a  build-up  factor. 
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Read  in  table 
of  a*  a 


DO  1,100 


photon 


last  card 
^een  reac* 


^  PUNCH 
Diatribution  of 
electrona  and 
their  energies 


senae  X 
awitch  2/ 
Non?/^ 

if  START 
AGAIN 


Pause 


'''DO  LoopX 
sat ia fled. 


(From  C) 


/  Has  N 
photon  left  the 
"X^ma  ter  la XyS 


/'if  sense  switch  3 
is  on,  PUNCH  no,  of 
collisions  suffered, 
V  final  coap,  final  \ 


Generate  random  no,  Rl 


Determine  from  R  ^x^ 
whether  the  collision  is  r 
.photoelectric  event  or/ 
X^rcxupton  scattering 


Photo. 


f/UNCH  "PHOTOELECTRIC" 


Compton  ^  (Go  to  b) 

Figure  74.  Flow  Chart  of  Main  Program 
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Continuation  of  Figure  74 


Take  value  of  WLEN  from  main  program 


Take  initial  energy  and  final  wavelength 
of  photon  from  main  program 


Compute  final  energy  of  photon  and  energy 
loss  (imparted  to  ejected  electron) 


Figure  77 •  Randon  Number  Generator  Subroutine 
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NUMBER  Of  PHOTONS  LEANING  MATERIAL 


l 

l 


1 - 1 1 - 1 - 1  |  I  I  I  I  1  l  |  1  !  I 

iO  99  90  49  40  .79  .70  4  9  40  .99  .90  .49  .40  .99  .90  .29 


COSINE  Of  ^ 


Figure  78.  Angular  Distribution  of  Photons 
Traversing  Two  Mean  Free  Paths  of  Carbon 
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The  second  problem  discussed  herein  is  done  in  spherical 
coordinates.  Let  a  point  source  emit  photons  of  1-mev  energy 
radially.  If  scattering  is  absent,  it  develops  that  the  trans¬ 
port  equation  can  be  solved  exactly  for  the  case  at  hand.  The 
analytical  solution  for  the  photon  intensity  is 

I  e-CTr 

I(r,t)  -  °  ,  f(t  -  |  )  H(r)  t  *  |  (166) 

4irr‘  0  ° 

and  l(r,  t)  *  0  for  t  <  — .  The  quantity  f(t  -  £)  is  the  time 

c 

dependence  of  the  source  which  is  arbitrary,  while  H(r)  is  the 
Heaviside  unit  function. 

A  Monte  Carlo  simulation  of  the  above  absorption  problem 
was  done  to  check  the  coding  in  the  spherical  case.  The  results 
are  shown  in  Figure  79*  Included  therein  is  the  intensity  with 
absorption  removed,  and  with  both  absorption  and  spherical 
spreading  removed.  If  in  equation  (l66)  both  sides  are  multi¬ 
plied  by  eCTr,  the  result  is 

eorl  =  f(t  -  |  )  H(r)  (167) 

4tttz 

In  the  absence  of  scattering  f(t  -  £  )  does  not  change  shape. 

Am  v> 

The  spatial  dependence  of  eUTI  should  then  be  proportional  to 
r“^.  The  close  agreement  between  the  slope  of  the  e^I  curve 
(absorption  removed)  and  a  curve  with  a  slope  of  -2  on  log-log 
paper  is  evident  in  Figure  79*  Finally,  if  both  sides  of  equa- 
tion  (167)  are  multiplied  by  47rr  ,  the  right-hand  side  of  (167) 
is  a  constant.  This  fact  is  indicated  in  Figure  79  in  which  it 
is  shown  that  the  curve  labeled  "absorption  and  spreading  removed" 
is  a  constant. 

From  the  above,  it  can  be  inferred  that  problems  in  spherical 
geometry  can  be  successfully  undertaken:  i.e.,  the  coding  is  such 
that  the  Monte  Carlo  simulation  is  correctly  done. 

Parenthetically,  it  is  added  that  the  first  of  the  two  prob¬ 
lems  discussed  herein  was  done  on  the  UNM  IBM  1620  computer  for 
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NUMBER  OF  PHOTONS  (SCALE 


Figure  79.  4,000  Photon  Histories  -  Pure  Absorption 


181 


■jl 


500  photon  histories  and  consumed  about  1  1/2  hours  of  computing 
time.  Including  "debugging"  on  the  Kirtland  AFB  CDC  1604,  the 
second  problem  discussed  took  about  18  minutes  for  4,000  photon 
histories. 
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SECTION  VIII 


GENERAL  CONCLUSIONS  AND  RECOMMENDATIONS 
By  W.  W.  Grannemann 

The  experiments  on  Carbon  composition  resistors  and  the  air 
effects  surrounding  them  show  that  the  radiation  effects  on  this 
combination  increase  until  a  pressure  of  about  150  mm  Hg  is 
reached.  Then  the  radiation  effects  decrease  until  a  pressure 
of  1  mm  of  Hg  is  reached;  the  magnitude  then  remains  essentially 
constant  with  decreased  pressure.  In  the  pressure  range  10  to 
100  mm  Hg  the  theory  indicates  that  the  time  dependent  behavior 
of  the  number  of  free  electrons  plays  the  dominant  role  in  the 
air  conductivity  around  the  resistor.  The  value  of  resistance 
changes  as  a  function  of  dose  rater  and  the  percentage  change  of 
resistance  is  related  to  the  original  value  of  the  resistance. 

The  major  transient  resistance  changes  under  gamma  radiation  are 
due  mostly  to  air  ionization  effects.  This  further  emphasizes  the 
importance  of  potting  resistors  and  any  open  circuit  wires,  as 
well  as  the  importance  of  designing  circuits  with  as  low  an 
impedance  as  possible  for  the  job  they  are  designed  to  do. 

The  thin  film  Hall  effect  devices  were  found  to  be  highly 
radiation  resistant  to  transient  gamma  pulses.  The  theory 
developed  indicates  that  the  type  of  circuit  used  for  the  Hall 
effect  device  can  make  considerable  difference  in  the  radiation 
tolerance  of  it.  Experimentally,  no  measurable  change  in  Hall 
voltage  was  observed  for  dose  rates  of  10'  r/sec.  Prom  theoreti¬ 
cal  calculation  it  appears  that  the  thin  film  Hall  effect  devices 
will  operate  well  in  radiation  dose  rates  of  30^  to  1010  with  failure 
occurring  somewhere  near  101*  r/sec.  It  was  also  found  that  dis¬ 
placement  voltage  effects  under  irradiation  can  appear  to  add  to 
the  transient  Hall  effect  voltage  so  that  for  a  radiation  resist¬ 
ant  Hall  effect  device  the  alignment  of  the  Hall  probe  is  extremely 
critical.  It  appears  that  thin  film  Hall  effect  devices  can  be 

recommended  for  use  in  a  transient  gamma  radiation  environment 
o 

up  to  107  r/sec  and  with  special  circuit  configurations  this  value 
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might  be  raised  one  or  two  orders  of  magnitude,  but  it  is  recom¬ 
mended  that  the  Hall  effect  devices  be  checked  on  higher  dose 
rate  machines  to  verify  the  theoretical  predictions. 

Titanium  dioxide  diodes  have  been  found  to  be  highly  radia¬ 
tion  resistant  both  to  neutrons  and  to  gamma  pulses.  No  transient 
radiation  effects  were  observed  for  gamma  pulses  of  10'  r/sec. 

Prom  experimental  data  it  appears  that  the  titanium  dioxide  diodes 
would  go  ccnsiderably  beyond  the  two  orders  of  magnitude  improve¬ 
ment  over  the  normal  junction  diode.  It  is  recommended  that 
titanium  dioxide  diodes  also  be  investigated  further  at  higher 
radiation  intensities  to  determine  accurately  the  magnitude  and 
shape  of  the  transient  current  pulse  resulting  from  the  radiation 
pulse;  and  further  development  and  improvement  of  electrical 
characteristics  of  titanium  dioxide  diodas  is  recommended.  These 
thin  film  devices  using  TiC^  diodes  hold  considerable  promise 
for  radiation  resistant  diodes. 

The  absorption  method  for  obtaining  X-ray  spectrums  has 
been  satisfactorily  used  on  the  600  KV  flash  X-ray  machine. 

The  present  method  requires  knowing  the  efficiency  of  the  energy 
absorbed  in  the  detector  and  the  reproducibility  of  the  X-ray 
pulse  from  one  pulse  to  the  next.  The  accuracy  of  the  computa¬ 
tion  is  dependent  on  the  For  this  reason  the  method  may 

require  some  modification  to  be  useful  on  the  super  flash  X-ray 
machines.  It  is  recommended  that  this  method  of  spectrum  deter¬ 
mination  be  considered  when  the  need  arises  for  spectrum  measure¬ 
ments  on  super  flash  X-ray,  even  though  it  may  be  desirable  to 
obtain  a  complete  spectrum  from  one  individual  pulse. 

Measurements  have  been  made  on  dielectric  materials  showing 
it  is  now  possible  to  make  Hall  effect  and  magnetoresistance 
measurements  on  some  of  the  usual  dielectric  materials  such  as 
mylar,  polyethylene,  and  similar  dielectric  materials  used  in 
capacitors.  This  makes  possible  a  measurement  of  mobilities  of 
the  charge  carriers  which  are  predominantly  holes  in  these  cases. 
It  is  also  possible  to  calculate  the  resistivity  as  a  function 
of  dose  rate  from  these  measurements .  The  largest  bulk  effect 
was  the  radiation  induced  conductivity  change. 
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A  Schering  Bridge  was  used  to  measure  the  effective  capacity 
change  of  diode  junctions  under  irradiation.  The  bridge  arrange¬ 
ment  gives  a  transient  picture  of  the  effect  of  charge  redistri¬ 
bution  in  the  diode.  It  is  planned  to  use  this  method  of 
measurement  on  transistor  elements  and  diode  elements  in  micro- 
circuits  in  the  coming  project. 

A  Monte-Carlo  code  for  transient  radiation  effects  has  been 
developed.  The  code  predicts  the  number  of  electrons  produced 
as  a  function  of  energy  for  a  transient  gamma  radiation  pulse 
incident  on  a  material.  Included  in  this  code  is  information 
on  the  time  dependent  absorption  and  scattering  of  the  gamma 
rays.  Once  the  number  of  electrons  produced  in  the  material 
is  known  as  a  function  of  time,  the  change  in  conductivity  as 
a  function  of  time  can  then  be  calculated.  The  code  can  also 
be  used  to  determine  gamma  dose  after  the  gamma  rays  have  passed 
through  some  intermediate  material  before  striking  an  electronic 
device. 
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APPENDIX  I 

AFWL  600  KV  PLASH  X-RAY  FACILITY 
By  LeRoy  Meyer  and  W.  W.  Grannemann 

i.  AFWL  600  KV  Flash  X-Ray  Facility 

The  AFWL  flash  X-ray  system  consists  of  a  600  KV  Fexitron 
unit  manufactured  by  Field  Emission  Corporation.  The  charac¬ 
teristics  of  the  unit  are  summarized  in  Table  XV.  A  double 
walled  rf  shield  room  houses  the  experimental  measuring  equip¬ 
ment.  By  this  means  the  rf  noise  level  is  reduced  by  120  db. 

The  X-ray  beam  is  fired  through  a  special  window  in  the  shielded 
room.  This  window  is  shown  in  Figure  80.  The  beam  area  and 
dose  rates  are  shown  in  Figure  8l.  The  room  layout  is  shown  in 
Figure  82.  The  energy  spectrum  of  the  system  is  shown  in  Figure 
85  for  400  KV  X-rays.  A  graph  showing  the  conversion  of  gamma 
flux  equivalent  to  1  roentgen-hour  as  a  function  of  gamma  energy 
is  shown  in  Figure  84. 


Table  XV.  Operating  Characteristics  of  the 
600  KV  Fexitron  Flash  X-Ray  System _ 


Output  voltage 

600  KV 

Output  current  at  peak 

output  voltage 

2000  amp 

Pulse  width 

0.2  fi  sec 

Pulse  rise  time 

50  nanosec 

Dose  rate  maximum 

2  x  107  r/sec 

i 

? 

i 
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APPENDIX  II 
TEST  EQUIPMENT 

By  LeRoy  Meyer  and  W.  w.  Grannemann 

The  test  equipment  that  was  used  in  the  experiments  is 
described  below.  The  major  items  are: 

1.  The  vacuum  system 

2.  DC  micro-volt-ammeter 

3.  Amplifier 

4.  NJE  High  voltage  power  supply 

3.  Harrison  labs  power  supply,  model  80$A 

6.  Harrison  labs  power  supply,  model  855B 

7.  Atomic  Laboratory  Electromagnet 

1.  Vacuum  System 

A  special  vacuum  system.  No.  1308,  built  by  Consolidated 

Vacuum  Corporation  of  Palo  Alto,  California,  is  capable  of 

—A 

reaching  the  10  torr  range  using  liquid  nitrogen  coolant  in 
the  baffle.  The  controls  are  pneumatically  operated  from  the 
control  panel  located  outside  the  lead  shielding  enclosure. 
Heating  capabilities  and  a  variable  leak  valve  are  also  built 
into  the  system. 

Two  vacuum  gauges  are  used.  They  are  the  magnevac-type 

GMA-140  and  a  CVC  ionization  gauge  type  GIC-110A.  The  magnevac 

gauge  has  a  range  of  from  1  micron  to  500  mm  and  the  ioniza- 

-3  _Q 

tion  gauge  has  a  range  of  from  1  x  10  '  torr  to  2  x  10  ^  torr. 

The  bell  jar  guard  has  a  window  for  the  X-ray  beam  to  pass 
through.  A  covering  of  aluminum  foil  is  also  used  over  the  bell 
jar  for  electrostatic  shielding. 

2 .  DC  Micro-Volt-Ammeter  Model  MV-07C 

The  model  MV-07C  is  an  extremely  stable,  sensitive,  wide 
range  multiplier  and  amplifier  with  full-scale  range  as  low  as 
10  microvolts  and  10  micro-micro-anperes  and  as  high  as  1  KV 
and  1  ma. 

Available  from  Millivac  Instruments  Inc.,  110  Altamont  Ave., 
Schenectady,  New  York. 
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3.  Amplifier 


The  1121  Tektronix  amplifier  has  an  overall  gain  of  100 
with  a  maximum  input  of  +  10  mv.  The  input  signal  can  be  attenu 
ated  from  1  x  to  500  x  in  9  calibrated  steps.  The  band  pass  of 
the  amplifier  is  5  cps  to  17  me.  The  noise  level  is  «  50  Mvolts 
4.  NJE  High  Voltage  Power  Supply 
Model  S-327 

500-5000  VDC,  0-10  ma. 

5«  Harrison  Labs  Power  Supply,  Model  809A 

The  809A  is  a  constant  voltage/constant  current,  transistor 
power  supply  which  has  an  output  of  0  to  36  volts  and  0  to  10 
amps.  It  is  equipped  with  a  continuously  adjustable  current 
limit  control  which  allows  the  maximum  output  current  to  be  set 
at  any  value  up  to  the  maximum  current  rating. 

6.  Harrison  Labs  Power  Supply,  Model  855B 

The  855B  is  a  constant  voltage/constant  current  transistor 
power  supply  which  has  an  output  of  0  to  18  volts  and  0  to  1.5 
amps.  It  has  a  continuously  adjustable  current  limit  control 
with  less  than  30  millivolts  drift  for  8  hours. 

7.  Atomic  Laboratory  Electromagnet 

Cat.  No.  79641 

Available  Atomic  Labs.  Inc. 

3086  Claremont  Ave. 

Berkely  5,  California 

This  electromagnet  can  produce  a  variable  magnetic  field 
from  0  to  3  kilogauss  in  a  1/2"  air  gap.  It  also  has  available 
hollow  pole  pieces  which  were  used  in  the  radiation  experiments. 
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APPENDIX  III 


THEORY  OF  THE  SOBERING  BRIDGE 
By  L.  T.  Boatwright 

Among  the  numerous  bridge  network  schemes  available,  the 
one  suggested  by  H.  Schering  is  one  of  the  better  methods  for 
evaluating  small  capacitors  at  low  and  high  voltages  (Reference 
48) .  The  effective  capacitance  and  the  equivalent  series  loss 
resistance  of  an  imperfect  capacitor  can  be  measured  with  high 
precision.  Furthermore,  the  sensitivity  of  the  Schering  net¬ 
work  is  much  greater  for  changes  in  capacitance  than  for  changes 
in  effective  resistance  as  shown  below.  In  transient  radiation 
studies,  instrumentation  techniques  are  required  for  recovery  of 
low  level  signals  in  the  presence  of  considerable  interference. 

It  appears  that  an  alternating  current  bridge  operating  at  fre¬ 
quencies  above  10  megacycles  may  offer  some  advantages  over  the 
usual  method  of  viewing  the  phenomena  directly  on  an  oscilloscope, 
The  Schering  Bridge  Network  has  been  described  extensively 
in  the  literature  (References  48-53);  however,  bridge  utiliza¬ 
tion  for  the  observation  of  transient  radiation  phenomena  is 
comparatively  recent  (Reference  5*0 •  The  basic  Schering  Bridge 
circuit  is  shown  in  Figure  85*  and  R^  are  the  effective  capa¬ 
citance  and  series  loss  resistance  respectively  of  the  capacitor 
to  be  tested;  is  a  loss -free  standard  capacitor;  R2  and  R^  are 
non-reactive  resistances,  and  is  a  variable  air  capacitor.  In 
addition  to  0^,  other  components  in  the  bridge  arms  may  be  vari¬ 
able  in  order  to  obtain  balance.  Bridge  balance  occurs  when 
there  is  no  detector  current,  or  the  potential  at  C  and  D  is 
the  same.  In  practice,  exact  balance  may  not  be  possible  and 
a  null  is  obtained  for  minimum  current  through  the  detector. 

Z^.  and  Zg  are  the  detector  and  signal  source  impedances  respec¬ 
tively.  The  analysis  below  assumes  sine  wave  excitation  of  the 
Lett 
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bridge . 
a. 


Ri + 


135c, 


U60) 

(169) 
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r  - 


(170) 

(171) 


Detector 


Figure  85*  Schering  Bridge 

Loop  equations  are  written  for  the  circuit  of  Figure  85  and  are 
solved  for  the  detector  current: 


(21  +  *4  +  z5^Ii  “  Z5I2  “  z4Ij 


(172) 


51!  +  (*2  +  *3  +  *5^2  ' 


-z.I 


(173) 


41!  -  zjx2  +  (zj  +  *4  +  z$) 


-zi.l 


(174) 


If  equations  ( 172) (173)  and  (174)  are  solved  utilizing  deter¬ 
minants, 

1 0  -  Z-  I 


0  *2  +  *3  +  *5 


*3  +  *4  +  *6 


J  .  E[yS  *  «4(,2  *  *3  +  ■»>] 


(175) 
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and  equation  (177)  for  the  detector  current  becomes 


f  ^  -  «2»4  -i 

L  zJiJCaJ  +  +  *R)J 


E 


(180) 


or 


[a>(cNc3R3  ~  ci2riR2^  +  ~  c1R2^]  ^ 

^2  "**  ^5  ^  *  j  ) 


Id  ~  I  °* 


(18X) 


(182) 
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If  transient  radiation  of  a  capacitor/  for  example/  causes  an 
increase  in  the  detector  current '(normally  zero  at  balance)/ 
this  can  be  ascribed  to  a  change  in  capacitance/  conductance/ 
or  dielectric  constant/  of  which  the  latter  can  be  related  to 
the  loss  angle  9  (9  +  tan-1  coR^c^) .  Equations  (l8l)  or  (l82) 
reveal  that  the  detector  current  is  directly  proportional  to 
both  the  frequency  and  the  applied  voltage.  If  these  are  held 
constant  the  rate  of  change  of  current  with  respect  to  any 
parameter  "p"  is 
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aiD 

~3p 


da 

JE. 


-  a 


dA 

-ftj?  ojE 
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When  the  bridge  is  balanced/  z^z^  *  Z2Z4  or  ciR2  "  c4R3i 
R^c^  m  CjR2  and  tan  9 
balance  a  -»  0  and 


coc^R^  *  coCjRj.  Consequently/  near 


3Id 

~3p 


[sip] 


OiE 
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If  all  other  factors  are  considered  constant  except  those  of 
the  test  item  (c^,  R^)  and  p  ■  c^  and  p  ■  R^  then 
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In  - 


Z1  +  z4  +  z5 

-  z_ 


-  z. 


0 

0 


-  z, 


-  z. 


z5  +  z4  +  z6 


-*[-»3(*l  +  *4  +  z=)  -  z4zg] 


(176) 


The  detector  current  I  =  I-  -  I, 

o  z  1 


Jo  * 


(177) 


where: 


*l  +  *4  +  *5 
*5 
-  ... 


"  *5 

z2  +  zj  +  *5 


-z. 


If  equation  (178)  is  expanded 


-  z. 


z5  +  z4  +  z6 


(178) 


*  -  *5*6  fzl  +  *2  +  *3  +  *4]  +  *5[(*1  +  *2)(*3  +  *4>] 


+  *6[(*i  +  *4)(*2  +  *3}]  +  t*l*2(*3  +  *4)+  *3*4(ll  +  *2)] 


If  the  bridge  ie  to  be  arranged  for  greatest  sensitivity  to  a 
given  out -of -balance  condition,  z^  and  z^  are  required  to  be 
much  larger  than  s 2  and  z^  (Reference  52).  If  this  is  the  case, 
4  can  be  simplified  to: 


A  **  *i“4(*2  +  *5  +  z^) 


(179) 
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APPENDIX  IV 
CABLE  STUDY 

By  Arthur  Golubiewski  and  W.  W.  Srannemann  I 

Many  studies  of  the  effects  of  nuclear  radiation  on  coaxial 
cables  have  been  carried  out  (References  55-63*  for  example) . 

Nearly  all  transient  radiation  studies  have  been  performed  with  1 

I 

pulsed  reactors  on  linear  accelerators  on  cables  whose  electrical  j 

length  is  short  compared  to  duration  of  the  radiation  pulse.  j 

The  puxpose  of  this  study  is  to  determine  the  effects  of  a  ! 

flash  X-ray  pulse  on  a  cable  whose  electrical  length  is  long  com¬ 
pared  to  the  duration  of  the  pulse.  A  mathematical  model  is 
proposed  which  relates  the  expected  output  pulse  to  gamma  induced 
photoconductivity.  It  is  assumed  that  the  cable  has  received  a 
sufficient  radiation  dose  to  eliminate  initial  effects.  j 

It  is  assumed  that  before  irradiation  the  cable  can  be  con-  ] 

i 

side red  "lossless".  The  incremental  equivalent  circuit  is  then  I 

an  L-C  circuit: 


L 


Figure  86 


Upon  irradiation,  gamma  photoconductivity  alters  the  incre¬ 
mental  equivalent  circuit  so  that  there  is  a  conductance  in 
parallel  with  the  capacitance: 


Figure  87 


Assuming  that  the  radiation  pulse  is  square,  one  can  take 
the  conductance,  G,  to  be  sero  until  the  gamma  rays  strike  the 
cable,  whereupon  it  increases  to  a  constant  value  for  the  dura¬ 
tion  of  the  pulse,  decreasing  again  to  sero  at  the  end  of  the 
pulse. 
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Mathematical  Analysis 

By  taking  the  Laplace  transform  of  the  general  transmission 
line  equations 


a2i(x,») 

ax2 


-  (Ls  +  R)(Cs  +  G)  l(x#s) 

c  aedco)  -  l(cs  +  g)[i  (x,o)1 
55c 

(La  +  R) (Cs  +  G)  e(x,s) 

aT(x,o) 

L  — igjj - C(L8  +  R)  E(x,0) 


(190) 


(191) 


where 

x  »  distance  along  the  line 

s  *  the  transform  variable  with  respect  to  time 
E  ■  voltage  across  the  line 
I  «  current 

L  -  inductance  per  unit  length 

C  *  capacitance  per  unit  length 
R  *  resistance  per  unit  length 

G  ■  leakage  conductance  per  unit  length 


By  using  the  circuit  below,  the  following  initial  conditions 
can  be  established. 


Figure  88 
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l(x,0) 
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E(x,0) 


(195) 


Boundary  condition  is: 
E(0,s)  “  — 


(196) 


if  »  Zo  is  used  to  simplify  the  analysis.  In  addition,  to 
simplify  the  analysis  the  resistance,  R,  in  the  incremental 
equivalent  circuit  should  be  taken  to  be  zero.  The  propagation 
constant,  y,  simplifies  to  JUS"  (Cs  +  &},  and  the  characteristic 
impedance  to 


Jr 


>8 

Cs‘+  <5 


Ls 

y 


Inserting  the  initial  conditions,  simplifying,  and  solving 
the  general  line  equations  gives 


a2E  _(**»I  -  LS(Cs  +  G)E(x,  s)  »  -CLsE 


dx4 


(197) 


dx 


-  LS (Cs  +'G)l(x,s)  a  -L(cs  +  G)  It  (198) 


CE 

E(x,  s)  -  ALeyx  +  B1e"yx  + 


(199) 


l(x,s) 


Aje^  +  B2e“yx  +  -~ 


Using  the  first  boundary  condition  in  equation  ( 199 ) 


E(o,s) 


E 


CE 


-t  -  Ai  +  Bi  +  urns 


( 200) 


(201) 


Eo  (  ¥  ■  5sPT  8  ) 


A1  +  ®1 


(202) 
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Taking  a  basic  line  relation  and  sinplifying  yields 


l[si  (x,s)  -  tj  -  *  iSg.feii,3.)  (203) 

Applying  equation  (203)  to  equations  ( 199)  and  (200)  yields 

LS(A2eyx  +  B2e~yx)  +  Ll±  -  =  -A^e**  +  B^e"^ 

^  1  (204) 

Equating  coefficients  of  like  exponents  one  gets 


LsA2  = 

-y\. 

(205) 

LsB^  m 

yBi 

(206) 

a2  - 

“A1  zs 

(207) 

B2  “ 

Bi  fe 

(208) 

If  it  is  remembered  that  Zol(£,s)  «  E(<6, s)#  equations  (199)  and 

(200)  can  be  combined  to  yield: 

I  CE 

iS.(A2eyt+  B2<3-yl  +  fli)  -  A^4  +  Bf-''1  +  (209) 


Multiplying  through  and  substituting  equations  (207)  and  (208) 
one  obtains 


-  Ajeyt  +  ^ 


+  V' l*  CTTHJ  (2«>) 


By  solving  for  A^ 


CE 

o 

-  Cs^nr 


2A. 


f=o _ 

"  cmr 


-yt 
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And  substituting  (21l)  in  (202) ,  and  solving  for 

LI,  CE 

_ . 

y 


B. 


E  (i  -  - — 
o's  Cs  + 


x  r  LIi  CE0  1  yt 
-  7.  L  -y-  -  cnnrJ  * 


(212) 


Rearranging,  one  gets 


B. 


B  = 


C 

Cs  T  <5 


L  e~Yl  ,  CE~yt  1 

7  Try-  +  2(Cs"  +  "c5) '  "  J 


(213) 


Substituting  the  expressions  for  A^  and  B^  in  equation  (199)  one 
obtains  the  equation  for  voltage  in  the  cable  during  irradiation: 


E 


_  1  r  i  o  "L-y(-t-x),  b,  f^l  C  \^-yt 

(x,s)  -  7  L  —  -  Cs  T  8je  +  EoLls  -  Cs  4  S>e 

(21*) 


*  7  (CsC'+' S  -  77 -y)e'y^+*^]  +  — -n 

L  8  +  : 


c 
c 

Since  the  signal  fed  to  the  oscilloscope  is  taken  off  the  load 
resistor,  ZL,  one  can  fix  x  *  l  and  simplify  equation  (214)  some¬ 
what: 


LE 


CE 


E(t's)  =  7Z7y  -  7T2!s  T-<3) 


+  E 


[(?  -  VTTs)e~yl 


*  7  (e8,C,4-  a  -  7^y)a'2rt]  +  5TT! 


(215) 


To  facilitate  transformation  of  equation  (215)  into  the  time  domain, 
E(1,b)  is  broken  up  into  two  parts 


E(t,») 

-  E1(-t,s)  +  E2(^,s) 

(216) 

where 

LE  CE 

Bx(*.») 

“  7Zj  +  g^es-T-s) 

(21?) 

LE  CE_ 

O  _i_  o 

_ _ _  +  . . . 

2Zj J  IBJCt  ■45)'  2<-Cm  + 
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(218) 


E 

o 

r* 


[ 


ITT  1 

^  *1 


1 _ 

Vi~V8  +  G/c 


8  + 


ir  ] 


■J  E1  9ivea  the  voltage  at  the  load  from  the  inception  of  irradia- 

"I  tion  until  the  reflected  pulse  given  by  E£  arrives  at  the  load. 

When  the  reflected  signal  arrives  at  the  load,  the  output  is  the 
sum  of  E^  and  When  the  radiation  ceases,  both  relations  no 

longer  apply.  By  consulting  tables  (Reference  64)  to  .find  the 
inverse  transform 


■*  1 


•"  7  *  i  (*)t 


yi"7»TT  °‘7 

and  letting  a  in  equation  (219)  *  ^  one  obtains 

©1(^t)  *  _£  [  Jl/c  1  -  t 


(219) 


K‘ 


Gt 

Vfff  t)  +  e  ?  3  (220) 


Returning  to  E^C-t,  s) 


V4'-)  -  eo[<f  - 

(221) 

By  multiplying  through  and  expressing  y  in  terms  of  circuit 
constants 


E2(t,s)  - 


-  JTS\Cm  +  B)  l  -JZiXV*  f  a)  l 

e  e 


s  + 


-2^Ls(ts  +  GJ  l 
e  L  tt 


2(»  +  2) 


7  *L  vntCS  +  8) 


(222) 
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The  exponential  factors  have  the  general  form  (Reference  64) 


-(-)  ,/s2  +  2sp  +  p2  -  cr* 
e 

which  has  the  inverse  transform: 


(223) 


e-(p*/v)  #(t  .  £)  +  g  (224) 

where 

* 

p  *  since  R  *  0 

0  "  7E  +  lr  +  -  2?  *inc*  R  ■  0 
*  »  A7  -  ($)*  t>s 

I1(orz)  -  -i  (ioz) 
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e  (£, t)  can  ba  represented  in  terms  of  convolution  integrals 
which  cannot/  in  turn,  be  expressed  in  closed  form.  One  then 
has,  finally, 

-  -  0  l 

e2(t,t)  -  Eo  j  u(t-r)  [•  C7V  «(r  -  4Z 


at  '  7C  T  Ti(-  fg)tT2  -  <'2lc) 

W  (T2  -  t2  LC)^2 


i e  b 


§7 


«(t  -  t^£c  )  - 


gt 

gi  ~  77 

7c  e 


| g  )(T2  -  ^LC)172 

(r2  -  t2  uc)1^2 


]  dT  +  7  l  8 


(t  -  t) 


*  r  V®  at  .  -  7C  T  M"  §iy)(T2-4t2LC)1/2 


3(r-2t,/I C)  -  gi  e 


(r2-4t2Lc)1//2 


5  - 

-  S  B  TJ5— ]  *  • 


-  §(t-T) 


^.(^vsr-gi  .-^T 


Ij,(-  7ff)(f2  -  W^lc)1^2 


(r*  -  4</lc) 


2Tjt»\  1/2 
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It  should  be  noted  that  the  convolution  integrals  in 
equation  (275)  <3°  not  have  to  be  solved  if  the  electrical 
length  of  the  cable  is  longer  than  the  gamma  ray  pulse,  since 
the  reflected  signals  will  no  longer  be  given  by  equation  (275)* 

For  a  pulse  width  of  0.2  Msec  and  a  signal  velocity  of  658 
c 

x  10  ft/sec,  a  cable  at  least  111  feet  long  should  be  used  to 
avoid  *2 (^t). 

Experimental  Setup 

Seventy -two  and  1/2.  feet  of  RG58  C/U  coaxial  cable  was 
spiral  wound  on  a  plywood  board.  This  was  hung  so  that  the 
plane  of  the  board  was  perpendicular  to  the  axis  of  the  beam 
from  the  flash  X-ray  tube  and  was  positioned  to  achieve  uniform 
radiation.  The  cable  was  terminated  with  a  51  ohm  resistor  in 
series  with  a  1  microfarad  blocking  capacitor.  Bias  was  applied 
to  the  cable  by  means  of  a  battery  connected  to  the  other  end  of 
the  cable.  The  signal  was  amplified  by  a  Tektronix  1121  A  ampli¬ 
fier  and  then  displayed  on  a  Tektronix  55*  dual  beam  oscilloscope. 
The  trace  was  photographed  by  a  Polaroid  camera  mounted  on  the 
oscilloscope. 

Results 

pecause  of  the  high  level  of  pulser  generated  r-f  noise 
pickup  in  the  cable,  gamma  induced  effer  s  could  not  be  deter¬ 
mined.  With  the  pulser  and  X-ray  tube  in  a  screen  room,  it  is 
felt  certain  that  r-f  noise  will  be  reduced  to  a  level  far  below 
gamma- induced  signals  in  the  cable. 
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APPENDIX  V 

DIRECTIONAL  RADIATION  STUDY 

By  Goebel  Davis  and  W.  w.  Grannemann 

■  \  • 

1.0  introduction 

Recent  interest  has  developed  in  orientational  effects 
associated  with  the  transient  response  of  electronic  components 
when  exposed  to  gamma  or  X-rays. 

These  effects  are  concerned  with: 

i)  the  response  for  omnidirectional  versus  mono- 

;  r 

directional  radiation*  and 
b)  the  response  for  different  orientations  of  com¬ 
ponents  relative  to  the  incident  radiation. 
Preliminary,  crude  calculations  using  deposited  energy  as  the 
response  criteria  indicate  that  these  effects  can  be  a  factor 
of  significance. 

2.  Analysis 

For  this  analysis,  a  rectangle  of  material  (Figure  89) 

with  dimensions  a,  b,  and  c  respectively  in  the  x,  y,  and  z 

direction*  is  exposed,  respectively,  to  X-ray  radiations  in 

the  x,  y,  and  z  directions.  The  intensities  of  these  incident 

radiations  are  I  (0),  I  (0),  and  I  (0)  respectively,  and  are 

x  y  z  2 

assumed  to  be  all  of  the  same  value,  IQ  (energy  per  cm*  per 
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s 


If  one  approximates  the  spectrum  of  the  X-ray  source  by 


a  single  frequency 
the  emerging  X-rays 

source,  then 

are: 

the  approximate  intensities  of 

Ix(a)  - 

lx(0)e^a  = 

I  e^a 

Ioe  ' 

- 

Iy(0)«-Mb  - 

IQe”fib  ,  and 

* 

Iz(0)e_,Je  - 

I  e"^c 

Ioe  ' 

where 

2 

M/p  *  mass  absorption  coefficient  (cm  /g),  and 
p  *  density  of  the  material  (g/cm^) 
respectively. 

The  energies  deposited  are  proportional  to  the  respective 
frontal  areas  and  the  change  in  intensities  of  the  incident- 
radiations  in  passing  through  the  material.  Let  these  deposited 
energies  be  Ex,  Ey,  and  Ez  respectively.  Then,  approximately, 

Bx  ~  Iobc(1  ' 

Ey  ~  I0ac(l  -  e"ub),  and 

*,  ~  l0«b(l  -  «'uc). 

From  the  form  of  the  above  equations,  it  can  be  seen  that  usually 
E  ,  E  and  E  would  not  be  of  the  same  value.  This  implies  a 

**  y  * 

different  response  for  different  orientations  of  the  components. 
Also,  these  equations  show  that  the  response  to  tridirectional 
radiation  of  intensity  IQ  (in  each  direction)  is  different  from 
the  response  to  monodirectional  radiation  of  intensity  3IQ. 

Hence,  the  response  to  omnidirectional  radiation  is  different 
from  that  for  monodirectional  radiation,  even  though  the  total 
incident  radiation  levels  are  the  same. 
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5.  Numerical  Example 

As  a  sample  numerical  calculation,  it  is  assumed  that  a  *  10 
b  «  10  ,  c  *  10~*  cm,  and  that  the  X-ray  energy  is  all  at  0.5 

mev,  and  the  material  is  silicon.  Then  p/p  «  0.107  cm  /g,  as 
p  *  2.42  g/cm^#  p  *  0.107  cmVg  x  2.42  g/cm^  *  0.259  cm”1 
(Reference  29)*  If  these  values  are  used, 

E  ~  2.60  x  10”T 

E  ~ 
v 

E 

z 

where  the  units  on  Ex#  Ey  and  Ez  are  arbitrary.  It  should  be 
noted  that  the  relative  values  of  these  energies  are  consistent 
with  the  cpnclusions  of  the  previous  section. 

4.  Experimental  Test 

One  experimental  test  of  the  response  of  a  1N491  diode  was 
made  using  the  TREES  flash  X-ray.  During  the  test  the  transient 
response  output  of  the  diode  in  a  balanced  bridge  arrangement 
was  displayed  on  an  oscilloscope.  The  diode  was  first  exposed 
to  radiation  straight  to  the  side  (exposing  the  maximum  area) 
and  then  progressively  exposed  at  tilt  angles  of  50°,  60°,  and 
90°.  The  peak  of  the  transient  output  tended  to  follow  the 
relation  A  +  B  cos  0,  where  A  and  B  are  constants  and  9  is  the 
angle  in  tilt. 

It  is  not  known  at  present  whether  part  of  the  response  is 
due  to  stray  ionization  around  the  diode  case  or  due  to  the 
diode  itself. 

5*  Limitations  and  Future  Studies 

A  more  sophisticated  future  study  should  include t 

a)  The  effect  of  energy  not  dissipated  locally  within 
the  material. 


1.60  x  10”7,  an<j 
2.59  x  10"7 


:V*  *j*  ,  f.  If, 


b)  The  use  of  a  general  X-ray  energy  spectrum. 

c)  The  derivation  of  energy  deposited  as  a 
function  of  a  general  angle  of  incidence  for 
the  impinging  radiation. 

d)  Possible  configurations  other  than  just  a 
rectangle. 

e)  Additional  experimental  verification. 
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